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ABSTRACT 


A pilot plant scale double effect evaporator 
was operated under direct digital control using an IBM 
1800 Data Acquisition and Control System. 

Feed blending was accomplished by the use of 
flow ratio control. Averaging and cascade control were 
applied to liquid holdups in both effects. 

Product concentration was controlled by a feed- 
back control system which manipulated the steam flow rate. 
Feedforward control was added to compensate for feed flow 
and concentration disturbances. In the case of feed flow 
disturbances, the compensation was based optionally on 
the feed flow or first effect bottoms flow signals. Only 
standard direct digital control software was employed 
to implement the feedforward controller, such that the 
method used is expected to be generally applicable to simi- 
lar processes. 

An inferential control system was developed and 
applied that was based on calculated values of product 
concentration rather than directly measured values. 

A series of 22 experiments were conducted to 
evaluate the performance of the control systems developed. 
The use of feedforward control alone gave rise to small 
off-sets in product quality that were eliminated by the 
addition of feedback correction. The combined controller 


maintained essentially constant product quality for changes 
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in feed flow and concentration 20% of their steady state 
value in magnitude. The inferential control system resul- 
ted in slightly better product concentration control than 
the conventional feedback system due to the fact that the 
simplified model response actually led, or anticipated, 
the process response. 

Several advantages of a digital computer over 
conventional analog control instruments were illustrated.. 
They include the possibility of communication between the 
direct digital control program and user written programs, 
the ease of implementing various control schemes by inter- 
connection of control loops, the wide range of controller 
constants available, and the ease with which control con- 


figurations may be modified. 
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1 INTRODUCTION 


eae Digital Process Control 


In recent years, digital computers have become 
commonplace in many industries. They have assumed a 
variety of duties, such as stock taking, data logging, 
and process control. Although they are capable of con- 
trolling processes directly, without the use of analog 
controllers, to date only about one quarter of all com- 
puters installed in the chemical industry really perform 
direct digital control. All others are used for applica- 
tions such as data logging or supervisory setpoint con- 
trol (4). This situation is likely to change in the future, 
for the advantages of direct digital control are numerous. 
Computer control offers more effective use of process in- 
formation and an almost unrestricted choice and combina- 
tion of control algorithms. Measurements and setpoints 
are available in engineering units and may therefore be 
determined and manipulated with high accuracy. The nece- 
ssity for human presence is reduced (8). 

The rapid increase of direct digital control 
applications has created a need for standard programs that 
can be readily adapted to handle a wide variety of control 
applications. It is now sarieye4d (30) that it should 
ultimately be possible to satisfy ninety percent of all 


user program requirements with standard off-the-shelf 
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packages. These standard programs will include direct 
digital control algorithms, feedforward control algorithms, 
heat and material balances, adaptive control, dynamic 
optimization of process behaviour, operator communications, 
log, scan, and alarm routines. The potential of the soft- 
ware development is illustrated by the fact that in 1966 
only fifty percent of software required by users was avail- 
able in form of ready made packages. 

This need for research and development led the 
University of Alberta to install an IBM 1800 Data Acquisi- 
tion and Control System in 1967. The installation is avail- 
able for research in the field of digital process control 
and as a service facilitvororsother research projects, It 
has been interfaced with an existing pilot plant-Scale evapo- 
rator in order to demonstrate some of the computer control 


applications. 


dpe Pilot Plant Evaporator 


The double effect concentrating evaporator had 
been built in this department three years ago. Prior to the 
installation of the computer, the evaporator was operated 
with analog controls. It was used for process control 
studies by Andre (2) and Wilson (31). 

Andre designed and supervised the construction. 
of the evaporator. He also tested various possible schemes 
for controlling the process. He used a forward feed mode 


of operation and concluded that the product concentration 
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could best be controlled by manipulating the steam flow. 
From practical cadunecaae considerations it followed that 
the feed rate would be determined by requirements else- 
where in the plant and was therefore not available for 
manipulation. The objective of Andre's work was the 
development of a mathematical model for the evaporation 
process. 

Wilson redesigned several sections of the unit 
to improve the operation. They included the recirculation 
in the second effect, the condensate collection, the vacuum 
control, and the steam supply. He implemented an analog 
control unit that adjusted the steam flow via feedforward’ 
control in response to feed concentration and feed flow 
disturbances, He concluded that feedforward compensation 
represented an important asset of the control system and 
recommended further work in this direction. After the 
termination of Wilson's experimental work, the evaporator 
was moved to a new building. On reassembling the apparatus, 
the major components remained unchanged, but several modi- 
fications designed to improve its performance and versati- 
lity were implemented. The feed system was rebuilt to per- 
mit blending of individual feed streams. Additional piping 
was installed that provides for more than twenty different 
modes of operation. Orifice plates and valve stems were - 
changed to allow for higher throughput rates. With the 


newly installed process-computer interface, the computer is 
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4 
able to read all essential process variables and exercise 
control over them, A switching system was also added to 
permit the choice between analog control and direct digital 


control of individual process variables. 


1.3 Objectives of the Present Study 


The present study was to provide a link between 
the previous work carried out in this department, and the 
field of digital computer control. The existing process 
equipment used in prior experimentation was retained, but 
the analog controllers were replaced by digital control 
algorithms, 

The specific objectives were to: 

a) develop a digital control system for the 

evaporator 

b) illustrate digital control configurations 

such as ratio and cascade systems 

c) implement digital feedforward control 

d) aenonberate inferential control 

The digital control system for the evaporator 
was considered in three parts. The first part was the 
new feed blending system which could be treated as a 
separate entity. The second part consisted of the evapo- 
ration process as such, comprising vacuum, level, and 
interstage flow controls. The product concentration con- 


trol formed the third part. 
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On implementing the digital control system for 
the evaporator, it was proposed to make effective use of 
control configurations that involved a network of indi- 
vidual control loops, such as cascade controllers and flow 
ratio controllers, Their successful application was ex- 
pected to demonstrate an additional advantage that direct 
digital control systems offer over analog controls. 

Wilson's experimental work had shown that definite 
improvements in product quality control can be achieved 
through the use of feedforward compensation for load dis- 
turbances. At the same time, the increasing use of digi- 
tal control systems demands general purpose software for 
various kinds of duties, Both these factors taken together 
made it desirable to develop a digital feedforward controller 
that would be assembled solely from available smaller pieces 
of standard software and would therefore be applicable to 
most processes of a similar nature. 

In certain industrial applications, it is desir- 
able to exercise control over a process variable without 
having to measure it physically (24), (14). Mathematical 
simulation of the process presents an alternate method for 
obtaining the value of such a variable. It was proposed 
to apply this type of control to the evaporation process. 

It would consist of a mathematical model of the process that 
would be given values of process inputs, and calculate 


the value of the product concentration at fixed intervals. 












ny 
P = ieee tlie wie teen P » es - fn 6 ao +B {- 
TO? @Szeye ieTsNGs 16822510 Of | 


a oe ani vor eee ee & 
AG SEY OVIZDSTISC Saga OF 9940q0 rg: eaw eo par 


oe thot ; O AITOWISA & boy lov rj Vasiz P eaolietughings L m3! 
+ 


“x9 eaw norisotiqyus Ivteassove sleAT earaeiecinens 
' b 


tooTlh tard sastnsvbs Isnoliibbs os ssateendiey oh 


baal so? sotjeaneaqnoa brsewaod best: 20 ‘eau aah 
gib to: eau gatesexyoni ofd . ombd omae PAR IA. soca 
102 avawsio0e ssogtaeq Jsetateg ebasaeb ewetaye foxan0a. J a3 ? 


aT : ae ¥ al ——e aie q 
rsriteg03 seed exotosi sseis diod ~éebeeb to ebata ewe OLTBY 


















ralioxgnos biewicibest Isaigih a qolavalesie idextasb 34. Dear 


se0siqg i9 m sfcalisyve mori -yilolee babiapamn: od bivow 2 and 
og sid gqs } stotsued? Dbivew baa ornwoioe traaboxse % to. 


beg Ss " 


_orvia 2eimte s te “asaesooxy 3 hal 
-"ytzob el st .efholbteotiags -Isixaaubrk atetzeo ol oie 
twotitw sidaizev aascc iq. 8 I9V0 (oxsa09, aauainniatie 


»@I)., ae Yi isstayag: <2 gona ae . 


to? bodsea: sismistis ana edness Palate anneeay: ‘ad 30. aaa 


r » + oe A i4 » WM 
$8ol3809e018aM 


ss 

s8ogorqg Baw II sidelrevy s dowe 2 sii ov ae anh : — 
- , 

| 7 
#es007q nolisyodevs «efio of : cy ak iqqs.o3 
r s3q NOLIBTOG 9 02 lor fIn0D . ° ons ke ele gs Ge 
ey 7 


a 


_ Z ar gy 
tess eesc osq on3 to Iabom pankioqeeny S Yo. abe i itn 










pee bos ,sdugat Bsssotg 30 eaulew ae 


ion | 







_ ; x. 
7 





A teas - 


{ rrostik o BS ‘ 72 —T a g add 
" olin , on? 28: £10 awnem 09. stot tail re 
Ie, wv 


7 


» 


This calculated, or inferred, value would be supplied to 


the digital feedback controller in place of the measured 


value. 


1.4 Experimental Work 


The experimental program for this study was 
designed to demonstrate the response of the evaporation 
System to changes in feed flow and concentration when 
operated under conventional feedback, feedforward, com- 
bined feedback-feedforward, or inferred feedback control. 

A series of preliminary experiments was also 
necessary to determine satisfactory steady state through- 
put rates, to calibrate all instruments, to verify the 
effectiveness of the proposed control schemes, and to 
demonstrate satisfactory performance of the computer con- 
trol system. 7 

During the experimental program, the analog 
control units were only used for recording of variables. 
To control the operation, the digital control facilities 
were used throughout. Start-up and shut-down werefacili- 
tated by a set of output stations that permitted manual 


setting of control valves. 
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2 LITERATURE SURVEY 


Zab Evaporation Theory 


Both Andre (2) and Wilson (31) have treated the 
theory of evaporation in detail. Andre has placed special 
emphasis on the aspect of heat transfer in evaporators and 
the possibilities of mathematical modelling. He developed 
a dynamic model that adequately represented the evaporation 
process, Wilson's report is oriented towards the operational 
aspects of evaporation equipment. He used a steady state 
material balance to provide the mathematical basis for his 
feedforward control unit. 

In a recent article by Vacop (28), equations have 
been @rived that describe the dynamic behaviour of boiling 
liquid in the tubes of a forced circulation evaporator. 

The special dimension of this treatment is that it considers 
pressure variations along the tubes. 

Various types of level control systems for vessels 
containing liquid have been described by Buckley (7). Among 
them is the averaging control system that has been applied to 


both effects of the evaporator in this study. 


2.2 Data Acquisition and Control Computers 


The general field of computers and process control 
is reviewed annually in Industrial & Engineering Chemistry 


(29), General trends of computer applications are discussed, 
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8 
and articles published on this subject over a one-year period 
are listed. 

A very practically oriented discusssion of the prob- 
lems that industrial computer applications present has been 
given by Control Engineering (34). The series provides in- 
formation on all aspects of computer usage for process con- 
trol, such as software, hardware, software-hardware interac- 
tion, various kinds of languages used in the development of 
software, as well as routine type programming techniques. 
Despite the continuing trend towards standardization of pro- 
cess control software, emphasis is repeatedly placed on the 
necessity for the engineer to be familiar with the operation 
of the computer in order to be able to make full use of its 
capabilities. 

A report on a compressor station controlled by a 
digital computer has been published by Friedli (11). It 
contains a qualitative description of the functions assumed 
by the computer. 

In a process control environment, it is generally 
necessary to assign different levels of priority to different 
types of programs. This ensures that process monitoring is 
not delayed or interrupted by non-process computer users. The 
concept of interrupts based on priority assignments has been 


treated by Gelder (13). 
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2.3 Computer - Process Interface 


Articles can be found that describe the techniques 
employed in converting digital signals into analog signals 
and vice versa. This signal conversion is an essential part 
of the computer - process interfacing equipment. Analog-to- 
digital signal conversion has been treated by Daley (9), and 
digital-to-analog conversion by Harris (15). An article by 


Abernethy contains information on multiplexers (1). 
2.4 Control Software 


In the development of control software, vendor and 
user usually complement each other. The computer vendor 
supplies the basic programs required for a particular type 
of application. The user supplements them with special pro- 
grams that suit his own specific purpose. In the present case, 
the main program representing the implementation of the direct 
digital control system has been supplied and described by the 
vendor (22). Additional descriptive material as well as 
Supplementary programs are being added as new applications 
arise (35). 

, A different application of the same computing system 
as that used in the present study has been described by 
Butler and Merryweather (8). The process in question is a 
glass fibre plant that is being controlled by single variable 
digital feedback controllers. It is stated that more elab- 


orate modes of control, such as cascading of several loops, 
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will be added to the system in the future. The application 
of the direct digital control system to a pilot plant scale 
evaporator, and the addition of feedforward and inferential 
controllers to the feedback controls, mark the present work 
as an additional step following the study by Butler and 
Merryweather. 

Bakke (5) has developed an algorithm that can be 
used in addition to multiple mode digital control to provide 
improved regulation for processes with inherent dead time 
periods. Implementation involves communication between the 
digital controller and a user written program at every samp- 
ling instant. This need for regular communication also 
arose in the present study on implementing the inferential 


controller, 
2.5 Feedforward Control 


A classical definition of feedforward control is 
contained in an article by Luhrs(21). It is a control action 
that converts data on process disturbances into corrective 
action before actual changes in the controlled variable 
occur. Thus feedforward is a control acting on news rather 
than history. The principle of feedforward control has been 
used in the power industry for the last forty years. 

Since the advent of computers, it has gained wider acceptance 
in the industry generally. 

The derivation of a feedforward transfer function 


for a process with one disturbance variable and one controll- 
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ed output variable has been presented by Luecke and McGuire 
(20). The first part of this article contains the direct 
analog to the method used in the present study to determine 
the relationship between the steam flow rate and feed con- 
centration disturbances, The authors continue from there to 
provide compensation for model errors and to include provi- 
Sions for limiting process conditions. 

A general method of designing feedforward con- 
trollers for processes described by dynamic time-invariant 
linear models has been presented by Bollinger and Lamb (6). 
The matrix approach is used to represent the process and the 
controller mathematically. The study is theoretical, no 
experimental application having been attempted. 

Feedforward control achieved with the use of 
analog components has been described by Wilson (31) and Michon 
(23). Wilson's experience with the evaporation process has 
been a point of departure for the present study. Michon has 
reported on the successful implementation of feedforward con- 
trol in electrolytic chlorine manufacture. In this case, 
the pressure of chlorine vapour downstream of the electroly- 
tic cells was regulated by using information on the condition 
of the cells, such as current density and liquid flow rates. 

A recent publication of the Foxboro Company des- 
cribes various applications of feedforward control as imple- 
mented with analog instrumentation (12). A feedforward con- 


trol system for a double effect evaporator is presented in 
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12 
which the feed tank level signal generates a series of control 
actions designed to adjust both steam and feed liquor flow 
rates, 

Reports by Woodley (33) and Tierney et al (26) 
describe applications of digital feedforward control. The 
object of Woodley's work was to regulate the temperature of 
the liquid leaving a heat exchanger. The feedforward con- 
troller was activated by disturbances in the liquid flow 
entering the exchanger. In his implementation, Woodley call- 
ed for feedback and feedforward calculations at specified 
intervals and combined the results. The concept is the 
Same as that used in the combined feedforward - feedback 
controller for the evaporator in the present study. 

In the work of Tierney et al, the product con- 
centration of a stirred tank reactor was controlled. Dis- 
turbances in feed temperature and feed concentration were 
considered, Controller equations, even for the feedback con- 
troller, had to be programmed especially for this process. 
For the present evaporation study, the use of the standard 
direct digital control package has eliminated this need, as 
the package provides algorithms for different modes of con- 
trol, and only the values of the controller constants have 


to be supplied. 
2.6 Inferential Control 


From the literature reviewed, it becomes apparent 
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13 
that the idea of inferential control has found acceptance in 
a variety of industries, 

Schaefer has described this type of control as 
applied to a lime stone quarry (24). The raw material is 
excavated at different locations and dumped onto a conveyor 
System via various crushers. The quality of the material at 
all excavating stations is observed frequently and entered 
as data into the computer. The combinedquality of all 
material reaching the point of convergence is then calcula- 
ted. A feedback signal system prevents input from points 
with off-specification material until such time that mater- 
ial from other input points balances the off-set, The study 
of the evaporation process followed essentially the same con- 
cept, with the additional feature that eventually, trouble- 
Some measuring devices can be dispensed with. 

Inferential control has also been treated by 

Harbaugh (14). Although the application is different from 
that of the present study, an interesting parallel can be 
drawn, as the concept and the implementation are the same, 
The system described is a steel rolling mill, i.e. a pure- 
ly mechanical operation. Input variables are measured and 
their values transmitted to the computer. The computer will 
use the mathematical process model to calculate the corres- 
ponding product quality. This calculated product quality 
is compared to the setpoint and the error used for control 


action, 
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A textbook approach to inferential control has 
been undertaken by Kelley (19). He presents the mechanism 
producing this type of control as a Ziebolz controller, 
which he defines as an automatic regulating device whose 
control action is based on information obtained from a 
Simulation of the process under study. The main difference 
between the closed loop feedback controller and the Ziebolz 
or inferential controller is that the former is activated 
by instantaneous error functions, whereas the latter der=- 
ives its control action from predictions obtained from a 
mathematical model. Kelley insists on the fact that the 
Ziebolz controller is more effective than a conventional 
feedback controller because it works with the future of the 
variable under control, instead of with its present or past. 
The problem of controlling a process by this method is 
essentially reduced to a problem of simulation. Kelley's 
treatment remains qualitative, as no specific experiments 
are presented. It can be regarded as a starting point for 
future work on this subject, which is, as the author states, 
still in the infancy of its development. 

Several industrial organizations are using or 
developing control schemes similar to those used in this 
work, However, for competitve reasons most of this work 
has not been published or is reported in only very quali- 


tative terms. 
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3 PROCESS EQUIPMENT AND CONTROLS 
3.1 Introduction 


The double effect evaporator that formed the ex- 
perimental apparatus is represented schematically on Figures 
3.1 and 3.2. The detailed equipment configuration. is shown 
on Figure B-19 of Appendix B. This final configuration has 
evolved gradually over the past two years, in a continuing 
effort to extend the capability and flexibility of the equip- 
ment. The extent of the modifications can be appreciated 
by comparing Figure B-19 to the flow sheet of the same evapo 
rator as used by Wilson (31). Apart from the changes notice 
able on the figure, the equipment has been moved to a new 
building and interfaced with the computer. With the comple- 
ted interface, all essential process variables can be reg- 
ulated on an optional basis either under local analog, 
analog supervisory, or under direct digital control. During 
the present study, the analog controllers were only used in 
some preliminary runs, The equipment as assembled in the 
new building is described in the operating manual available 
in this department (32). 

It was also decided to abandon the use of sugar 
solutions because of their objectionable sticky properties 
and their tendency to ferment. AY experiments were con- ~ 
ducted with triethylene glycol solutions. Triethylene gly- 
col is a colorless liquid that meets all the safety criteria 
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Legend to Figures 3.1 and 3.2 


<> process stream 

V service stream 

[eee] digital controller 

C2) signal transmitted to computer 

fame focal control 
process link between Figures 3.1 and 3.2 
instrumentation link between Figures 3.1 and 3.2 

CC concentration controller 

CT concentration transmitter 

DVPC differential vapor pressure controller 


DVPT differential vapor pressure transmitter 


FC flow controller 

FFC feedforward compensation 
FT flow transmitter 

LC level controller 

LT level transmitter 

PC pressure controller 

an pressure transmitter 

RC ratio controller 


TEG triethylene glycol 
TC temperature controller 


dt temperature transmitter 


Table 3.1 
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to be considered in a student laboratory environment. Its 
volatility is very much less than that of water, and it does 
not exhibit a boiling point rise in the temperature and con- 


centration ranges considered (27). 


3.2 Design Changes 


3.2.1 Feed System 


The evaporator feed system had been completely 
redesigned when the apparatus was reassembled in the new 
building. As had been pointed out by Andre (2), in indus= 
trial installations the feed rate to an evaporator is usually 
determined by conditions prevailing elsewhere in the plant, 
and is therefore not available for manipulation. The feed 
System as used in the present study is shown as a separate 
unit in Figure 3.1. The network of controllers that has 
been successfully employed illustrates the flexibility of 
the direct digital control system. 

The equipment configuration and controls shown 
allow for independent changes in concentration, flow rate, 
and temperature to be introduced. Seven controllers were 
required cs regulate these three variables. Detailed des- 
criptions of the individual digital controllers are con- 
tained in Appendix D. 

The feed consisted of two streams: ‘triethylene 


glycol solution (F 12) and water (F 11) that could be 
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blended in such proportions as to yield the desired concen- 
tration (C 1) of the total feed flow (F 8). The solution 
flow controller received its setpoint from the total flow 
controller. The ratio controller used the measurement of 
the solution flow and the ratio constant supplied by the 
concentration controller to determine the setpoint for the 
water flow. Setpoints for the total feed flow rate and 
concentration were entered by the operator. Provision was 
also made for operating the concentration controller manu- 
ally in cases where the measurement obtained from the refrac 
tometer was found faulty. Heat exchangers were available 
to control the temperature of the two individual streams. 
Satisfactory operation was achieved with a constant feed 
temperature of 195°F, which ensured the prevention of boil- 
ing in the feed piping even under severe load disturbances. 
The flow controllers were first tuned ‘ndividu- 
ally during the testing period. The control constants and 
filters were later adjusted as need arose under actual 
operating conditions, such that disturbances could be 
introduced without resultant cycling. As the feed temp- 
erature setpoint was kept constant at all times, the two 
temperature control loops had to respond only to load changes 
resulting from steps in flow rate. Over the extended experi 
mental period, the momentary temperature disturbance was 
kept to + 4°F by the two temperature controllers. A list 


of controller constants is contained in Appendix D. 
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3.2.2 Evaporation System 


On reassembling the apparatus, provisions were 
made for operating the evaporator in at least twenty diff- 
erent modes. Although this particular feature did not 
affect the present study, it made the unit adaptable for 
various kinds of investigations. 

The closed loop operation that became feasible 
with the new feed system involved recirculation of the 
product to the feed tanks. In order to avoid elevated 
liquid temperatures in the tanks, a cooler was designed and 
installed that cooled the product to room temperature. 

To reduce the noise in the recordings of the 
bottoms flows from both effects, the orifices were repo- 
Sitioned so that they are now located upstream of the flow 
control valves. 

It was also felt that the evaporator in its 
present condition could handle much higher throughput 
rates than those used fy Andre(2)andWilson (31). Consequent 
ly, the orifice plates and valve stems in nearly all lines 
were changed to allow for higher liquid and steam flow 
rates, The equipment was then operated at an increased 
capacity, which created conditions more typical of indus- 
trial practice. It was also expected that the increased 
rates would reduce the total hold-up time of the system. 
Appendix B contains a list of orifice diameters and valve 


stem sizes as well as calibration curves for all recorders. 
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Other new installations included a drum in the 
steam supply line to dampen out pressure fluctuations in 
the header, cartridge filters in the recirculation lines 
from the feed pumps to separate accumulated solid particles, 
and additional piping in the product line to permit pump- 


ing of product to three different tanks. 


3.3 Evaporator Controls 


3.3.1 Averaging Level Control 


The control system for the evaporation process is 
shown on Figure 3.2. In the case of liquid levels, at 
least two choices of control criteria are available. First- 
ly, a level can be controlled at a constant value. The 
outlet flow from the vessel will then have to vary in the 
Same manner as the inlet flow. Secondly, the level can be 
allowed to rise or fall within the physical limits of the 
process equipment to absorb disturbances in the inlet flow 
and only slowly transmit them to the outlet flow. The liquid 
levels in both effects of the evaporator were controlled by 
manipulating the outlet flows. 

Simulation studies carried out in this department 
had shown that the product concentration is affected to a 
considerable extent by sudden changes in the interstage flow 
rates (10). It was therefore decided to use averaging Tete 
control of the type described by Buckley (7). It consists 


of a level control loop with proportional band equal to the 
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allowable total fluctuation in liquid level and no or very 
weak integral action. With this type of control, the outlet 
valve will be completely shut when the level is at the 
lower limit, and completely open when it is at the upper 
limit. Buckley's mathematical treatment shows that maximum 
regulation of outflow against inflow or pressure disturb- 
ances can be achieved with a cascade level+flow control 
System. Since in the present case, outflow regulation 
was critical to product concentration control, both level 
controllers were cascaded to flow control loops, which 
resulted in satisfactory smoothing of outlet flows. Actual 
response curves of both levels and bottoms flows resulting 
from feed flow disturbances are presented in Appendix D. 

In the earlier stages of the experimentation, the 
first effect liquid level control presented difficulties 
the real cause of which was only found much later in the 
program, The sightglass used for level readings was connec- 
ted to the vessel by quarter inch tubing. Apparently, 
vapour would condense and collect in this tubing, creating 
a pressure difference between the vessel and the sightglass. 
As this pressure difference reached a certain value, the 
water column was blown out of the tubing, and the resulting 
sudden pressure increase in the sightglass caused the indi- 
cated liquid level to fall instantaneously by approximately 
1% inches. This sequence of events was repeated every 


thirty to thirty-five minutes. Once the cause was apparent, 
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the sightglass was removed, and the level transmitter connec- 
ted directly to the top and the bottom of the vessel. Al- 
though the resulting signal was slightly noisier, the improve- 


ment in level control was obvious, as may be seen on Figure 


Seas 
3.3.2 Basic Controls 


The vacutmin the second effect and the cooling 
water flow to the condenser were both regulated by simple 
feedback loops. A differential vapour pressure controller 
was available to remove accumulated air from the second effect 
Steam chest. The condensate from the first effect steam 
chest and from the rundown tank was removed by local level 
controls, 

A cascade controller admitted the condensate from 
the second effect steam chest to the rundown tank. In the 
original installation, the orifice was located downstream of 
the control valve. The pressure in the chest ran at 4 to 5 
psig, while the rundown tank was kept under 15 in Hg vacuum. 
This relatively high pressure drop caused severe flashing 
of the almost saturated liquid that made flow recordings 
noisy. The valve and the orifice were interchanged, such 
as to place the orifice upstream of the valve. This arrange 
ment resulted in considerably improved flow recordings, as- 
is illustrated on Figure 3.4. 

Another problem arose from fluctuations in the 


vacuimheader. The vacUumservice in the building was main- 























a ' , 
71 s) 
= i ‘5 
ws —. % 
Po we wu 
at. 
7 


tire ? -% 7 a <x be 
‘ 3 aaa ee ‘ —— = 
-9enros tos3imenets ftsvel ant. Sen bevomes mal Ben. guage 


“s : - / 7 x ah at 4 bs 
ASA .fe2eeyv oft Xo moSsod sds bas gor ons os a 29240 ® 
= ’ “ ‘ 


5 
-svorgmi sis ,s9latoa.yi3dgile vaw fengte: gnisiemox ‘eda 
este] — ¢ vam en , amokvdo saw towss09, lever i 


78 Jostis . 5958 Fe 

[ues x Asod & ones zaeasbnos — at iho 

FINO Swse9TG Tweqsy- {sts ugT9 3325 A. | -8qg0F soa ~ 

$06226 | * 2 99 sovtt vie. hazel UMD Os evomox ‘bs satire: 
sits, JevEY ed3 mov? 52 senebaoo ont » Seeder nse3a 


0 bevomet abw Aned pwobrur ers maoz? ee tee oni: > 


pe >» “salon 


’ . 2? 


r 


Ott stsenebno> siz beIsimbs sel homo, ncnmaan’ 
dia a ainsi mwobtuyt eis oF taon fo mete apsdie a a 
ob bedacol saw 4oitixo ods fobs iisvenk Ientgt: - 

2 of # ga tzario orls nk omudeerg eit vovinw J foxsaeo 9 
uuosy gH ot 25 zebra 298 I esw wines mwobotry' ods “otk: Ebel 
enines | teve2 basse qoxb STUReSTgG ‘ight’ lovisatew ah 
egnibiooes voll sbam “teria blupti bosexusoe -" oa: = 


ar 


=< 


Ta 
‘gous , -egneno1ss ni eaxow aaiitxo od bas oviev set 
egfe tr eldT .avisy 983 20 smoxzeqy eatRivo add 69s sq 0: 
88 ,egaibrooss wolk Sevorqat elders tana ‘e = jiueet % 
7 > 
; < can .. oa 





2) 





ee i 


soiiaeaenote aeestepeni spe ortntatnasteicimriis a 
406 LORPALG ANG Wy YONG = eee 








— 








a ee ee CM a ete 
aE ST Ban hak = ones ES 






















































































d Levels 


iqui 


Recordings of L 


ht glass on L14 


wwiLthsae 


above 


ight glass on L14 


without s 


below 


Figure 3,3 


| es 
. 
é 
‘ 
fi 


. a 1 : i ae 

> : : a : a a 

rf aie fe +f Re . 
+f "iF ‘at AI nt ir % ‘i 


Bye 















.y 


es 
“it yb ail a 


eo 


- om me 
. “a 
—_ -« 


i in nee re ved 
#3 ein y itaeitmaiele) 
‘ a mi 7 x 
| ene Ld | | | 
hy 1 iq A ee 2 Ths 
- } it” at phar tyi “3 4 ' 
‘4 T ae il a " 
y vr " y siore 
vt r ay 


Bais 
‘ , Molt We aN 
‘ bi a a 


—< 4 { ~ . i a 
WT he Ge | ih the i ne 
Wan este ea 2 Ee 








— ‘ 
Sm] Til Fare 
F eutne 4p 


ie ap 


~ 
: =~ 
| ‘ies 3 f vd ona - 
“9 246479 Iigte iy 
i 


ov 

? x y 
» 2enig siyhe 3am athe sid fad 
Se ene 


, 
. wen - Bee 











a : 





26 





Se a a i ree ; aap cea 
SS ice 2 Ss Saree ees 3 GRAPHIC GONTAOLS COAMPORATION —5 fay : FRE SS SAREE SS) | EET RIL x : 
aS CONTACTS CORFORATION 1OQ3 8 = 100 -icuarr s00rk 


SSS = 




















DN een ee ee 
SSeS oS So aimee weceomeroes Sevmeionee 











































































































ee 





Hh i 
: ; a pa) 
Seen teen dan 





Figure 3.4 





Recordings of Condensate Flow F5 
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Recordings of Vacuum P22 
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Figure 3.5 
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tained by a pump that operated under on-off control, supply- 
ing a vacuum of 20 to 24 in Hg. The fluctuations caused by 
this pump affected the operation of the plant to an extent 
that had not been anticipated, as they propagated through 
the entire process. The problem was solved by constructing 
a pulsation damper and inserting it between the vacuum header 
and the condenser. The damper consisted of the 6.1 ft3 
steam condensate tank with 15 feet of quarter inch tubing 
on each side. The tank was made available by leading the 
condensate from the first effect steam chest directly to 

the sewer. Figure 3.5 illustrates the results achieved. 

The values of several additional variables were 
available to the computer to be used in material and energy 
balance calculations. They included the condensate flow 
from the condenser and the temperatures of all essential 


streams, 
3.3.3 Product Concentration Control 


The product concentration was controlled by mani- 
pulating the steam flow. The basis of this control system 
was a cascade concentration-flow controller as shown on 
Figure 3.2. Additional features of the concentration con- 
trol system, such as feedforward compensation and inferen- 
tial control, will be described in Chapter 5. | 

The constants for the feedback controller were 


determined empirically during a test run. Values reported 
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29 
by Wilson (31) were used as a starting point. The prelimi- 
Nary experiment showed that these constants gave rise to 
an extended period of oscillations in response to a step 
type load disturbance. It was possible to reduce the osci- 
llations considerably by decreasing the gain. The reset 


time was then increased to speed up the elimination of off- 


sets, 


3.4 General Operating Description 


At the outset of the experimental program, a 
Steady state of reference was established that was used as 
‘a pivot for all experiments. A set of data for this reference 
State is presented in Appendix A. For all experiments con- 
ducted, a forward feed mode of evaporator operation was 
employed. Reference is made to Figures 3.1 and 3.2 for 
the following description. 

The feedflow consisted of two merging streams, 
water and triethylene glycol solution. Mixing took place 
in the pipe between the point of convergence and the first 
effect. The proper feed concentration was established by 
manipulating the ratio. of the two flow rates. The two feed 
streams were heated individually by passing them through 
steam jacketed heat exchangers. The vapour produced in the 
first effect was condensed in the steam chest of the second 
effect and its condensate admitted to the condensate collec- 


tion system, The bottoms flow from the first effect reached 
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the second effect at the suction side of the circulation 
pump. The liquid in the second effect was circulated 
through the vertical tubes at a velocity of 4.2 ft/sec. 
Product was withdrawn at the discharge side of the circu- 
lation pump and pumped to storage tanks via a cooler. 
The vapour produced in the second effect was condensed in 
the condenser and admitted to the rundown tank whose level 
was kept constant by pumping the condensate back to the 
water storage tank. The first effect consisted of a cal- 
andria type heat exchanger with a total surface area of 
9.88 £t2, The second effect was a vertical tube exchanger 
with a surface of 4.1 ft2. 

More detailed descriptions of the equipment can 


be found in references (32), (31), and (2). 
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4 COMPUTER SYSTEM 


4.1 Computer Equipment 


4.1.1 Summary of System Components 


The data acquisition and control system that was 
available for this work is represented schematically on 
Figure 4.1. Apart from the central processing unit, the 
components most vital to this study were the analog inputs, 
the analog outputs, and the process operator's console. The | 
card read-punch unit, the printer, and the plotter were used 


for purposes of data processing. 
4.1.2 Process Operator's Console 


The process operator's console consisted of a panel, 
a typewriter, and a keyboard. The panel contained a set of 
push buttons that facilitated operator communication with 
the computer. Actual data input and output proceeded via 
the typewriter-keyboard unit. 

The functions of the process operator's console 
can be summed up by stating that it presented the physical 
means of communication between operator, computer, and pro- 
cess. More specifically, it served to select and address 
process variable records, to display status information on 
process variables, and to supply or change parameters used | 
for controlling or recording process variables. Process 


variable records will be discussed in Section 4.3.3. The 
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33 
execution of disk resident user programs could be initiated 
from the process operator's console. 

All analog signals transmitted from the computer 
to the process had to pass output stations located near the 
process operator's console. These stations could be used to 
break the signal transmission circuit and to select valve 
settings manually. They also displayed the control signal, 


in percent, going to all computer operated control valves. 
4.1.3 Process Input-Output Handling 


The functions of the analog input equipment is the 
collection of analog data and its conversion for presen- 
tation to the digital computer. The analog data in this 
application are current signals emitted by standard indust- 
rial transducers. The hardware to convert this current 
signal into a digital value consists of a signal conditioner, 
an analog filter, a multiplexer, and an analog-to-digital 
converter. The signal conditioner converts the current in- 
to a voltage signal. The analog filter reduces noise in 
the voltage signal arriving at the converter. The flying 
capacitor multiplexer contained in the analog filter cir- 
cuit allows voltage signals to be transmitted to the analog- 
to-digital converter. Since this converter operates at high 
speeds, it can be time-shared among a number of input voltages 
through an analog relay multiplexer, With this system, the 
computer can read up to one hundred analog inputs per second. 


The analog-to-digital converter used allows for a fourteen 
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Sequence of Signal Conversions 


Device 
variable measured, e.g. flow 
measuring device, e.g. orifice 


transducer, e.g. dp-cell 


signal conditioner, e.g. resistor 


analog filter 
time constant 0.33 sec 


signal converter 
analog-to-digital (14 bits) 


program manipulation 
e.g. range expansion 


digital calculations 
signal converter 


digital(10 bits) -to-analog 


current output station, facility 
for reading and manipulation 


current-to-air converter 
final control element 


e.g. valve 


Table 4,1 


Output Signal 


Lbs/min 
in; H2Z0 
10 - 50 ma 


1 - 5 volts 
La= 5 cace 
6553-32767 
0 - 32767 


0 - 32767 


10 volts 


Orc 
i] 


4 - 20 ma 


3 - 15 psig 


Lbs/min 
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bit resolution. 

The function of the analog output equipment is the 
conversion of digital values into current signals acceptable 
to analog instruments such as current-to-air converters. 
Digital data is stored in output channel registers connected 
to digital-to-analog converters. The digital values are con- 
verted into voltages with a resolution of ten bits. The 
voltages so obtained are amplified and then converted into 
4-20 ma current signals by the current output stations. 

A typical sequence of signal conversions is shown 
on Table 4.1. Detailed descriptions of the equipment involv- 


ed are contained in publications of IBM Corporation (17, 18). 


4,2 Operating System 


4.2.1 Time Sharing Executive System 


There are two distinctive features of real-time 
computer applications as compared to all other applications. 
Firstly, a continuous need exists for communication between 
the computer and the process under control. Secondly, the 
computing system has to keep pace with the process, i.e. it 
has to update records, measurements, and output signals at 
regular intervals. The Time Sharing Executive System is a 
group of programs that permit the control of processes while 
at the same time providing an off-line monitor system for 
data processing and scientific computations (17, 18). 


The framework of the Time Sharing Executive System 
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36 
is provided by the executive program which resides perman- 


antly in core storage. It consists of the most important 
programs, sub-routines, and storage areas that are essential 
for real-time applications, Most other portions of the 
System are stored on disk and are only brought into core 
when required for specific functions. 

Since the process under control will not utilize 
all of the computer's time and resources, a certain amount 
of background duties can be performed by the computer while 
monitoring the process. The notion of time sharing implies 
the sharing of computer time, storage, and input-output faci- 
lities among on-line and off-line programs or users. 

The core storage size is not sufficient to contain 
all instructions required for the execution of all functions 
at any one time. It is therefore necessary to break down the 
total set of instructions into smaller pieces that reside on 
disk and are available for immediate transfer to core storage 
when required. These smaller pieces of a set of instructions 
are stored on disk in executable core image format and are 
referred to as core loads. 

The Time Sharing Executive System consists essent- 
ially of two groups of programs; namely, control programs, 
and processing programs. In general, control programs super- 
vise the execution of processing programs and user programs. 
The three control programs within this system are the 
Temporary Assembled Skeleton, the System Director, and the 


Nonprocess Supervisor. The System Director is the nucleus 
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37 
of the skeleton executive, It exercises control over all 
operations connected with real-time operation. It handles 
the transfer of control between on-line and off-line opera- 
tions. It governs the transfer from one user specified core 
load to the next. This function was important to the imple- 
mentation of the inferential control system described in 
Chapter 5, where a number of core loads had to be executed in 
a given sequence. The System Director recognizes, controls, 
and directs the servicing of interrupts. Interrupts are an 
essential feature of an on-line real-time control system 
operated on a time sharing basis. Different priority levels 
are assigned to all types of programs according to the rela- 
tive importance that is attached to the time of their execu- 
tion. The priority scheme thus established basically causes 
a given series of programs to be executed selectively, rather 
than sequentially. With respect to the present study, it was 
important that duties related to the control of the process 
would be performed by the computer with preference over other 
duties to prevent disruption of the process. 

The System Director also supervises the operation 
of interval timers. These clocking devices allow for a 
given program to be executed at fixed intervals, a feature 
that has been used to advantage in the implementation of the 


inferential controller. 
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4.2.2 Programming Considerations 


As the skeleton executive occupies a considerable 
part of core storage, the area available for the execution of 
user programs, referred to as variable core, was restricted. 
This factor had to be considered when the special programs 
needed for this study were written. In general, the programs 
were divided into small logical units that could be stored as 
core loads and executed in the proper sequence. 

Since there was only a limited number of priority 
levels available on this system, many process related pro- 
grams shared the same level. These included the program 
monitoring the process operator's console. During the time 
of the experimental program, efficient coordination of various 
operations involving these programs became imperative. As 
an example, it was not possible to use the process operator's 
console while graphs were being produced on the plotter or 
material balance calculations were being performed. 

Facilities were available to accumulate process 
data in files established in the disk storage unit. With this 
feature, graphs of process variables as functions of time 
could be produced on the system plotter without manually 


recording the data. 
4.3 Direct Digital Control Program Package 
4.3.1 Direct Digital Control 


Direct Digital Control is an advanced form of 
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computer process control. Computer generated signals are 
used to control the process directly, without the inter- 
vention of analog controllers. 

The most important characteristics of direct digi- 
tal control can be summarized as follows. 

Controller adjustments can be carried out with 
ease from the process operator's console. New controllers 
can be added to an existing system by simply inputting the 
appropriate parameters from punched cards. Measurements, 
setpoints, and controller outputs can be determined with high 
precision. An almost unrestricted range of controller con- 
stants is available. Advanced control algorithms can be 
applied through existing communication links between the 


digital control software and user written programs. 
4.3.2 Direct Digital Control Program 


The direct digital control program is a system 
program that resides permanently in core storage. It monitors 
all data acquisition and control functions of the computer. 
The initial activation of the program proceeds via the main 
eetsare and the card reader. Its execution is initiated by 
interrupts at intervals specified by the user. In this 
particular installation, it was executed every second, 

Special communication programs are available that service the 
process operator's console and print alarm messages. Commu- 


nication with the direct digital control program is provided 
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for by utility programs. With the aid of these utility 
programs, the user can load the process variable table or 
parts thereof from input devices into core storage, and 
obtain the contents of the table through output devices at 
any time. The table will be described in Section 4.3.3. 

A typical set of functions handled by the direct 
digital control program are the following. 

Measurements of process variables are obtained 
at user specified intervals, compared to given alarm limit 
values, filtered according to user requests, and stored in 
the process variable table. Valve commands are calculated 
using control algorithms chosen by the user. Commands are 
transmitted to the process via analog output equipment. 
Calculated controller outputs are stored in the process 


variable table for use in other variable records or programs. 
4.3.3 Process Variable Table 


The direct digital control program is organized 
around a single table of data referred to as process vari- 
able table. This table contains a variable record, also 
referred to as loop record, for each variable that is con- 
trolled, manipulated, or recorded. A variable record con- 
tains coded instructions, parameters, and logical flags that 
define the action to be taken by the direct digital control 
program with respect to the variable in question, The pro- 


cess variable table can be brought into core storage via the 


0 sy as ‘ ns . 





















yoiliay sessed’ to bks erfs c3IW - ammaegong & 
to eldat sidatteay eaescorq 903 baod 89 senw" “ng 3 


soo. otnk 295 veh tugnt mont cadihial 


° 


— 


48 esoiveh iuqsuo dguords efdad: ond 20 einssaeo ah kt sade 


oa2 nit beditoeeb sd iliw aides. ed 
y! . g acl emolssawt to! Jee ‘faotqys A 


uniwollot sz oe ms 12079, forssoes fostakb 
Ag 


ba o{detrav sea00rg, to. etaogexuesoM ce ae 
P a i, 


.% 


aro2 7 & ,asvasupez seseu oF guthro598 & bored 12% neiaic 8 i a 
- i we a Pate a 
be _avs ebnacmoo svisV . oldest, sidattsy eae oma “ 


. . : a 
i “ 


7 .zeeu odd yd nesodo em d3isogi 5 fossa eo 
. ef = 


rr pups Jjyqsauo a iste eiv — ait o2 aaa ae 


ba? Do +eaqm02 el avistrn | bathiaags, roan 3B . s 
Ms "i 


4 
a ’ r = 5 cP a 
; et aldaizvay torlto alesse wr aids: 
.< .* mf 
5 oe 
a é 
Yee? ~ fila bent? ob oeunwll cS 4 
LUZ IGBITAY 2239023 rs 
‘ q he 
_ - PY 


ossg0rq foxtmop statb taextb siiT 
-ixey a2sc0rq sa 09 bervoter atab * v0 aida algase » Bnw03 28 


c ,bsocer sideltav a eniedaop ods oid?  ekda es ; 


yy. 
is sidsizevedoss ‘Sok aioe qool. vr Ox 
“00 Sbroost sidsaivtev A Lbobr006% 10%, boasts q a aa Stier 


—_ <n sere y, anit 6 


3a2 ¢gsnid feolgol bas .exaasmab3, 








me eae {fesigtb S297Lb ods: si asain’ od eile 


“af omy : wit. ~ .pokseoup me afdets S ee B ll +2 : W 19530" 
: | ial | et, ] 


“ 


wer = se 
2 sida: etidal 
7 gm - 


41 
card reader in a prescribed format. Most entries in the 


table may be altered at any time from the process operator's 
console, 

An example of the process variable table has been 
included in Appendix D. At execution time, the direct digi- 
tal control program scans the process variable table to deter- 
mine the action to be taken with respect to every single pro- 
cess variable at this particular instant. The parameters 
in the table will define the action required. In particular, 
they will specify the frequency at which each variable 
requires the program's attention. In addition, the table 
contains all information required by the program for servic- 
ing each process variable, such as the source of data, the 
range of measurements to be expected, alarm limits, control 
parameters, setpoints, filter constants, and the destination 
of the output. | 

Although the process variable table is loaded into 
core storage and maintained by utility programs, it can be 
altered at the discretion of the user from the process 
operator's console. A special set of sub-routines is pro- 
vided that allows user written programs to access the pro- 
cess variable table, to extract information from it, or to 
supply information to it, All parameters contained in the 
table may be displayed on the console in hexadecimal or deci- 
mal notation, or in engineering units. The various means 
of communication with the table are illustrated on Figure 


4.2. 
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Three basic kinds of variable records have been 
employed in the present study; namely, data acquisition 
records, control records, and data accumulation records. A 
data acquisition record references a variable that is to be 
measured only. At time intervals specified in the record, 
the program obtains a measurement and stores it in the record 
from where it is available to other records or user pro- 
grams. A variable that is to be controlled or manipulated 
is represented in the table by a control record. This record 
contains specifications for the required control action and 
for the destination of the calculated output. The output 
may go directly to a final control element, or it may be 
placed into the record of another variable. Data accumulation 
records may be used to collect a series of measurement values 
from any process variable, Certain options may be exercised 
as to the destination of the accumulated data when the record 
is full. The data may be disregarded or written into a file 
On disk, or the record may be taken out of service, thus 
Saving the data it contains. 

The calculation of control output values is based 
on error signals. If a variable is referenced by a control 
record in the process variable table, the program will per- 
form control calculations as specified in the record. The 
basic control modes that may be selected are proportional, 
integral, derivative, or a combination thereof. The algo- 


rithms the program uses are listed in Appendix C. In addi- 
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tion, a number of special routines are available that permit 
averaging of several variables, ratio control, and material 
integration. An output bias term is provided that permits 
results of advanced control calculations to be added to the 
controller output calculated by the standard routines. This 
feature has been used in the implementation of feedforward 


control, 
4.3.4 Distinctive Features 


Two characteristics of direct digital control that 
require specifications not encountered with analog controllers 
are the availability of digital filters and the discrete 


sampling techniques, In the present study, filters were 
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employed for two purposes, Firstly, they smoothed measure- 
ments to reduce noise or fluctuations. Filter constants were 
selected based on consideration of process time constants and 
then refined empirically. Secondly, filters were used in 
conjunction with control calculations to obtain a specific 
predetermined shape of a controller output curve. This 
latter application will be expanded on in Section 5.2. 

The two kinds of filters that the direct digital 
control program provides are exponential filters and Union 
filters, the latter named after their inventor. The expon- 
ential filter reduces a step change in the raw measurement 
to an exponential curve in the filtered measurement, thus 
acting as a damper. 

The Union filter uses a velocity Limiting procedure. 
It allows the rate of change of the measurement to increase 
as the error increases, At steady state, it prevents the 
measurement received by the computer from passing a speci- 
fied limiting distance from the setpoint in any one scan 
period, In this way, spikes in measurements are effectively 
reduced, This filter was used in the nested control configu- 
ration of the feed system, where it prevented noise from 
propagating. 

The algorithms for the filter actions are listed 
in Appendix C, 

It is expected that digital filters will become 


increasingly important in future work, as more practical 
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experience will more closely define their assets and their 
range of applications, 

As has been mentioned, the frequency at which 
variable records are polled by the direct digital control 
program is specified by the user. In the present study, these 
intervals were determined according to the expected rate of 
change of different variables, and according to the destina- 
tion of the measurements. They ranged from one second for 
flow controllers up to 64 seconds for data accumulation records. 
In systems of nested control loops, the inner loops were 
polled more frequently than the outer loops. 

The program was basically capable of reading up to 
100 analog inputs per second. Since the computer was simul- 
taneously serving other process users who required analog 
input readings, it became desirable to distribute the work 
load for the evaporator evenly over long periods of time. 
The polling sequence resulting from these considerations is 


shown in Appendix D. 
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5 DIGITAL CONTROL OF PRODUCT CONCENTRATION 
5.1 Mathematical Model 


A mathematical model was needed for both the feed- 
forward and the inferential controller that were to be in- 
corporated in the digital control system. The model was 
intentionally kept simple, in order not to obscure the con- 
cepts of implementation. Material balances and experimental 
steam economy formed the basis for this model, which was 
developed with reference to Figure 5.1. In this model, it 
was assumed that the feed entered the first effect at a con- 
stant temperature, and that the levels in both effects were 
constant at all times. The former assumption held true through- 
out the experimental program, as feed temperature controls 
were available. The latter assumption did not hold true 
during the experimental period, as with the successful im- 
plementation of averaging level control, both levels were 
allowed to vary. As a consequence, it was expected that the 
model would give a slightly faster response than the process. 
The dynamics of heat transfer in both effects was assumed to 


be negligible compared to the other process time constants. 


Overall material balance 


WSC 


WL = O1 + O2 + B2 (since WST 
01C) (5-1) 


and OL 


First effect solute balance 


WI*CI = B1XC1l + H «gol (5-2) 
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TIME IN MINUTES 





$9 
Second effect solute balance 
BI¥C2 = B2*CO + H2*dCO (5-3) 
Steam economy based on an overall balance 
O1 + 02 = WST*STEC (5-4) 


For equally distributed evaporation (experimentally 
verified) 


O1 = 0.5*WST*STEC = 02 (5-5) 
The 0.86 minute transportation lag between the effects has 
been neglected in this model, so that 

C2 = Cl (5-6) 


at all times. 


Once the model was developed, it was necessary to 
select a time interval at which the output of the model would 
be evaluated and used for control purposes. By simulating 
the process and plotting response curves obtained for vari- 
ous computation intervals, it was shown that the response 
curve does not change if intervals of less than one minute 
are used. The results of the simulation are shown on Figure 
5.2 for sampling intervals ranging from ten minutes to less 
than one minute. The curves represent the response to a 20% 
step change in feed concentration. Considering these results, 
it was decided to use a one-minute computation interval for 
all simulations. The effect of sampling the measured vari- 
ables at one-minute intervals is equivalent to introducing 
a pure time delay of approximately thirty seconds. This 
is the same order of magnitude as the neglected process 


time delay. 
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50 
Equations (5-1) to (5-6) can be combined to yield 
two differential equations expressing the concentrations in 


both effects as functions of time. 


22 =| wr *c1 B ar Z 0. sksTECHuST *C2 (5-7) 
dco 


nz | WI - 0. SkSTECKNST| *C2 


dt 
. byt |W ~ STECHUST| *CO (5-8) 


Equations (5-7) and (5-8) constitute the mathe- 
matical model of the process that was used to design both 


the feedforward and the inferential controller. 


5.2 Feedforward Controller 
5.2.1 General Description 


It was desired to provide feedforward compensa- 
tion for feed flow as well as concentration disturbances. 
The digital controller had to recognize these load distur- 
bances and correct the steam flow to such a degree as to 
prevent them from propagating through the process. The en- 
tire feedforward controller was made up of three variable 
records that could be manipulated easily from the process 
operator's console. No additional programming was required. 


The controller is shown schematically on Figure 5.6. 


a7 
be > 
- oe? 














i. 
- = *% a + a 
> of. neo (e-¢) oF Ci-2 y eftoljs aopa - fn ( 
- a i Ee 

it anoitsracsonn> art goteastqxs enoliaups Iatsnardiab ¢ yw 


* ~ - ~ >. =" a Ns 
32 k0 anokionpi ¢a 23 sento azo 


P __ 2 ry 
0G ae 7) teres « i | Li aoe z 5 
~ FP a ‘it bel Gaal ¢-3 .2 
_— 

| £ 09 

“ oper 1. aa 

WAQRTS*e 0. = Tay = 

: 2 
4 


- 
— a 
) Drawrotbest $12 
. Lane. a aD 
- g 7 
— 
x 
c . fsxen90 1) 8.2 - 
a 
f « 
s I é : 79°71 baw 4 
_ 7 
Ib: nolisxinsonoes es sw a8 WO Yi'b 1oi oz 
: < Cita, 7 a ~ fats 
’ . ' . ° i a. a - aclT 
Z : ‘ Ae ; Arlee a, ee tee i I ) Dan +9ilIc 70 oo Ses Jigib s 


~ “=” ‘6 —% 7 
douse }- wold msetea oat 4ne1105 brs esoped 
-— a, : a 2 : 


so90%q sft wort yilese betaluqinam amt “ed Sives 


betiunpe: esw gnimmseagomd J anolyibbs: On . formes & <a 
WR: nies ne oe 
.3.2 s7ught fe at) ey waccEs a4 al ael fom sa08 oo 4 
zs a 4 a 
| ; ; noe » , - 7 
: [ae re ws a A v 
em te 


y 
7 


aa 


5.2.2 Compensation for Feed Concentration Disturbances 


The relation between the feed concentration and the 
steam flow is a first order lag system due to the liquid 
hold-up in the first effect. Equations (5-7)..and (5-8) were 
used as a basis to derive this relationship. The mathe- 
matical development is presented in Appendix C. The result- 


ing relation in terms of perturbation variables is as follows: 





GWSE oe OWL 2. WL. - 055 eeSTRGm aU WST * CL! 
cat. HL (inde GaO LECH SCC MEEGTEG: *9CO me 
+ | STEC * CO WIo-n0 sok STEC *-WST ‘ 
| ste, 50 * 0.5 * STEC * C2 - STEC ¥ CO | * WST 
(5-9) 


It remained to assemble a digital controller that would 
produce the same output as a solution of this equation. At 
least two methods were available to arrive at this controller. 
Firstly, the equation could have been programmed and solved 
on the computer, either analytically, or numerically, and 
the result transferred to the control record concerned. In 
the analytical case, the solution is available explicitly in 
general terms. Appropriate constants are substituted by the 
program at regular intervals, and the result is calculated. 
In the numerical case, the differential equation is used as 
such in the computation. Successive approximations yield 
the value of the dependent variable at the end of a given 
interval. 


As for the second method, the direct digital con- 


tot nots wig. + aah 


- SVSD [eotdam 


Lugs o a0 rq 
goa 3 bal exgo 7 


eaB9 Iso Sromin * edd al 








| 2 


ok mated 


sd # 86 bpeaw: 
! 


+ 


gat 


¥ 
5 


a2 


trol program contained provision for exponential filtering 
of a process signal. This feature made it unnecessary to 
resort to the use of a separate, special program in the 
present case. For large time constants, the exponential 
filter approximates a first order system, and the filter 
constant is a function of the process time constant and the 
sampling interval. The mathematical definition of the ex- 
ponential filter is given in Appendix C, Table C-1. 

To demonstrate the equivalence of the two methods, 
a test run was performed during which a step change in feed 
concentration was simulated and the filter output recorded. 
The output obtained from the filter as a function of time 
agreed clearly with the solution of the given differential 
equation. Figure 5.3 shows the comparison of the two 
methods. After this test, the consideration of FORTRAN pro- 
gramming was abandoned, and the feedforward controller was 
implemented as shown on Figure 5.6. 

As the process model was linearized to arrive at 
the equation for the feedforward controller, the accuracy of 
the compensation was expected to diminish as the process 
moved away from the initially defined steady state. The 
method described provides for very simple redefinition of 
the steady state for the controller. If a new steady state 
value of the feed concentration is entered into the variable 
record from the process operator's console, a new reference 


state is defined for the controller. Control action will 
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begin as soon as the measurement moves away from this latest 
value. 

The pertinent constants for the feedforward con- 
troller were determined from equation (5-9). After Laplace 
transformation and substitution of the appropriate constants, 


this equation takes the following form: 


WST 3 = K = -0 eS 14 
CI © Taeee ae! Dass StL (5-10) 


The dynamic compensation represented by equation (5-10) was 
implemented by using one DDC loop. This loop 

(a) measured the feed concentration CI 

(b) generated the perturbation variable CI' 

(c) multiplied CI' by the gain factor K 

(d) dynamically compensated (filtered) 

(e) generated the deviation in steam flow WST' needed 

to eepeneace for CL) 

The perturbation in steam flow WST' was added to the steady 
state value in a separate steam flow control loop. The im- 


plementation is described in Appendix D. 
5.2.3 Compensation for Feed Flow Disturbances 


The feedforward compensation for feed flow distur- 
bances was determined from considerations of the level con- 
trol scheme employed. The ee atioushi between the feed flow 
and the steam flow will depend on the type of level control 
used in the first effect. If a tight level control is used, 


then the bottoms flow will respond immediately to a disturb- 


























-o. brawroibest off +02 a2uKRtehoo Inenizizeq 6dT 
anatq rot2A §, (0-2) nolbseups move bonlarnsteb ox9w xeilora— 
ststxqoyqgs ez to notjustiedue bas nolssorrotenss3 


nriot gniwollo? end esiag not nups ein ‘a 


pit. 
(or~2) eer. <a sh : ee 


Of-2) motsaups vd betnees rast nolisansgmes gimeayb edt 
goo! ahiT gool 900 smo gaieu ya beamensiqnt ~) 
soijaxzanesonos bead aiid bssuesem (3) < 

TD sidaisay sotisdruirsq end bestessneg (a). 
79398) sisg of yd ‘TO bekigttiom {o) 


(bexesiit) bstseneqmos yileoimenyd . (6) - 


hoon 'T2W wolt maete of aolteiveb sii betaveneg (s) 
ID 302 sJseneqaap oF is 


vbseta oft of bebbe eaw 'TeW wolt masse af nokgadtutzeq. sar 
ae 
ofl gool loxgnos wolt msste statage: # Ar sulav siste 


.d xibaeggA at bedisesab at notssinemslg 
: , _ 


ssonsdyuteld wolt best 20? moldaaneqned £.8.2 \ Se 


¥ Si X oe 

-tugalb wolt bsst 102 nolssensqmos. brswrokbes? ait | a: 
-fio9 {evel sit Io svodsaxeblenos, moxt beakmneseb anw f 
wolt beet eft neswied gidenoltssisz siff -boystqua ouadoe aa 
fox3noo level. to aqys eds no basqeb iitw wolt Loeadall pe: i 
















oo oa ne 
idl ME pees Apes: aigha x 21 ssonte a : 
' = i) ; is - < 
* s. Pal p a oe iY 7 a - ody a, 
iad : wer ah atboom. Aviat 7 Eye mit aa peer basic 


ast a ae >t +N. ik 






55 
ance in feed flow, and the steam flow will be related to 
the feed flow by a constant. In the case of averaging 
level control, the first effect is used as a damper, and 
the response of ene bottoms flow is exponential. This 
dampening effect introduces a first order lag between the 
feed flow and the steam flow that had to be simulated in 
the feedforward controller 

Two choices were available for the implementa- 
tion, as is shown on Figure 5.6. In the first scheme, the 
feedforward controller read the feed flow disturbance. 
First order dynamic compensation as well as a gain factor as 
defined by equation (5-11) was applied, such that the out- 
put from this controller would be an exponential curve. 

The second scheme, that based on the simplified model, 
would produce the same end result, was to base the feed- 
forward eempencatl on on the bottoms flow from the first 
effect. In this case, the controller had to be given a 
gain factor only, since the first order lag was inherent 

in the physical relationship of B1' to WI'. Both possibil- 
ities were experimentally applied. 

As it had not been attempted to describe mathe- 
matically the influence of averaging level control on the 
propagation of load disturbances through the process, the 
exact relation between steam flow and feed flow was unknown. 
In order to determine appropriate gain factors for the 


feedforward controller, constant steam economy as calculated 
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from a material balance was assumed. Starting from an 
initial steady state, the gain factor was calculated on 
the basis that the fraction of feed evaporated in the pro- 
cess be the same at the initial and at the final steady 
state. The time constant introduced between steam flow 
and feed flow by the changing hold-up in the first effect 
was assumed to be the same as that existing between the 
feed concentration and the steam flow, i.e. approximately 
equal to the average hold-up time in the first effect. 
These considerations led to the following trans- 


fer functions (10): 
WST 3 CsyeTo* 0.40 (5-11) 
WL Ss al Bo 8 Loel 
WST ts) (5-12) 
Bis Cs 


The constants were supplied to the controller in the same 


Teak © 


manner as in the case of feed concentration. As can be 
seen on Figure 5.6, the flexibility of the digital con- 
troller provided the choice of using either one of the 


two methods. 
5.3 Inferential Controller 
5.3.1 Process Simulation 


Inferential control consists of controlling a 
process variable whose value is not measured, but inferred. 


In this study, it has been applied to the concentration of 
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the product flow from the evaporator. A mathematical model 
of the process was developed and programmed for the computer. 
System routines were available which permitted communication 
between user programs and the process variable table. For 
the simulation program, the process model given by equations 
(5-7) and (5-8) was linearized using perturbation techniques. 
The linearized version was then Laplace transformed such that 
the model in terms of perturbation variables could be written 


in the following form: 


a CSTIC oe CSTOW — 
C2(a) = 3- cstip. * Cl(s)t-Sestgq7 * WICs) 
& ee rn * WST (s) (5-13) 
Or esTacu. CSTTW 
EOsGye 2s te rGsTan C28) ab BececsTing = WL)Cs) 
CST5W 
+ —sorsp * WST'(s) (5-14) 


It was found convenient to represent this process 
model in block diagram form as shown on Figure 5.4. Each of 
the six boxes on this diagram represents one term of equation 
(Se) 3)0r-( 5-14). 

Equations (5-13) and (5-14) were solved analytically, 
term by term. A program was then written to evaluate C2 and 
CO by substituting appropriate values into the analytical 


solution of these equations. 
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5.3.2 Basic Implementation 


The constants required for the linearized model were 
calculated using steady state values obtained from the process 
variable table at the start of each run. All programs are 
listed in Appendix E. The constants were calculated by the 
program DIANE, and the solutions of the equations were eval- 
uated by the program COLOR. This latter program was called 
into core from disk storage twice during the specified control 
interval by a repetitive calling routine. It performed two 
different functions on an alternate basis. At the beginning 
and at the end of each control interval, the values of pro- 
cess input variableswere obtained ernie process variable 
table. At the midpoint of each interval, the product concen- 
tration was calculated from the model and its value trans- 
ferred into the meas neement of the concentration control 
loop. 

Considering the control interval from time nT to 
(n+1)T, the measurement of the process inputs was done at 
time nT. At time (nt%)T the program calculated the value of 
the product concentration, CO((nt1)T) based on the inputs 
measured at time nT, which was sent to the DDC control pro- 
gram. This was done to spread the computing load over the 
full interval and partially compensate for the delay due 
to sampling. 

In the analytical solution of the differential 


equations, the value of the time was constant and equal to 
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the interval between two successive calls to the program. 
Thus, the result obtained at every call represented an exact 
solution of the equations at this particular instant. AI1l 
constants as well as the process vector consisting of inputs, 
outputs, and intermediate quantities on Figure 5.4, were 
written into a special section of core storage called the 
skeleton common area. The vector was updated at every call 
to the program, and the latest values used for calculations 
at the next calling instant. 

As no program interrupt levels were available at 
the time the experimental work began,queuing routines were 
used to bring the program into core for execution. Thus, 
although the program was assigned the highest possible pri- 
ority in the queue, it could not interrupt other process pro- 
grams executing at the time of the call. The length of the 
delays caused was not known, but it was not expected that it 
frequently exceeded the time lapse between two successive 
calls: 


5.Jde3. Brogram-Structure 


The information flow within the inferential control 
system is shown on the logical flow diagram of Figure 5.5. 
The program consisted of a number of interconnected core 
loads. The functions performed by these core loads are 
summarized in Table 5.1. There are several features of this 
system that merit further discussion. 


Firstly, the division into parts made it easy to 
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Figure 5.5 Information Flow Diagram 
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Functions 


Initiated by operator 


1s 


Operator enters flag l 


a) immediate exit 

b) calculate new constants for inferential 
controller 

c) use previous constants 


Operator enters flag 2 


a) continue by calling DIANE 
b) turn off timer and terminate program 


Initiated by INFER 


Obtains parameters from POC 


Obtains process inputs from PVT 


Obtains present value of product concentration 
from PVT 
Calculates constants for process model 


Calculates initial process vector 
Operator enters flag 3 
a call EDITH 


b) exit 
c) repeat execution of DIANE 


Initiated by DIANE or INFER 


Operator enters flag 4 to select model to 
be used 
Flag 5 is set to read input values 


Requests time delay between process and model 
from POC 


. Initializes timer to the desired calling 


interval for COLOR 


Table 5.1 Program Summary 
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Programs Functions 
COLOR Initiated by timer 


Timer is automatically reset by the timer 
routine 


1. Obtains present values of process inputs 
or (according to flag 5) 


Simulates process to obtain product concentra- 
tion 


2. Writes value of product concentration into a 
disk file 


or (according to flags 5 and 6) 


Transfers value of product concentration 
to feedback controller 


3. Resets flag 5 


Table 5.1 continued 
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modify the program. Any section could be removed and replaced 
by a different version without affecting the other sections. 

Secondly, provision was made for the use of differ- 
ent process models. An alternate model could be called 
simply by changing a logical flag number. It is expected 
that this feature will be made use of in future work with 
this equipment. 

Thirdly, multiple entry points to the timer called 
program COLOR were provided. This feature allowed for 
different functions to be performed at different calls from 
the timer. 

Fourthly, the most important feature was the possi- 
bility of repetitive execution at fixed intervals. Once 
the timer was initialized, it queued the simulation program 
at specified intervals until it was turned off. 

Fifthly, the program COLOR would write the results 
of the simulation into a file on disk from where they would 
be transferred as product concentration measurements to the 
feedback controller. Obviously, the value so transferred 
did not have to be the last value in the file, but could be 
specified to be any value from the file history. It was 
therefore possible to base control action on past information 
obtained from the model. This feature was applied in ex- 
periments INFI2 and INFII2. The model had ‘been found to 
respond to disturbances seven minutes faster than the process. 


In an attempt to cancel this lead, the controller was supplied 
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with measurements seven minutes after they had been calcu- 


lated by the model. 


5.4 Options for Control Action 


The controller shown on Figure 5.6 was structured 
in such a way that any one component could be used individ- 
ually or in combination with one or more of the other com- 
ponents. The inferential control system essentially substi- 
tuted a mathematical model for the process. that would supply 
the measurement for the feedback controller. Either the 
real process measurement or the inferred measurement could 
be used for control purposes, and feedforward compensation 
could be applied in either case. Moreover, the inferential 
control program shown on Figure 5.5 could be executed without 
actually performing control. In this case, the calculated 
product concentration was written into a disk file for plot- 
ting purposes, but was not transferred to the controller. 
This scheme permitted comparisons between model and process 
behaviour at times when the model was not involved in con- 
Erol 

These features made the controller an extremely 
versatile and well functioning unit. Reference is made to 


Appendix D for details of controller implementation. 


5.5 Direct Digital Control Arithmetic 


The operation of the direct digital control pro- 
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gram and the nature of variable records in this system are 
such that arithmetic and numerical operations may be per- 
formed by judiciously interconnecting a number of variable 
records. The method of using an exponential filter algo- 
rithm in the implementation of feedforward control was one 
example of this use of direct digital control operations. 
In general, "standard" variable records may be used to per- 
form several operations, such as multiplication, division, 
addition, subtraction, or numerical integration. 

The combination of actual control action with 
this type of manipulation is an important asset of direct 
digital control systems. It permits the implementation of 
feedforward and other advanced control actions without 
additional programming. The efficiency of this method of 
manipulation depends on the type of application and on the 
number of variable records involved. 

The following is an illustration of several 
possible manipulations: 

1. Proportional control algorithm 


output = Kp * meas, -_setpt. + bias 
p 


a) subtraction: bias = 0, Kp =l1 

b) addition: Kp = 1, setpt = 0 

c) multiplication: setpt = 0, bias = 0 
2. Proportional - integral control algorithm 

output = Kp * [ meas. - setpt. | En 

Kp * KI * [tieas. - setpt | + REST 




























4 
\y a : 
. i 
rae S. 
i ' 7 a 
malt ee 
aa : 


ca 
; Py a 

gsacmeteye aids at ebtoset oidatzev do etu3en ‘ody bas aesg ce 
-yaq ad yam enoltazsqo feo baonun bus oijamisits. and4 cous 
e{dsiisy to sedmip 8 gntjosnnolTeIoF vigvoLoibut yd borzot 


1t3 Tebteeneqxe os gateu: to boritem odT .ebx039% 
lo poLsstnametqad ofa ot mitts yi 


a 


a 


ad 


ocala 2a: 
ono asw Lozvanes braewrdrbes? 
fovtnoo Lavigth gootkb. lo-esu eins 30 olqmaxe 
abyooot ofdaiirsv ‘hesbnete” ,lareasg ot 


itfom ag dove-, eootse7eqo fasavee arroy 


qoletvih .noksaorrg. 


Enos 3 Sf SQ 
ra9 03 beav 80 Yaa 


fentroauna 10 . nokzon7adye , no bs tbbs 


,mozsserav: 


losdeo> {agdon to notsanidaes eft 


a 


7 ny 
is IW HOLS 2s 
< 


roectoymt ms et aobtaduqhesm to sqva ekdy 


2o97T4D 10 49 
26 noltsanarolgt et@ astnssq: +t .,ameteya forsnos Ist igid 


foxgnos beonavbs t9Af30 bos ‘brawrotbest 
a 


ty0hs Iw ea°oIsos 


to yanetoritis oT _ grbmergoxg Tano ps thbs 
to eqy? eit 16 ebrreqsd colssluqiiem 


to borizo0 


sf3 no bos coltesligde 
boviovat ebroost aidpizav to sedeud 
fevoyves 26 noksertaplib me at antwol lot sit coe 


aE 
mitttosiIs laxtnoo Ianokrrogost “yi gles - 
‘enld + [ .tqee - pees * a ~ suqsuo : ar 
ay ~~ * = ae ie 


{ = qa ,0 = enid -nokioatadue ts 

; 0 = tqg08 .i = qu :notsibbs- {a 

0» sald ,0 = aq20 -aotsaoiiqislye ‘eo 
aris tiogis - foxanog. Sergetat + 3 


pe a 


in >», 


re it | aie! 
+ Pili 





an 


68 


where REST is equal to the sum of Kp * KI *|meas 
- setpt | evaluated at all previous control 
intervals. 

a) same manipulation as in l. 

b) integration: setpt = 0, KI = sampling 
interval delta t, Kp = O to eliminate the 
proportional term (a value of Kp = 1 is 
then automatically selected for the integral 
term) 
output = x meas * delta t 


Proportional - derivative control algorithm 


output = constant “| wmeas) (filtered) + 
t 


proportional action as per 1. ] 
filtered derivative action: setpt = 0, 


Kp = 1 or other values 


output = —a(meas) (filtered) 


Averaging control algorithm 

Output = Y word (i) of loop (Jj), for all speci- 
fied pairs of (i,j) values 

Different poll and phase times can be used to 
perform iterative or sequential calculations. 

In the early stages of this study, the differ- 
ential equation defining the action of the feed- 
forward controller was solved by a group of 


six variable records that were arranged in such 
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a way as to perform a numerical integration. 

6. Minor modifications to the data accumulation 
program would provide for implementation of 
pure time delays. 


7. Model simulation 


A model of the form K * ih 
cs +1 Ss 


could be implemented by two DDC loops. The 
first loop would perform the multiplication by 
the gain factor and the filtering action equi- 
valent to the given first order system. The 
second loop would perform the integration. 

8. Other feedforward models 
The simulation of transfer functions of the 


type eo A eerrecy Ge ae 
Cizg st T3 sot 


is possible with the use of DDC loops. The 
efficiency of this method would have to be 
evaluated for a given application. 

These examples serve to illustrate further the 
versatility of the direct digital control system and indicate 
that DDC programs can be used for other applications where 
the dynamic compensation or control scheme differ from those 


implemented in this study. 
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6 EXPERIMENTAL PROGRAM 


6.1 Summary of Experiments 


The complete experimental program is shown on 
Tables 6.1 and 6.2, and the complete results obtained from 
each experiment are given in Appendix A. 

To ensure that all experiments were performed under 
Similar process conditions, a steady state of reference was 
established. This reference state is defined in Appendix 
A, All load disturbances introduced into the process were 
20% of the steady state value in magnitude, and used this 
Steady state as a pivot. 

The program was planned such that it would demon- 
strate the performance of all parts of the digital product 
concentration controller. These included the feedback con- 
troller responding to both measured and inferred signals, as 
well as the feedforward compensation for feed flow and con- 
centration disturbances. 

The feedforward compensation for flow disturbances 
was based on the feed flow signal in all experiments labelled 
FF II and FFB II, and on the first effect bottoms flow signal 
in experiments labelled FF IV. 

Of the experiments using the inferential controller, 
runs INF I2 and INF II2 used a process model containing a 
time delay. For runs INF I and INF II, no time delay was 


implemented. 
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6.2 Experimental Procedure 


All experiments listed on Tables 6.1 and 6.2 were 
conducted in three extended periods of operation. At the 
beginning of each of these periods, the process was brought 
to the same steady state. All disturbances were introduced 
from the process operator's console. 

Data accumulation records were available to retain 
measurements of the most important process variables at speci- 
fied intervals. These included the product concentration, 
the steam flow, the liquid levels in both effects, and the 
bottoms flows from both effects. From the accumulation 
records } data were collected in a disk file for later plott- 
ing. At the end of each run, all required diagrams were pro- 
duced directly on the system plotter. The disk file was then 
reset for the next run. 

At the beginning and at the end of each run, a 
material balance was performed around the major sections of 
the process, that determined the steam economy and the errors 
of closure. The results of this balance could only approxi- 
mate the condition of the process at that particular in- 
Stant, as the measurements used in the balance calculations 
were instantaneous values that could contain noise compon- 
ents. For the later part of the experimental period, a more 
sophisticated material and energy balance program developed 
elsewhere in this department was available. This program 


obtained several consecutive measurements of each variable 
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and used an average in the subsequent calculations. It also 
contained the option to adjust the process vector such that 
the error of closure disappeared. 

The refractometer recording the product concentra- 
tion was not reliable at all times. To produce a continuous 
chart output that was correct within the error of measure- 
ment, samples were taken manually at intervals of approxi- 
mately seven minutes and analysed on a bench refractometer. 
Whenever the concentration obtained from this analysis 
differed from that shown on the recording chart, the instru- 
ment was adjusted on the spot so that it would indicate the 
correct value. The discrete steps visible on some of the res- 
ponse curves in Appendix A are a result of this adjusting 


procedure, 


6.8 Supplementary Software 


To permit the proper functioning of the procedure 
outlined in the previous section, several additional programs 
had to be written, They have been kept as general as poss- 
ible and are expected to be used in future work with this 
equipment. At the end of each run, the program BEGIN was 
called to reset the disk file that collected data from accu- 
mulation records. The diagrams at the end of each run were 
produced using programs THIS and THAT, The material balances 
were performed by FDFWD., Listings of these programs are cor 


tained in Appendix E. 
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Feed Flow Rate 
1bs/min 
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7 DISCUSSION OF RESULTS 


7.1 Experimental Data 


7.1.1 Scope of Experimentation 


As has been stated at the outset of this report, 
it was assumed that the feedstream to the evaporator was 
not available for manipulation. It was the function of the 
control system to correct for disturbances in the feed 
stream by manipulating the steam flow. 

The experimental program was planned to permit the 
evaluation of the following digital control systems: 

basic feedback control 

feedforward control 

combined feedback and feedforward control 

inferential control, 

Both the flow rate and the concentration of the 
feed stream were disturbed by applying positive as well as 
negative step changes. 

The results obtained have been summarized on 
Figures 7.1 to 7.4. For all experiments, the initial and 
final steady state operating conditions were maintained as 
close as practicable to the design conditions, and all step 
changes introduced were of the ae magnitude. Therefore, 
direct comparisons of response curves are possible. All 
data obtained during the experimental program, including 


initial and final steady state conditions and transient res- 
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Response to a Step up in Feed Concentration 
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Response to a Step down in Feed Concentration 


LEAT: 8.09 NDI a cet RC ny SE Ee | | 


OFT 


PO 


ee ee ee ER Te 


— 
© 
n 


00 


product concentration in per cent 
e = 
: ° 
A ios repost rtepe mente temumiieerriiigsl aie teases teamed 
ye | 
| 
| 
| 
| 
| 
| 
| 
ae am aa 





SS cnEIEEEIEE ESERESESSERERIEEEE euemememeeen 


| Py et = 0 We ge ia eae 
0 20 40 60 80 100 


time in minutes 


Figure 7.2 


to 


ta 
Bn 
ag 
2 
» 
oO 
8 
Q 


brogyee 





78 


Response to a Step up in Feed Flow Rate 
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ponse curves for all runs, are included in Appendix A. 

All runs were terminated as soon as a new steady 
state was reached. As is apparent from the diagrams, the 
feedforward runs required considerably less time to reach 
steady state than the feedback runs. 

Direct comparisons of the present results with 
those obtained by Andre (2) and Wilson (31) are not possible, 
due to physical changes, different steady state conditions, 
and differences in disturbances applied. General agreement 
can be shown in most cases. 

Figures 7.5 and 7.6 show the results of two open 
loop experiments. As had been anticipated, the model res- 
ponded faster than the process. Andre had reported a simi- 
lar situation which he attributed to the imperfectly mixed 
system. In the present case, two additional reasons are 
suggested for the discrepancy. Neither the delay in the 
feed piping nor the lag created by the averaging level con- 
trol system were considered by the model. How well this 
simple model was suited for control purposes was determined 


by the experiments using the inferential controller. 
7.1.2 Feedback Control 


The general shape of the response curves obtained 
from the feedback runs (FB) shows evidence of the tuning 
criterion employed. The feedback controller was tuned for 


rapid damping and gave rise to analogous response curves : 
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83 
for all feedback runs. 

The maximum deviation is larger for step changes 
resulting in a rising product concentration, as can be seen 
on Figures 7.1 and 7.4. This result is due to the varying 
sensitivity of product concentration to steam flow. As the 
concentration rises, as in runs FBI2 and FBII3, the change 
of concentration per pound of water evaporated increases. 
The difference in maximum deviation between runs FBII3 and 
FBII2 is much more pronounced than between runs FBI2 and 
FBI3, because the averaging level control system employed 
gave rise to varying liquid levels in the case of flow dis- 
turbances. A step down in feed flow causes the liquid hold- 
up to decrease, which entrails a larger change in concentra- 
tion per pound of water removed. 

The time required to regain steady state follow- 
ing a concentration disturbance was less than following a 
flow disturbance. This difference is due to the length of 
time taken by the liquid levels to return to their steady 
state values after a flow change. 

The experimental results reported by Wilson show 
that in the case of a concentration change, cycling lasted 
for 180 minutes, and in the case of a flow change, steady 
state was not regained after 200 minutes. The higher through- 
put rates used in the present experimentation are considered 
a major factor in reducing the response time, as they 
decrease the hold-up time of the process. In addition, with 


the controller constants used here, cycling was not exper- 
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84 
ienced. Andre did not report on concentration disturbances. 
In general, the shapes of his response curves for flow dis- 


turbances are similar to those obtained in the present work. 
7.1.3 Feedforward and Feedback-Feedforward 


The improvements in product quality control for 
this evaporator that can be achieved by the use of feed- 
forward compensation had been indicated by Wilson. In the 
present work, it has been attempted to move toward general 
purpose software in the implementation of feedforward control. 
As the present experimental work confirmed the results obtain- 
ed by Wilson, it is felt that as feedforward compensation is 
possible with standard software, it will become increasingly 
common in industrial processes. 

The finan improvement due to this type of con- 
trol can be seen fram Figures 7.1 to 7.4. Concentration 
disturbances were compensated for completely by the feedfor- 
ward controller. As no correction for off-set was required, 
the combined feedback-feedforward controller produced the 
same product quality as the feedforward controller. In the 
case of flow disturbances, the compensation was not complete, 
such that a slight off-set in product quality occurred that 
was eliminated by feedback correction for both positive and 
negative step changes. 

Two methods of providing feedforward compensation 


for flow changes were applied. For runs FFII2 and FFII3, 
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85 
the steam flow adjustment was determined by the relation 
WST (s) = (Kl / (ts + 1))* WI, (s) 
For runs FF IV1 and FF IV2, the relation 
WST (s) = K2 * Bl (s) 


was used, where 


WST = steam flow 
WL = feed flow 
Bl = first effect bottoms flow 


The response curves for runs FFII2 and FFII3 are 
analogous. The off-set occurring in these runs was due to 
an insufficient gain. If the constant Kl, that had been found 
from theoretical considerations, were increased by an appro- 
priate amount, the off-set would disappear. 

As for runs FFIV1 and FFIV2 the steam flow was 
directly proportional to the first effect bottoms flow, 
Figures 7.7 and 7.8 have been included here for illustrative 
purposes. They show the behaviour of the bottoms flow for 
these two experiments. A careful examination of these fi- 
gures plus recordings of the other variables did not pro- 
vide any basis for the observed difference in response 
shown on Figures 7.3 and 7.4. Therefore, it must be conce- 
ded that the cause of the off-set occurring in experiment 
FFIV2 is unknown. 

Small variations in process behaviour such as these 
suggest the continuing need for feedback correction in addi- 
tion to feedforward compensation. This point is illustrated 


very clearly by experiments FFBII2 and SSUES Me which feed- 
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back correction eliminated the off-set. 

The experiments have shown that both methods of 
feedforward compensation are effective and are valid alter- 
natives. The combined feedback-feedforward controller has 
made the process under study effectively immune to disturb- 
ances in feed conditions, and there appears to be no doubt 
that the method of feedforward compensation developed in 
this study can find applications in a variety of similar 
processes. 

The feedforward compensation achieved by Wilson 
with analog components did not appear to produce results 
as satisfactory as the digital feedforward controller used 
here. All feedforward runs reported by Wilson contain fluc- 


tuations in the product response. 
7.1.4 Inferential Control 


The inferential controller has been used in a 
series of experiments labelled INF on Figures 7.2 and 7.3. 
This type of control is employed in situations where it is 
desirable to control a process variable without physically 
measuring it. It was the intention of this work to apply 
the concept to a process that was more complex than those 
described in the literature (24, 14). The results obtained 
yielded more information about the usefulness and flexibili- 
ty of inferential control than had been anticipated. 


It had been observed from open loop experiments 
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shown on Figures 7.5 and 7.6 that the model responded 

seven to eight minutes faster than the process. On imple- 
menting the inferential controller, provisions had been 
been made for arbitrary selection of the length of time that 
would elapse before the information obtained from the model 
was supplied to the controller. The four experiments of 

the series INF clearly show the effect of varying the time 
correspondence between model and process. 

In runs INFI2 and INFII2, the controller obtained 
its information from a model that contained a time delay term 
which improved the agreement between the open loop response 
curves of the process and the model. In runs INFI and INFII, 
the information came from a model that did not contain the 
time delay and thus responded seven minutes faster than the 
process. The influence of this time difference is apparent 
from the diagrams. If allowance is made for the slight off- 
set that actually did occur in the case of INFII, it is 
apparent that the undelayed model returned the process to its 
original steady state, whereas the delayed model did not. 
This result was not expected, nor can it be explained theor- 
etically, since a time delay usually affects the speed of 
response of a given system, but not the final accuracy of the 
response curve when compared to the setpoint. Since the 
delay applied to the model was the only known difference 
between experiments INFI and INFII on the one hand and INFI2 


and INFII2 on the other, no satisfactory explanation is 
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available for the off-set that resulted. As far as the time 
delay itself is concerned, it is natural to raise the question 
of the optimal time difference between model and process res- 
ponses. It is felt that the limited number of experiments 
using the inferential controller have helped illustrate the 
importance of this question, rather than answer it. In this 
particular case, the controller obtained its information from 
the simplified model which did not contain a time delay tern, 
but experimental runs had shown that it anticipated the pro- 
cess by approximately seven minutes. It is hoped that the 
results obtained will generate further work on this aspect 

of inferential control. 

Returning to experiment INFII, it is felt that the 
slight off-set that occurred at the end of the run is a con- 
sequence of the model limitations. A cross reference to 
Figure A-16 of Appendix A shows that the model had returned 
to steady state. As is seen from the open loop run on Figure 
7.6, the model did not exactly reproduce the process response. 
On programming the controller, the process model has been 
included as a replaceable unit, such that in future work more 
complete models can be substituted. 

The question obviously will arise as to the rela- 
tion between the accuracy of a model when compared to an 
open loop process response, and the quality of control that 
can be achieved with it in an inferential control system. 


One phenomenon that has become apparent in this work is that 
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although the model response anticipated the process and thus 
did not reproduce the open loop process response curve, the 
inferential control runs INFI and INFII resulted in better 
product quality control than the conventional feedback 
runs FBI3 and FBII2. The lead of the model with respect to 
the process is suggested as explanation for this phenomenon. 

It is apparent that in an inferential control system, 
model parameters become in a sense control parameters, and 
it does not follow that a model that exactly reproduces the 
open loop process response is also the best one for control 
purposes. 

On one occasion during the experimental program, 
the process was operated under inferential control for a 
period in excess of four hours. During this period, various 
disturbances were introduced into the process. The results 
are shown on Figure 7.9. The arrows indicate the instances 
at which step changes were applied. For the first run, infer- 
ential feedback alone was used. At the end of this run, feed- 
forward compensation was added. The step changes were the 
following: 

1. feed flow 5.5 to 5.0 lbs/min and feed concen- 

tration 3 .07toea.s7 

wersecpoint 7,6 to 9.07 
. feed concentration 3.3 to 3.6% 


feed flow 5.0 to 5.3 1lbs/min 
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At the beginning of the first run, the refracto- 
meter did not function, as can be seen on the diagram. As 
soon as this was noticed, it was corrected, which explains 
the sudden steep rise of the process response curve. The 
malfunctioning did not have any effect on the operation in 
this case, as the model output was used for control. The 
process output was only recorded. The occasional deviations 
of the model from the process in the later part of the ex- 
tended run are all within the error of measurement. 

This experiment was intended to test the perform- 
ance of the combined controller under a set of operating 
conditions that was realistic in the general industrial sense; 
namely, an extended run with frequent load disturbances, 
without necessarily reaching steady state between two succ- 
essive changes. It is felt that the controller performed 
satisfactorily, which admits the belief that the control 
system developed is generally applicable to processes of a 


similar nature. 


7.2 Equipment Performance 


The process equipment performed satisfactorily 
throughout the entire program. Andre had reported that 
frequent switches of feed tanks were necessary which intro- 
duced undesired disturbances into the process. With chen 
ratio controlled feed system used in this study, this type 


of problem did not arise. The product was pumped back into 
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the feed tank and feed concentration monitored by the ratio 
controller, such that the operation was continuous. 

Wilson had faced the problem of sticky valves on 
various occasions due to the properties of the sugar solu- 
tions he used. The triethylene glycol used in this work did 
not present difficulties of this nature. 

The higher throughput rates reduced the hold-up 
time, such that generally, the process regained steady state 
faster than was the case with Andre and Wilson. 

The steam economy generally achieved in industrial 
double effect evaporators is 1.6 (16). In this work, it ran 
from 1.5 to 1.7 which is one indication of realistic operat- 
ing conditions. 

The effect of the averaging level control is best 
illustrated by Figures 7.7 and 7.8. Due to the dampened 
response of the bottoms flows, disturbances propagated only 
gradually and control of product quality was facilitated. 

The errors of closure for all material balances 
are summarized on Table 7.1. Errors are expressed as per- 
centage of inlet flows. All calculations were carried out 
using raw data obtained from the process at the time of the 
experiments concerned. As a rule, the steady state data at 
the start and at the end of all experiments were provided by 
a list of instantaneous readings that was printed by a she 
gram written especially for this purpose. Only in a very 


few cases were the noise components contained in these in- 
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Table 7.1 Percentage Errors in Material Balance Calculations 
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stantaneous readings so obvious that it was necessary to 
resort to the Foxboro recording charts for a more repre- 
sentative steady state value. This applied mainly to the 
condensate flow from the condenser. No control was exer- 
cised over this variable, and fluctuations were quite 
severe. A typical recording is shown on Figure 7.10. 

The errors of closure for the overall balance 
were consistently positive, which must be attributed to 
slight shifts in instrument calibrations. 

Table 7.2 shows the results of energy balance 
calculations. The heat transfer coefficients obtained for 
the first effect are considerably higher than those reported 
in previous studies. This is attributed to the higher 
throughput rates which gave rise to an increased heat 
flux and higher induced circulation velocities of the boil- 
ing liquid in the heating tubes of the first effect. The 
higher velocities are considered to be the main factor res- 
ponsible for the increased heat transfer coefficients. In 
addition, the higher throughput rates resulted in pressures 
between 4 and 7 psig in the first effect. Andre and Wilson 
had operated at a slight vacuum. The boiling point of the 
solution in the first effect reported by Wilson was 190°F. 
against 220°F. for the present study. The higher boiling 
point resulted in a lower liquid viscosity which in turn > 
gave rise to higher heat transfer coefficients (3). The 
average coefficient for the second effect was slightly 


lower than that reported by Wilson, which is attributed 
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to reduced circulation velocities that were necessary to 
obtain satisfactory level readings in the separator. A 
comparison of the heat transfer coefficients reported in 
various studies with this equipment is shown in Table 7.3. 
The spread of the coefficient values for the 
first effect is larger than for the second effect. The 
experiments for which coefficients of 900 and over have 
been obtained were all conducted during a period when the 
vacuum service in the building did not supply the suction 
required for this process. It was necessary to lower the 
setpoint of the vacuum in order to maintain good control. 
This vacuum change carried through to the pressure in the 
first effect, amplified by a factor of two. The resulting 
boiling point changes were thus more pronounced in the 
first than in the second effect, such that the change in 
coefficient values was noticeable in the former case and 
not in the latter. All energy balances were closed by con- 
sidering the heat losses to the surroundings, calculated by 
taking the difference between heat input and heat output. 
Average values of heat losses were as follows: 
Total heat input 4960 BIU/min 
Losses: First effect 208 " " 
second effect! «179.0 ©." 
interstage 
piping 129 vt'fart! 
condenser Ai Aaash ot 


Total losses 559 BIU/min or 11% of total heat 


input. 
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7.3 Direct Digital Control Performance 


At the outset of this study it was mentioned that 
presently only 25% of digital computers installed in the 
chemical industry perform direct digital control. It is 
hoped that the experimental work presented here will con- 
tribute its modest share to demonstrate the versatility of 
direct digital control as applied to continuous processes. 

In the present study, various complex configura- 
tions, such as cascade loops and ratio control loops, were 
applied successfully to the pilot plant scale evaporator. 
The advantages of these configurations have been demonstra- 
ted in the feed system and in the level controls used. 
There is no doubt that they can be applied effectively in 
other industrial installations as well. 

The piney ena etee operating system of the com- 
puter allowed for the process to be controlled independent- 
ly from other computer users. With the exception of some 
problems in the intial phase of the program, no difficul- 
ties were encountered as a result of the time-shared opera- 
tion. 

Since all measurements and setpoints could be 
made available in engineering units from the process 
operator's console, a steady state of reference could easi- 
ly be defined and achieved, and load changes could be 
introduced with ease and high accuracy. 


Although many aspects of direct digital control 
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were considered in this work, its capabilities have not 
been completely made use of. Digital controller tuning 
has not been employed in its true sense, and the possi- 
bilities offered by digital filters have only been ex- 
ploited in a few specific cases. 

Last but not least, it is expected that the 
present work has contributed to the development of stan- 
dard software for direct digital control applications. 
The methods used to implement the feedforward as well as 
the inferential control system are not restricted to 


this particular equipment, but can be employed generally. 
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8 CONCLUSIONS 


1. The advantages of a digital computer over conventional 
analog control instruments as illustrated by this work 
include the following: 


Control configurations were easily modified. 


Controller constants were implemented with high 
accuracy and changed readily. A practically un- 


limited range of constants was available. 


User written programs communicated with the direct 
digital control program. This feature allowed for 


the application of special algorithms and methods. 


A wide range of control schemes was implemented 
by interconnection of individual control loops. 

2. A method of implementing feedforward control using 
"standard" DDC programs was demonstrated. This specific 
method was only one out of several possibilities offered 
by this system. The standard DDC software permits accu- 
rate implementation of various kinds of feedforward models 
as well as other arithmetic and numerical manipulations. 
In the case of the evaporator, feedforward compensation 
was based optionally on feed flow or on first effect 
bottoms flow signals. The two schemes were equally effec- 
tive and were shown to be valid alternatives. The digi- 


tal feedforward controller essentially prevented load 
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disturbances from affecting the product quality, but 
continued need for feedback correction was indicated 

to correct for off-sets and undetected or unknown 
disturbances. 

The physical modifications of the equipment, the averag- 
ing level control scheme employed, as well as more suit- 
able controller constants contributed to the achieve- 
ment of improved control under the basic feedback scheme 
as compared to previous studies. 

With the aid of the digital computer, inferential con- 
trol was applied to the evaporation process using a 
very simple mathematical model. The performance of the 
inferential controller was influenced significantly by 
applying artificial leads or lags to the model with res- 


pect to the process. 
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NOMENCLATURE 
ADC Analog-to-digital converter 
Bl First effect bottoms flow 
B2 Second effect bottoms (product) flow 
CI Feed concentration 
co Product concentration 
cos Current output station 
CPU Central processing unit 
Cl Concentration in first effect 
C2 | Concentration at inlet to second effect 
DAC Digital-to-analog converter 
DDC Direct digital control 
DDC Loop Process variable record 
delta Denotes an increment 
H Computation interval 
Hl Liquid hold-up in first effect 
H2 Liquid hold-up in second effect 
1/0 Input -output 
K Constant 
KI Integral control constant 
KP Proportional control constant 
ma Milli-amperes 
meas Measurement 
MPX Multiplexer 


neg Negative 


s i 
aca & - 


f 


t9So0n09 


- 
= 
rr 


taiia joer EC ; OG 


IBATSYV egoo0sd qoo.l oad 


ns eszoosd 


=, %,° 
- te 


BI gGmM ov 


qjvo-suqn z - 


+ 


Jneseaco: 


Sa > ey hgh 


tenon Iortnes large? atk 


lanotd Pp ere 
eo soqmu= bit [ i 


ee 


LT7JNOO 


OAR 


szisb 












Ol 
02 
PD 
PI 
POC 


poll time 


phase time 


PVT 


ring buffer 


S 
setpt 
ss 
STEC 
c 
TEG 
TSX 


V-core 


First effect overhead flow 
Second effect overhead flow 
Proportional-derivative 
Proportional-integral 
Process operator's console 


Time interval between two successive 
calls to a variable record 


Time elapsed between the beginning of 
a poll time interval and the call to 
the variable record, 

phase time < poll time 

Process variable table 

Data accumulation record 

Laplace operator 

setpoint 

Steady state 

Steam economy 

Time 

Triethylene glycol 

Time-sharing executive system 
Variable core 

Feed flow rate 

Steam flow rate 

Time constant 


A prime denotes a perturbation 
variable 


A horizontal bar denotes a steady 
state value 
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Appendix A 


Steady State of Reference 


Sample Material and 
Energy Balance 


Steady State and Transient 
Data for All Experiments 


Transient Steam Flow 
Data 


A- 2 

A- 3 to A- 6 
A- 7 to A-49 
A-50 to A-56 
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Steady State of Reference 


Reference is made’ to Figures /A-1 B-19 3.1,jand 3.2 for the 


Vari- Vari- 
able able 
Name Value 
Flow Rates 
Fl 2.00 
F2 3.30 
F5 1A a8) 
F6 1.66 
F7 1.64 
F8 5.00 
F9 48.00 
F10 190.0 
F10 Ce 
Concentrations 
CL eV 
C6 9.65 
Pressures 
P20 5.9 
P22 -15.0 
Levels 
ree 9.10 
L14 16.00 
Temperatures 
LL 106.6 
T2 2215.0 
T4 180.7 
dk 246 <1 
T1LO 220.8 
AWeye 123.0 
Le? 156.3 
LS i Oe k 
T19 2249 
T28 Me 92 Py 
T29 65.0 
T34 Pa 5 
diy 195.0 


lbs/min 
1lbs/min 
1lbs/min 
1lbs/min 
1lbs/min 
1lbs/min 
lbs/min 
1lbs/min 
ft/sec 


percent 
percent 


psig 
in Hg 


in H20 
in H20 


deg. 
deg. 


deg. 
deg. 
deg. 
deg. 
deg. 
deg. 
deg. 
deg. 
deg. 
deg. 


Pe bf af ef be bf 


location of measuring devices 


Variable Description 


steam flow to first effect 
first effect bottoms flow 
first effect overhead flow 
product flow 

second effect overhead flow 
total feed flow 

cooling water to condenser 
second effect circulation rate 
circulation velocity 


feed concentration 
product concentration 


first effect pressure 
second effect pressure 


separator liquid level 
first effect liquid level 


condenser cooling water out 
first effect vapour 

solution to second effect 
steam condensate first effect 
steam to second effect 
condenser condensate 
separator vapour 

steam supply to evaporator 
liquid in first effect 

steam condensate second effect 
condenser cooling water in 
product from second effect 


feed to first effect 


Table A-l1 
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Sample Material Balance, Experiment FBI2 


Section Considered 


Total: Process 


Input F8 

Output F5 L753 
F6 jreags 
F7 ek 


Solute 
Input F8*Cl 
Output F6*C6 


First Effect 
Input F8 
Output F2 3.14 


F5 e7s 


Second Effect 
Input F2 
Output F6 LoL 


F/ Lleol 


Steam Economy 


(F5+F6) /F1 


Time Basis 1 Minute 
Mass Units Lbs 


Start Error 


° 


5.02 


4,75 Byres 


0.140 
0.136 258 


Table A-2 


1.46 
T69 


1.34 


End 


78 


09 


.L69 
.170 


98 


es 


46 


a ae: 


.68 


Error 
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Sample Heat Balance, Experiment FBI2 


Time Basis 1 Minute 
Heat Units BTU 


Section Considered Stream Start Loss End Loss 
rirse ertect Sinput. Fleat TL5 2250 1973 
F8 at T7 792 794 
3042 AM Ee) 
Output Fl at T5 384 334 
F2eate lug 580 640 
BotatFi2 2000 1680 
2964 782 2654 118 
Second Effect Input F5 at T10 2000 1680 
cabins Tec gal F2 at T4 468 508 
2468 2188 
Output F5 at T28 282 236 
F6 at T34 Z22 284 
F7 at T12 1705 1513 
2209 2595 . 2033 155 
Condenser Input: -F/-atetl2 1705 iMag hs, 
F9 at T29 L738 1760 
3478 8275) 
Output F/7 at T1ll Uh Fag 109 
F9 at Tl 3080 2960 
3207 2 00) 204 


Interstage Loss F2(T19-T4) Lr U31 


Heat Transfer Coefficients BTU/hr*°F*£t2 


First Effect (FlatT15-FlatT5)*60 661 665 
: Le 
: ke 
Second Effect ane sethep 4.93 384 


Table A-3 
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Sample Heat Balance, Experiment FBI2 


Time Basis 1 Minute 
Heat Units BTU 


Section Considered Stream Scare 

Total Process Input Fl at T15 2250 

F8.at T7 792 

F9 at T29 L773 

Output llFl at 5 384 

FE5 at T2656 282 

Bo,at 134 Zee 

Pyeate LL i ag 

Borat Tl 3080 

4095 

Losses First Effect 78 

Interstage 113 

Second Effect 259 

Condenser Pag el 

Balance: Error % 0.0 
Heat Transfer Areas First Effect 9.88 £t2 


Second Effect 4.10 ft 
Liquid Hold-up at Steady State 


First Effect 23 lbs 
Second Effect 28 lbs 


Table A-3 continued 


End 


1 pe by 
794 


1760 
4532 


334 
236 
284 
109 


2960 
3923 


118 
131 
135 
204 
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DATA FOR EXPERIMENT OL I 2 


VARIABLE DESCRIPTION 


FLOW RATES IN LBS/MIN 


STEAM FLOW TO FIRST EFFECT 
FIRST EFFECT BOTTOMS FLOW 
FIRST EFFECT OVERHEAD FLOW 
PRODUCT FLOW 

SECOND EFFECT OVERHEAD FLOW 
TOTAL FEED FLOW 

COOLING WATER TO CONDENSER 


CONCENTRATIONS WEIGHT PER CENT 


FEED CONCENTRATION 
PRODUCT CONCENTRATION 


TEMPERATURES IN DEGREES F 


CONDENSER COOLING WATER OUT 
FIRST EFFECT VAPOUR 

SOLUTION TO SECOND EFFECT 
STEAM CONDENSATE FIRST EFFECT 
FEED TO FIRST EFFECT 

STEAM TO SECOND EFFECT 
CONDENSER CONDENSATE 
SEPARATOR VAPOUR 

STEAM SUPPLY TO EVAPORATOR 
LIQUID IN FIRST EFFECT 


STEAM CONDENSATE SECOND EFFECT 


CONDENSER COOLING WATER IN 
PRODUCT FROM SECOND EFFECT 


TABLE A= 4 


VARI 
ABLE 
NAME 


Cl 
C6 


Tl 


T4 

ao 

a7 

T10 
Tia 
T12 
T15 
119 
T28 
T29 
T34 


INITIAL 


STEADY 
STATE 


2001 
3014 
le7l 
1646 
170 
4096 
4303 


2040 
8200 


103200 
220200 
181240 
238430 
189240 
219280 
119290 
16050 
313400 
220004 
194430 

70200 
2072680 


FINAL 
STEADY 
STATE 


2201 
3014 
1669 
1642 
1654 
4099 
43612 


3000 
1008 


103430 
220210 
18280 
23870 
19040 
220200 
118220 
16240 
312480 
220070 
195200 

70600 
20830 
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DATA FOR EXPERIMENT FB I 2 


VARIABLE DESCRIPTION 


FLOW RATES IN LBS/MIN 


STEAM FLOW TO FIRST EFFECT 
FIRST EFFECT BOTTOMS FLOW 
FIRST EFFECT OVERHEAD FLOW 
PRODUCT FLOW 

SECOND EFFECT OVERHEAD FLOW 
TOTAL FEED. FLOW 

COOLING WATER TO CONDENSER 


CONCENTRATIONS WEIGHT PER CENT 


FEED CONCENTRATION 
PRODUCT CONCENTRATION 


TEMPERATURES IN DEGREES F 


CONDENSER COOLING WATER OUT 
FIRST EFFECT VAPOUR 

SOLUTION TO SECOND EFFECT 
STEAM CONDENSATE FIRST EFFECT 
FEED TO FIRST EFFECT 

STEAM TO SECOND EFFECT 
CONDENSER CONDENSATE 
SEPARATOR VAPOUR 

STEAM SUPPLY TO EVAPORATOR 
LIQUID IN FIRST EFFECT 


STEAM CONDENSATE SECOND EFFECT 


CONDENSER COOLING WATER IN 
PRODUCT FROM SECOND EFFECT 


TABLE A= 5 


VARI 
ABLE 
NAME 


Cl 
C6 


INITIAL 


STEADY 
STATE 


190 
3014 
1.73 
1¢51 
1¢51 
5002 
46067 


2079 
8497 


97250 
216450 
18050 
234020 
189460 
216440 
11570 
159220 
313470 
217430 
195630 

70200 
178670 
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STEADY 
STATE 
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11280 
155230 
314680 
216460 
19350 

70200 
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DATA FOR EXPERIMENT FF I 2 


VARIABLE DESCRIPTION 


FLOW RATES IN LBS/MIN 


STEAM FLOW TO FIRST EFFECT 
FIRST EFFECT BOTTOMS FLOW 
FIRST EFFECT OVERHEAD FLOW 
PRODUCT FLOW 

SECOND EFFECT OVERHEAD FLOW 
TOTAL FEED. FLOW 

COOLING WATER TO CONDENSER 


CONCENTRATIONS WEIGHT PER CENT 


FEED CONCENTRATION 
PRODUCT CONCENTRATION 


TEMPERATURES IN DEGREES F 


CONDENSER COOLING WATER OUT 
FIRST EFFECT VAPOUR 

SOLUTION TO SECOND EFFECT 
STEAM CONDENSATE FIRST EFFECT 
PEGR TO FIRST-EFFECT 

STEAM TO SECOND EFFECT 
CONDENSER CONDENSATE 
SEPARATOR VAPOUR 

STEAM SUPPLY TO EVAPORATOR 
LIQUID IN FIRST EFFECT 

STEAM CONDENSATE SECOND EFFECT 
CONDENSER COOLING WATER IN 
PRODUCT FROM SECOND EFFECT 


TABLE A= 6 


VARI 
ABLE 
NAME 


Fl 
F2 
FS 
F6 
F7 
F8 
F9 


Cl 
C6 


Tl 
T2 
T4 
tS 


T10 
Til 
T12 
T15 
T19 
T28 
T29 
T34 


INITIAL 


STEADY 
STATE 


1690 
3016 
164 
1663 
1251 
4098 
46052 


2060 
8239 


99400 
215420 
17350 
232280 
189440 
215210 
11850 
146280 
312230 
215470 
187220 
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19270 


A-11 
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STATE 


1065 
3056 
1*52 
1292 
1631 
4097 
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DATA FOR EXPERIMENT FFB I 


TABLE A= 7 


VARIABLE DESCRIPTION VARI INITIAL FINAL 
ABLE STEADY STEADY 
NAME STATE STATE 
FLOW RATES IN LBS/MIN 
STEAM FLOW TO FIRST EFFECT rs 1.87 le67 
FIRST EFFECT BOTTOMS FLOW a - 3027 3036 
FIRST EFFECT OVERHEAD FLOW ro 1659 1638 
PRODUCT FLOW F6 1654 1697 
SECOND EFFECT OVERHEAD FLOW ai 1650 1033 
TOTAL FEED. FLOW F8 4099 5000 
COOLING WATER TO CONDENSER F9 47209 46094 
CONCENTRATIONS WEIGHT PER CENT 
FEED CONCENTRATION Cl 2079 3039 
PRODUCT CONCENTRATION C6 8049 860 
TEMPERATURES IN DEGREES F 
CONDENSER COOLING WATER OUT Tl 9940 96010 
FIRST EFFECT VAPOUR T2 216410 215050 
SOLUTION TO SECOND EFFECT T4 17570 175480 
STEAM CONDENSATE FIRST EFFECT T5 233440 231210 
FEED TO FIRST EFFECT T7 19010 19090 
STEAM TO SECOND EFFECT T10 215290 215240 
CONDENSER CONDENSATE Thi 119290 116400 
SEPARATOR VAPOUR tie 15080 149470 
STEAM SUPPLY TO EVAPORATOR Fis 313610 314010 
LIQUID IN FIRST EFFECT T19 216260 216020 
STEAM CONDENSATE SECOND EFFECT T28 19090 18750 
CONDENSER COOLING WATER IN T29 70200 70200 
PRODUCT FROM SECOND EFFECT 134 186480 188640 
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DATA FOR EXPERIMENT FB I 3 


VARTABLE DESCRIPTION 


FLOW RATES IN LBS/MIN 


STEAM FLOW TO FIRST EFFECT 
FIRST EFFECT BOTTOMS FLOW 
FIRST EFFECT OVERHEAD FLOW 
PRODUCT FLOW 

SECOND EFFECT OVERHEAD FLOW 
TOTAL FEED FLOW 

COOLING WATER TO CONDENSER 


CONCENTRATIONS WEIGHT PER CENT 


FEED CONCENTRATION 
PRODUCT CONCENTRATION 


TEMPERATURES IN DEGREES F 


CONDENSER COOLING WATER OUT 
FIRST EFFECT VAPOUR 

SOLUTION TO SECOND EFFECT 
STEAM CONDENSATE FIRST EFFECT 
FEED TO FIRST EFFECT 

STEAM TO SECOND EFFECT 
CONDENSER CONDENSATE 
SEPARATOR VAPOUR 

STEAM SUPPLY TO EVAPORATOR 
LIQUID IN FIRST EFFECT 


STEAM CONDENSATE SECOND EFFECT 


CONDENSER COOLING WATER IN 
PRODUCT FROM SECOND EFFECT 


TABLE A= 8 


VARI 
ABLE 
NAME 


Cl 
C6 


Th 
Tz 
T4 
T5 


T10 
Ted 2 
T1l2 
T15 
th9 
T28 
tae 
T34 


INITIAL 


STEADY 
STATE 


1666 
3029 
1658 
1296 
1¢31 
5200 
47266 


3039 
8256 


96240 
215430 
175240 
231200 
191220 
215430 
11630 
149260 
3214200 
216420 
18270 

70200 
188480 


A=15 


FINAL 
STEADY 
STATE 


1090 
3018 
1061 
1053 
1250 
4099 
47.207 


2079 
8664 


9920 
214020 
173440 
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214020 
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147240 
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70200 
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DATA FOR EXPERIMENT FF I 3 


VARIABLE DESCRIPTION 


FLOW RATES IN LBS/MIN 


STEAM FLOW TO FIRST EFFECT 
FIRST EFFECT BOTTOMS FLOW 
FIRST EFFECT OVERHEAD FLOW 
PRODUCT FLOW 

SECOND EFFECT OVERHEAD FLOW 
TOTAL FEED. FLOW 

COOLING WATER TO CONDENSER 


CONCENTRATIONS WEIGHT PER CENT 


FEED CONCENTRATION 
PRODUCT CONCENTRATION 


TEMPERATURES IN DEGREES F 


CONDENSER COOLING WATER OUT 
FIRST EFFECT VAPOUR 

SOLUTION TO SECOND EFFECT 
STEAM CONDENSATE FIRST EFFECT 
PECD 70. FIRST EFFECT 

STEAM TO SECOND EFFECT 
CONDENSER CONDENSATE 
SEPARATOR VAPOUR 

STEAM SUPPLY TO EVAPORATOR 
LIQUID IN FIRST EFFECT 


STEAM CONDENSATE SECOND EFFECT 


CONDENSER COOLING WATER IN 
PRODUCT FROM SECOND EFFECT 


TABLE A= 9 


VARI 
ABLE 
NAME 


Cl 
C6 


Tl 
T2 
T4 
T$ 
Va 
T10 
Tid 
T12 
*i5 
9 
T28 
T29 
134 


INITIAL 


STEADY 
STATE 


1682 
3035 
155 
1.75 
1250 
4095 
4881 


32039 
9219 


106450 
229600 
19030 
240280 
196600 
228470 

97450 
172250 
3216240 
229250 
207260 

75200 
12620 


A=-17 


FINAL 
STEADY 
STATE 
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DATA FOR EXPERIMENT FFB I 2 


VARTABLE DESCRIPTION 


FLOW RATES IN LBS/MIN 


STEAM FLOW TO FIRST EFFECT 
FIRST EFFECT BOTTOMS FLOW 
FIRST EFFECT OVERHEAD FLOW 
PRODUCT FLOW 

SECOND EFFECT OVERHEAD FLOW 
TOTAL FEED. FLOW 

COOLING WATER TO CONDENSER 


CONCENTRATIONS WEIGHT PER CENT 


FEED CONCENTRATION 
PRODUCT CONCENTRATION 


TEMPERATURES IN DEGREES F 


CONDENSER COOLING WATER OUT 
FIRST EFFECT VAPOUR 

SOLUTION TO SECOND EFFECT 
STEAM CONDENSATE FIRST EFFECT 
PEED 1O FIRST EFFECT 

STEAM TO SECOND EFFECT 
CONDENSER CONDENSATE 
SEPARATOR VAPOUR 

STEAM SUPPLY TO EVAPORATOR 
LIQUID IN FIRST EFFECT 


STEAM CONDENSATE SECOND EFFECT 


CONDENSER COOLING WATER IN 
PRODUCT FROM SECOND EFFECT 


TABLE A=-10 


VARI INITIAL 
ABLE STEADY 


NAME STATE 
rs 2200 
ae: 3e1i 
ro 1672 
F6 140 
ae | 152 
F8 496 
F9 49.33 
ip 3030 
C6 10677 
14 10930 
ce 23120 
T4 191400 
TS 244640 
: FI 19620 
T10 231600 
Tae 10560 
T12 173460 
T23 315260 
a9 231280 
T28 2090670 
T29 75200 
T34 122600 
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FINAL 
STEADY 
STATE 


2e21l 
3608 
1686 
1e22 
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5202 
57064 


20/0 
1086 


110430 
231470 
19130 
246060 
195020 
23140 
133050 
174e20 
314480 
231290 
210680 

75200 
12060 
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Figure A- 8 
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DATA FOR EXPERIMENT INF I 


TABLE Aw11l 


VARIABLE DESCRIPTION VARI INITIAL FINAL 
ABLE STEADY STEADY 
NAME STATE STATE 
FLOW RATES IN LBS/MIN 
STEAM FLOW TO FIRST EFFECT ri 1.81 2201 
FIRST EFFECT BOTTOMS FLOW ra 3030 3009 
FIRST EFFECT OVERHEAD FLOW F5 1058 1279 
PRODUCT FLOW F6 le7l 1036 
SECOND EFFECT OVERHEAD FLOW ei 1.55 1065 
TOTAL FEED. FLOW F8 4098 4099 
COOLING WATER TO CONDENSER F9 4291 43017 
CONCENTRATIONS WEIGHT PER CENT 
FEED CONCENTRATION ae | 3209 2050 
PRODUCT CONCENTRATION C6 8030 831 
TEMPERATURES IN DEGREES F 
CONDENSER COOLING WATER OUT ie 9920 102650 
FIRST EFFECT VAPOUR Te 215290 219470 
SOLUTION TO SECOND EFFECT T4 179230 18080 
STEAM CONDENSATE FIRST EFFECT T5 232040 237680 
FEED TO FIRST EFFECT T7 189690 1894670 
STEAM TO SECOND EFFECT T10 215050 219430 
CONDENSER CONDENSATE Lee! 11550 120610 
SEPARATOR VAPOUR T12 157420 159480 
STEAM SUPPLY TO EVAPORATOR T15 31410 313400 
LIQUID IN FIRST EFFECT T19 216030 219480 
STEAM CONDENSATE SECOND EFFECT 728 188420 19330 
CONDENSER COOLING WATER IN T29 70200 70600 
PRODUCT FROM SECOND EFFECT T34 206010 207060 
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DATA FOR EXPERIMENT INF I 2 


TABLE Aw12 


VARIABLE DESCRIPTION VARI INITIAL FINAL 
ABLE STEADY STEADY 
NAME STATE STATE 
FLOW RATES IN LBS/MIN 
STEAM FLOW TO FIRST EFFECT Fi 1085 2002 
FIRST EFFECT BOTTOMS FLOW a 3e21 Zell 
FIRST EFFECT OVERHEAD FLOW FS 1057 1076 
PRODUCT FLOW F6 le7l 1460 
SECOND EFFECT OVERHEAD FLOW F7 153 1051 
TOTAL FEED. FLOW F8 4098 502 
COOLING WATER TO CONDENSER F9 49699 4880 
CONCENTRATIONS WEIGHT PER CENT 
FEED CONCENTRATION Cl 350 2089 
PRODUCT CONCENTRATION C6 9067 8097 
TEMPERATURES IN DEGREES F 
CONDENSER COOLING WATER OUT tea 106200 10950 
FIRST EFFECT VAPOUR T2 229240 231010 
SOLUTION TO SECOND EFFECT T4 18940 189480 
STEAM CONDENSATE FIRST EFFECT T5 241200 244010 
FEED TO FIRST EFFECT T7 19610 195070 
STEAM TO SECOND EFFECT T10 229200 2320090 
CONDENSER CONDENSATE ee 98 «30 113400 
SEPARATOR VAPOUR Fiz 171270 17230 
STEAM SUPPLY TO EVAPORATOR Tis 31620 315250 
LIQUID IN FIRST EFFECT T19 2290670 23160 
STEAM CONDENSATE SECOND EFFECT T28 207200 209010 
CONDENSER COOLING WATER IN , T29 75200 75400 
PRODUCT FROM SECOND EFFECT T34 12740 127050 
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DATA FOR EXPERIMENT OL II 2 


TABLE A=-13 


VARIABLE DESCRIPTION VARI INITIAL FINAL 
ABLE STEADY STEADY 
NAME STATE STATE 
FLOW RATES IN LBS/MIN 
STEAM FLOW TO FIRST EFFECT Fl 2201 201 
FIRST EFFECT BOTTOMS FLOW Pa 3013 3078 
FIRST EFFECT OVERHEAD FLOW F5 1270 1665 
PRODUCT FLOW F6 1042 2009 
SECOND EFFECT OVERHEAD FLOW ard 1251 1655 
TOTAL FEED FLOW F8 5200 5050 
COOLING WATER TO CONDENSER FS 42097 43602 
CONCENTRATIONS WEIGHT PER CENT 
FEED CONCENTRATION a 3000 3000 
PRODUCT CONCENTRATION C6 9058 7061 
TEMPERATURES IN DEGREES F 
CONDENSER COOLING WATER OUT Le 10340 103490 
FIRST EFFECT VAPOUR Ke 22040 2220690 
SOLUTION TO SECOND EFFECT T4 18320 185610 
STEAM CONDENSATE FIRST EFFECT T5 23880 241450 
FEED TO FIRST EFFECT ws 19070 191210 
STEAM TO SECOND EFFECT T10 220220 222070 
CONDENSER CONDENSATE Eis 117210 117260 
SEPARATOR VAPOUR T12 163420 164620 
STEAM SUPPLY TO EVAPORATOR T15 313200 312460 
LIQUID IN FIRST EFFECT T19 220690 223040 
STEAM CONDENSATE SECOND EFFECT T28 19580 197010 
CONDENSER COOLING WATER IN T29 70200 70000 
PRODUCT FROM SECOND EFFECT T34 208460 208090 
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DATA FOR EXPERIMENT FB II 2 


TABLE A=-14 


VARIABLE DESCRIPTION VARI INITIAL FINAL 
ABLE STEADY STEADY 
NAME STATE STATE 
FLOW RATES IN LBS/MIN 
STEAM FLOW TO FIRST EFFECT Fl 1681 2020 
FIRST EFFECT BOTTOMS FLOW F2 3002 3041 
FIRST EFFECT OVERHEAD FLOW FS 1235 166 
PRODUCT FLOW F6 1638 1679 
SECOND EFFECT OVERHEAD FLOW aot J 1250 1067 
TOTAL FEED. FLOW F8 4050 5048 
COOLING WATER TO CONDENSER F9 48035 48.78 
CONCENTRATIONS WEIGHT PER CENT 
FEED CONCENTRATION Cl 2280 280 
PRODUCT CONCENTRATION C6 9204 9ell 
TEMPERATURES IN DEGREES F 
CONDENSER COOLING WATER OUT :e! 96620 101260 
FIRST EFFECT VAPOUR T2 216660 221670 
SOLUTION TO SECOND EFFECT T4 17840 18140 
STEAM CONDENSATE FIRST EFFECT T& 23360 241040 
FEED TO FIRST EFFECT T7 189420 189.290 
STEAM TO SECOND EFFECT T10 216640 221240 
CONDENSER CONDENSATE Til 114250 121-80 
SEPARATOR VAPOUR pp. 154200 157240 
STEAM SUPPLY TO EVAPORATOR T15 318250 314290 
LIQUID IN FIRST EFFECT T19 217¢30 222050 
STEAM CONDENSATE SECOND EFFECT T28 193420 197280 
CONDENSER COOLING WATER IN T29 65200 65200 
PRODUCT FROM SECOND EFFECT T34 195670 196200 
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DATA FOR EXPERIMENT FF II 2 


VARIABLE DESCRIPTION 


FLOW RATES IN LBS/MIN 


STEAM FLOW TO FIRST EFFECT 
FIRST EFFECT BOTTOMS FLOw 
FIRST EFFECT OVERHEAD FLOW 
PRODUCT FLOW 

SECOND EFFECT OVERHEAD FLOW 
TOTAL FEED FLOW 

COOLING WATER TO CONDENSER 


CONCENTRATIONS WEIGHT PER CENT 


FEED CONCENTRATION 
PRODUCT CONCENTRATION 


TEMPERATURES IN DEGREES F 


CONDENSER COOLING WATER OUT 
FIRST EFFECT VAPOUR 

SOLUTION TO SECOND EFFECT 
STEAM CONDENSATE FIRST EFFECT 
FEED 10 FIRST EFFECT 

STEAM TO SECOND EFFECT 
CONDENSER CONDENSATE 
SEPARATOR VAPOUR 

STEAM SUPPLY TO EVAPORATOR 
LIQUID IN FIRST EFFECT 


STEAM CONDENSATE SECOND EFFECT 


CONDENSER COOLING WATER IN 
PRODUCT FROM SECOND EFFECT 


TABLE Aw15 


VARI INITIAL 
ABLE STEADY 


NAME STATE 
a | 1279 
F2 2094 
FS 1630 
F6 1640 
rT 1650 
F8 4649 
F9 49 «22 
ay | 2069 
C6 9049 
Tl 98240 
T2 220290 
T4 182400 
Lie 237250 
TT? 18920 
T10 220650 
Lae! 109460 
1 he 159480 
TiS 319010 
T19 221260 
T28 197200 
T29 65200 
T34 1984620 


FINAL 
STEADY 
STATE 
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315470 
222440 
198420 
6500 
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DATA FOR EXPERIMENT FF IV 2 


VARIABLE DESCRIPTION 


FLOW RATES IN LBS/MIN 


STEAM FLOW TO FIRST EFFECT 
FIRST EFFECT BOTTOMS FLOW 
FIRST EFFECT OVERHEAD FLOW 
PRODUCT FLOW 

SECOND EFFECT OVERHEAD FLOW 
TOTAL FEED. FLOW 

COOLING WATER TO CONDENSER 


CONCENTRATIONS WEIGHT PER CENT 
FEED CONCENTRATION 
PRODUCT CONCENTRATION 
TEMPERATURES IN DEGREES F 
CONDENSER COOLING WATER OUT 


FIRST EFFECT VAPOUR 
SOLUTION TO SECOND EFFECT 


STEAM CONDENSATE FIRST EFFECT 


FEED TO FIRST EFFECT 

STEAM TO SECOND EFFECT 
CONDENSER CONDENSATE 
SEPARATOR VAPOUR 

STEAM SUPPLY TO EVAPORATOR 
LIQUID IN FIRST EFFECT 


STEAM CONDENSATE SECOND EFFECT 


CONDENSER COOLING WATER IN 
PRODUCT FROM SECOND EFFECT 


TABLE Aw16 


VARI 
ABLE 
NAME 


Ci 
C6 


INITIAL 


STEADY 
STATE 


184 
2081 
le72 
120 
1230 
4050 
49263 
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228290 
189260 
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195210 
228060 
9430 
17180 
317260 
229 240 
207210 
75400 
90200 


A=-31 


FINAL 
STEADY 
STATE 
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DATA FOR EXPERIMENT FFB II 3 


TABLE A~w17 


VARIABLE DESCRIPTION VARI INITIAL FINAL 
ABLE STEADY STEADY 
NAME STATE STATE 
FLOW RATES IN LBS/MIN 
STEAM FLOW TO FIRST EFFECT Fi 1.92 2230 
FIRST EFFECT BOTTOMS FLOW F2 2072 3061 
FIRST EFFECT OVERHEAD FLOW ‘ak: 163 1089 
PRODUCT FLOW F6 12616 1069 
SECOND EFFECT OVERHEAD FLOW F7 1655 1681 
TOTAL FEED FLOW F8 4e51 5051 
COOLING WATER TO CONDENSER F9 39051 50001 
CONCENTRATIONS WEIGHT PER CENT 
FEED CONCENTRATION Cl 3200 3000 
PRODUCT CONCENTRATION C6 1126 1118 
TEMPERATURES IN DEGREES F 
CONDENSER COOLING WATER OUT Lp 10540 109670 
FIRST EFFECT VAPOUR re 232420 231690 
SOLUTION TO SECOND EFFECT T4 19040 191610 
STEAM CONDENSATE FIRST EFFECT +S 245400 248620 
FEED TO FIRST EFFECT be 195200 195490 
STEAM TO SECOND EFFECT T10 231670 23160 
CONDENSER CONDENSATE T1l 98040 138480 
SEPARATOR VAPOUR T12 17320 173200 
STEAM SUPPLY TO EVAPORATOR TiS 316060 313620 
LIQUID IN FIRST EFFECT Ti? 232050 232290 
STEAM CONDENSATE SECOND EFFECT T28 210220 210290 
CONDENSER COOLING WATER IN T29 75000 75000 
PRODUCT FROM SECOND EFFECT 134 97210 95020 
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DATA FOR EXPERIMENT INF II 


VARIABLE DESCRIPTION 


FLOW RATES IN LBS/MIN 


STEAM FLOW TO FIRST EFFECT 
FIRST EFFECT BOTTOMS FLOW 
FIRST EFFECT OVERHEAD FLOW 
PRODUCT FLOW 

SECOND EFFECT OVERHEAD FLOW 
TOTAL FEED FLOW 

COOLING WATER TO CONDENSER 


CONCENTRATIONS WEIGHT PER CENT 


FEED CONCENTRATION 
PRODUCT CONCENTRATION 


TEMPERATURES IN DEGREES F 


CONDENSER COOLING WATER OUT 
FIRST EFFECT VAPOUR 

SOLUTION TO SECOND EFFECT 
STEAM CONDENSATE FIRST EFFECT 
FEED TO FIRST EFFECT 

STEAM TO SECOND EFFECT 
CONDENSER CONDENSATE 
SEPARATOR VAPOUR 

STEAM SUPPLY TO EVAPORATOR 
LIQUID IN FIRST EFFECT 


STEAM CONDENSATE SECOND EFFECT 


CONDENSER COOLING WATER IN 
PRODUCT FROM SECOND EFFECT 


TABLE A=-18 


VARI 
ABLE 
NAME 


Ci 
C6 


INITIAL 


STEADY 
STATE 


1085 
2075 
1662 
le21 
146 
4049 
47036 


300 
10298 


10160 
217290 
17740 
235090 
18940 
218200 
121460 
15280 
317200 
219200 
193620 

70200 
154270 


A=35 


FINAL 
STEADY 
STATE 
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1038 
10699 
5048 
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3200 
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DATA FOR EXPERIMENT INF IT 2 


TABLE A=19 


VARIABLE DESCRIPTION VARI INITIAL FINAL 
ABLE STEADY STEADY 
NAME STATE STATE 
FLOW RATES IN LBS/MIN 
STEAM FLOW TO FIRST EFFECT ri 190 2031 
FIRST EFFECT BOTTOMS FLOW FZ 2281 3044 
FIRST EFFECT OVERHEAD FLOW FS 1663 1.88 
PRODUCT FLOW F6 1.19 1250 
SECOND EFFECT OVERHEAD FLOW aad Leas 1290 
TOTAL FEED FLOW F8 4051 5048 
COOLING WATER TO CONDENSER F9 50603 45066 
CONCENTRATIONS WEIGHT PER CENT 
FEED CONCENTRATION | 3200 3600 
PRODUCT CONCENTRATION C6 10077 9296 
TEMPERATURES IN DEGREES F 
CONDENSER COOLING WATER OUT Le 106270 115670 
FIRST EFFECT VAPOUR 2 231260 23210 
SOLUTION TO SECOND EFFECT T4 191460 19160 
STEAM CONDENSATE FIRST EFFECT T5 245210 248460 
PEED TO FIRST EFFECT re 195420 19580 
STEAM TO SECOND EFFECT T10 231230 231290 
CONDENSER CONDENSATE ee | 97230 143020 
SEPARATOR VAPOUR Tid 174470 174620 
STEAM SUPPLY TO EVAPORATOR TiS 317+30 313630 
LIQUID IN FIRST EFFECT rae 23180 2320660 
STEAM CONDENSATE SECOND EFFECT T28 210210 211610 
CONDENSER COOLING WATER IN ee) 75200 75200 
PRODUCT FROM SECOND EFFECT T34 115200 107290 
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DATA FOR EXPERIMENT FB II 3 


TABLE A=20 


VARTABLE DESCRIPTION VARI INITIAL FINAL 
ABLE STEADY STEADY 
NAME STATE STATE 
FLOW RATES IN LBS/MIN 
STEAM FLOW TO FIRST EFFECT Fl 2019 1685 
FIRST EFFECT BOTTOMS FLOW F2 370 2090 
FIRST EFFECT OVERHEAD FLOW FS 1058 1034 
PRODUCT FLOW F6 1.77 1630 
SECOND EFFECT OVERHEAD FLOW F7 1.78 1055 
TOTAL FEED FLOw F8 5049 4e51 
COOLING WATER TO CONDENSER e9 49240 49024 
CONCENTRATIONS WEIGHT PER CENT 
FEED CONCENTRATION Ci 2080 280 
PRODUCT CONCENTRATION Ccé 9204 9206 
TEMPERATURES IN DEGREES F 
CONDENSER COOLING WATER OUT Tl 101670 9680 
FIRST EFFECT VAPOUR T2 221260 219400 
SOLUTION TO SECOND EFFECT T4 18140 18120 
STEAM CONDENSATE FIRST EFFECT T5 241240 235490 
FEEO fO FIRST EFFECT Tr 19030 1894620 
STEAM TO SECOND EFFECT T10 221430 218050 
CONDENSER CONDENSATE Tll 121480 11160 
SEPARATOR VAPOUR py: 15780 15930 
STEAM SUPPLY TO EVAPORATOR 715 314490 318430 
LIQUID IN FIRST EFFECT Ti9 222050 219450 
STEAM CONDENSATE SECOND EFFECT T28 19820 195200 
CONDENSER COOLING WATER IN 29 65200 65200 
PRODUCT FROM SECOND EFFECT T34 19590 196060 
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DATA FOR EXPERIMENT FF II 3 


VARIABLE DESCRIPTION 


FLOW RATES IN LBS/MIN 


STEAM FLOW TO FIRST EFFECT 
FIRST EFFECT BOTTOMS FLOW 
FIRST EFFECT OVERHEAD FLOW 
PRODUCT FLOW 

SECOND EFFECT OVERHEAD FLOW 
TOTAL FEED FLOW 

COOLING WATER TO CONDENSER 


CONCENTRATIONS WEIGHT PER CENT 


FEED CONCENTRATION 
PRODUCT CONCENTRATION 


TEMPERATURES IN DEGREES F 


CONDENSER COOLING WATER OUT 
FIRST EFFECT VAPOUR 

SOLUTION TO SECOND EFFECT 
STEAM CONDENSATE FIRST EFFECT 
rEeO 10: FikSt EFFECT 

STEAM TO SECOND EFFECT 
CONDENSER CONDENSATE 
SEPARATOR VAPOUR 

STEAM SUPPLY TO EVAPORATOR 
LIQUID IN FIRST EFFECT 

STEAM CONDENSATE SECOND EFFECT 
CONDENSER COOLING WATER IN 
PRODUCT FROM SECOND EFFECT 


TABLE Aw21 


VARI INITIAL 
ABLE STEADY 


NAME STATE 
Pi 2e31 
r2 3043 
FS 1685 
F6é 1053 
F7 le7l 
F8 5250 
F9 43038 
Cl 3200 
C6 1028 
T1 111.200 
T2 232210 
T4 191230 
‘. 248040 
TT 19630 
T10 231270 
ge 1328260 
he 17380 
co 31320 
T19 232050 
T28 211470 
T29 75200 
T34 10260 
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1073 


105280 
23170 
190640 
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231440 
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316270 
232400 
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98010 
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DATA FOR EXPERIMENT FF IV 1 


TABLE A=22 


VARIABLE DESCRIPTION VARI INITIAL FINAL 
ABLE STEADY STEADY 
NAME STATE STATE 
FLOW RATES IN LBS/MIN 
STEAM FLOW TO FIRST EFFECT rt 2031 1680 
FIRST EFFECT BOTTOMS FLOW re 3052 2079 
FIRST EFFECT OVERHEAD FLOW we, 2001 1.57 
PRODUCT FLOW F6 152 1019 
SECOND EFFECT OVERHEAD FLOW Fr 1263 1037 
TOTAL FEED FLOW F8 505] 4048 
COOLING WATER TO CONDENSER r9 4948 49036 
CONCENTRATIONS WEIGHT PER CENT 
FEED CONCENTRATION Cl 3200 300 
PRODUCT CONCENTRATION C6 1038 1038 
TEMPERATURES IN DEGREES F 
CONDENSER COOLING WATER OUT T1 10940 103610 
FIRST EFFECT VAPOUR T2 231¢70 2284200 
SOLUTION TO SECOND EFFECT T4 19060 189480 
STEAM CONDENSATE FIRST EFFECT T5 247090 240010 
FEEO 'O FIRST EFFECT T7 19640 194630 
STEAM TO SECOND EFFECT T10 231430 227040 
CONDENSER CONDENSATE yan 138490 93090 
SEPARATOR VAPOUR tie &- 172270 172610 
STEAM SUPPLY TO EVAPORATOR Tis 31290 aL feiO 
LIQUID IN FIRST EFFECT Tis 232060 228430 
STEAM CONDENSATE SECOND EFFECT T28 210090 206090 
CONDENSER COOLING WATER IN tae 75200 75000 
PRODUCT FROM SECOND EFFECT T34 94030 9140 
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A=45 


DATA FOR EXPERIMENT FFB II 2 


TABLE A=23 


VARIABLE DESCRIPTION VARI INITIAL FINAL 
ABLE STEADY STEADY 
NAME STATE STATE 
FLOW RATES IN LBS/MIN 
STEAM FLOW TO FIRST EFFECT ae 2019 1.80 
EIRST EFFECT BOTTOMS FLOW ae 3051 2090 
FIRST EFFECT OVERHEAD FLOW FS 1265 133 
PRODUCT FLOW F6 le7l 1040 
SECOND EFFECT OVERHEAD FLOW rt 1280 1250 
TOTAL FEED FLOW F8 5050 4050 
COOLING WATER TO CONDENSER F9 4878 49618 
CONCENTRATIONS WEIGHT PER CENT 
FEED CONCENTRATION Cl 2080 280 
PRODUCT CONCENTRATION C6 8059 8053 
TEMPERATURES IN DEGREES F 
CONDENSER COOLING WATER OUT Tk 102680 9730 
FIRST EFFECT VAPOUR T2 221670 219430 
SOLUTION TO SECOND EFFECT T4 18280 180290 
STEAM CONDENSATE FIRST EFFECT 1 240290 235000 
FEED TO FIRST EFFECT Vo 189660 188430 
STEAM TO SECOND EFFECT T10 2214250 219610 
CONDENSER CONDENSATE Til 12280 113020 
SEPARATOR VAPOUR T12 159430 158280 
STEAM SUPPLY TO EVAPORATOR TPS 315480 319200 
LIQUID IN FIRST EFFECT T19 222670 219480 
STEAM CONDENSATE SECOND EFFECT T28 19820 196200 
CONDENSER COOLING WATER IN T29 65200 65400 
PRODUCT FROM SECOND EFFECT T34 19780 197250 
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DATA FOR EXPERIMENT OL IV 


TABLE A=24 


VARIABLE DESCRIPTION VARI INITIAL FINAL 
ABLE STEADY STEADY 
NAME STATE STATE 
FLOW RATES IN LBS/MIN 
STEAM FLOW TO FIRST EFFECT Pi 1658 1.80 
FIRST EFFECT BOTTOMS FLOW F2 3001 2078 
FIRST EFFECT OVERHEAD FLOW ro 130 1052 
PRODUCT FLOW F6 1662 le27 
SECOND EFFECT OVERHEAD FLOW as le27 1640 
TOTAL FEED FLOW F8 4053 4048 
COOLING WATER TO CONDENSER re? 4949 48079 
CONCENTRATIONS WEIGHT PER CENT 
FEED CONCENTRATION a | 3018 3007 
PRODUCT CONCENTRATION C6 9256 1129 
TEMPERATURES IN DEGREES F 
CONDENSER COOLING WATER OUT te 9520 9740 
FIRST EFFECT VAPOUR T2 210210 219080 
SOLUTION TO SECOND EFFECT T4 17720 17720 
STEAM CONDENSATE FIRST EFFECT T5 226210 237090 
Feeo TO FIRST EFFECT te 191.200 19130 
STEAM TO SECOND EFFECT T10 209290 219650 
CONDENSER CONDENSATE te 110210 117630 
SEPARATOR VAPOUR Tie 15280 152290 
STEAM SUPPLY TO EVAPORATOR 125 318060 317610 
LIQUID IN FIRST EFFECT Loe 211210 220690 
STEAM CONDENSATE SECOND EFFECT T28 194020 19230 
CONDENSER COOLING WATER IN T29 65200 65400 
PRODUCT FROM SECOND EFFECT T34 192270 194020 
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VARIABLE DESCRIPTION 


FLOW RATES IN LBS/MIN 


STEAM FLOW TO FIRST EFFECT 
FIRST EFFECT BOTTOMS FLOW 
FIRST EFFECT OVERHEAD FLOW 
PRODUCT FLOW 

SECOND EFFECT OVERHEAD FLOW 
TOTAL FEED FLOW 

COOLING WATER TO CONDENSER 


CONCENTRATIONS WEIGHT PER CENT 
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PRODUCT CONCENTRATION 


TEMPERATURES IN DEGREES F 
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The following six pages contain graphs showing the 
transient behaviour of the steam flow for all expe- 
riments in which a positive step in feed flow rate 
was applied. They have been included here to illus- 
trate the controller output obtained for different 


types of controllers. 
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Appendix B 


Determination of Trans- 
portation Lags 


List of Orifice and 
Valve Stem Sizes 


Calibration Curves 
Process Flow Sheet 


Legend to Process Flow 
Sheet 
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Determination of Transportation Lags 


Feed Piping 
pipe diameter 3/8 
pipe length 27 ft 
pipe cross section 


equivalent lengths 


sl koh pays 


0.11 in 


of resistances (25) 


2 globe valves 16 ft each 
7 elbows 3 ft each 
orifice and refractometer 0.6 ft each 


total equivalent pipe length 81 ft 


flow rate 5.0 lbs/min specific gravity 1.0 


velocity 105 ft/min 


transportation lag 


Interstage Piping 
pipe diameter 3/8 


pipe length 11 ft 
pipe cross section 
equivalent lengths 


3 globe valves 
izoritice =0.6 


81/105 = 0.78 min 


ineivd. 


0.11 in2 
of resistances 


16 ft each 
ft 


total equivalent pipe length 60 ft 


flow rate 3.30 1bs 


velocity 69 ft/min 


transportation lag 


/min specific gravity 1.0 


60/69 = 0.86 min 


Table B- 1 
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List of Orifice and Valve Stem Sizes 


Variable Concerned Orifice Diameter Valve Constant Oe 


Inches 
Fl 0.382 2,25 
F2 0.125 1.50 
F5 0.100 0.56 
F6 0.100 0.56 
F7 0.150 
F8 0.160 
F9 0.327 nos 
F10 0.875 
F1l 0.125 0.20 
F12 0.075 0.125 
ri 1.13 
L13 0.56 
DVP21 0.012 
P22 0.012 
TT10 0.28 
TT11 0.28 


Table B- 2 
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Figure B-19 
Process Flow Sheet of Double Effect Evaporator @) 
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Legend to Figure B-19 


B-23 


Note: The DDC facilities allow for a wide range of con- 
trol schemes to be implemented readily. The con- 
troller configuration employed in this study is 
shown on Figures 3.1 and 3.2. 


Vessels: 


E-1 
E-2 
E-3 
E-5 
E-10 
E-11 
E-13 
E-14 
E-15 
E-16 
E-17 
E-19 


first effect main body 
storage tank 

service side condensate tank 
storage tank 

storage tank 

water tank 

separator 

second effect main body 
condenser 

process side condensate tank 
rundown tank 


damper 


heat exchanger 


cartridge filter 


vacuum manifold 


Table B-3 
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stream identification 


hand valve 


solenoid valve 


check valve 


pneumatic control valve 


pneumatic control valve with valve 
positioner 
orifice 


cold water header 


steam header 


vacuum header 
pressure gauge 


sight glass 


thermocouple 


flow transmitter 


level transmitter 


Table B-3 continued 


ey: pressure transmitter 


Ss concentration transmitter 


eS = ‘S) a horizontal bar indicates the 
existence of a panel mounted re- 


CE © corder for the variable concerned 


O C)} O an outer circle indicates that the 


variable can be read by the computer 


Valve Labels solenoid valves are identified as 
per reference (32) 


pneumatic control valves are identi- 


fied by the letter V followed by the 
variable regulated as per Table B-2 


Table B-3 continued 
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Appendix C 


Control and Filter Algorithms 
Calculations 
Feedforward Controller 
Inferential Controller 


Building Blocks on 
Figure 5.4 


Analytical Solution of 
Individual Equations 


C- 2 to C= 3 
C- 4 
C- 7 
C- 8 
C- 8 to C-10 
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Control and Pilter Algorithms 


Digital Controllers 


proportional = KP*E + B 
integral = KP*KI*E + REST + B 
proportional -integral V = KP*E + KP*KI* 
PeRi od ct B 
Digital jFilters 
exponential filter PVO = (1-FF)*PVI + FF*PVOP 


In(FF)*TI = -DT T 


Union filter PVO = PVI 
when | PVI-PVOP | =p 


PVO = PVOP + (sign of PVI)*D 
when | PVI-PVoP| > D 
D=M 
when | SP-PvoP| <M 
D = |SP-PvoP| 


when | SP-PVoP| >M 


+ recommended for use 0.8=FF=1.0 (see page D-11) 


Table C-l 
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DT 


FF 
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pata 
PVvO 
PVOP 


REST 


SP 
TI 


Definition of Symbols 


output bias 
poll time interval 


effective error = filtered input - setpoint 
+ error bias 


filter factor, positive constant, O<FF=1 
DT/TI 
proportional constant 


quasi filter factor, constant defining 
allowable distance | PVO-PVOP| 


raw value of process variable 
filtered value of process variable 


filtered value of process variable at 
previous poll 


sum of previous integral contributions 
to V 


setpoint 


integral (reset) time as used in conti- 
nuous controllers 


control loop output 


Table C-1l continued 
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Calculations 


Symbols as defined on Figure 5.1 


Feedforward Controller 


Derivation of Equation (5-9) 


Given: process model in form of equations (5-7) and (5-8) 
Required: relation WST/CI 
Solution: 

control criterion CO = constant 


dcO 
dt = 0 (C-1) 


consider changes in CI independently, assume 
WI = constant (C-2) 
substitute (C-1) and (C-2) into (5-8) 
O = WI*C2 - 0.5*STECKWST*C2 - WI*CO 
+ STEC*WST*CO (C-4) 


linearize (C-3) by using perturbation techniques 


co xCl aa ke 
iotee we JS LEA = ae AGORA Sern s (c-4) 


substitute (C-2) into (5-7) 


ace 2 at WI*CI - aTeWI*C2 + ST WST*C2 (C-5) 
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linearize (C-5) by using perturbation techniques 


91 es ae SS. Kaede 
H1xgS2_ = (WI)*CI' + (0.5*STECAWST - WI)*C2' 


+ (0.5*STECKC2)*WST' (C-6) 


substitute (C-4) into (C-6) to obtain (5-9) 


Derivation of Equation (5-10) 


Given: equation (5-9) and steady state values from 
Table A-1 


Required: transfer function WST(s) /CI(s) 
Solution; 
STEC =(O1°+.02) /WST = 1.67 
C2 = WI*CI/B1 = 4.85% 
substitute steady state values in equation (5-9) 


1 
qiST = -0.060*CI' - 0.189*WST' (C-7) 





apply Laplace transformation to obtain (5-10) 





Derivation of Equation (5-. 


Given: equation (5-10) and steady state values from 
Table A-1 


Required: transfer function WST(s) /WI(s) 
Solution: 


assume first order lag system with time constant 
equal to that of equation (5-10) 
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to determine gain factor, assume constant steam eco- 
nomy which is realistic according to table 7.2 


STEC = (O01 + 02)/WST = 1.67 

fraction of feed evaporated at initial state 
(Ol + 02)/WI = 0.67 

fraction of feed evaporated at disturbed state 


delta(Ol + 02)/delta(WI) 


0.67 
STEC*delta(WST) /delta (WI) 


it follows that 
delta(WST) = 0.67*delta(WI)/STEC = 0.40*delta(WI) 


where delta designates the deviation from steady state 


Derivation of Equation (5-12) 


Given: steady state values from table A-1l 
Required: transfer function WST(s)/B1(s) 
Solution: 


analogous to the previous derivation 
assume constant steam economy 


ratio of quantity evaporated to first effect bottoms 
flow at initial state 


(Ol + 02)/Bl = 1,01 
at disturbed state 


delta(0Ol + 02)/delta(B1) EROL 


STEC*delta(WST) /delta(B1) 
it follows that 
delta(WST) = 1.01*delta(Bl)/STEC = 0.61*delta(B1) 
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Inferential Controller 


Derivation of Equation (5-13) 
Given: equation (5-7) 
Required: equation (5-13) 
Solution: 


linearize (5-7) by using perturbation techniques 


dC2 ! 


qo = CSTIC*CL' + CSTOWKWI' + CST2A*C2' 
+ CST2W*WST' (C-8) 
where 
CST1C = WI/H1 
CST2A = CST1B = CST9C = (STEC*WST - 2*WI) /(2*H1) 


CSTOW = (CI - C2)/H1 
CST2W 


STEC*C2 /(2*H1) 


hie oso aL to (C-8) to obtain 
5-13 


Derivation of Equation (5-14) 
Given: equation (5-8) 
Required: equation (5-14) 
Solution: 


linearize (5-8) by using perturbation techniques 
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! 
gee = CST4C*C2' + CSTTW*WL' + CST5D*CO' 
+ CST5W*WST' (C-9) 
where 


CST4C = (2*WI - STEC*WST) /(2*H2) 

CSTTW = (C2 - CO) /H2 

CST5D = CST4E = CSTIF = (STEC*WST - WI) /H2 
CST5W = (2*STECKCO - STEC*C2) /(2*H2) 


(oRtyy a cool transformation to (C-9) to obtain 


Building Blocks on Figure 5.4 


The building blocks on Figure 5.4 represent the 
following expressions 


block 1 CST1C/(s - CST1B) 
block 2 CST2W/(s - CST2A) 
block 9 CSTOW/(s - CST9C) 
block 4 -CST4C/(s - CST4E) 
block 5 CST5W/(s - CST5D) 
block 10 CSTTW/(s - CSTTF) 


Analytical Solution of Individual Equations 


In the solution of equations (5-13) and (5-14), each 
term was treated separately and the results added in 
accordance with the block diagram shown in Figure 5.4. 





Example: block 1 


The equivalent differential equation of the first 
term in equation (5-13) is 


a CST1B*B + CSTI1C*CI (C-10) 


The solution was obtained via Laplace transformation, 
considering initial conditions, B (0). 


s*B(s) - B(o) = CST1B*B(s) + csTiceSt 


(C-11) 
ARCO FL CSTICKCT : 
B(s) Or s(s - CSTIB) (C-12) 
B(t) = B(o)*exp(CST1B*t) 
, CSTIGHER 


<CSTIB * ( 1 - exp(CSTLB*t)) (C-13) 


The value of the constants was calculated using 
steady state readings at the beginning of each run. 
The value of the initial conditions B(o) was then 
determined from the steady state value of CI obtained 
from the process variable table and the relation 
required by equation (C-10) as follows 


B(o) = - Seerte (6-14) 


The value of t was constant and equal to the compu- 
tation interval H. The value of CI was obtained di- 
rectly from the process variable table every time a 
computation was called for. The values of B(t) at 
each sampling interval during the transient were ewva- 
luated from a difference equation based on (C-13) of 
the following general form 


B[(nt1)H] = B[nk] * K1 + CI[nH] * K2 (C-15) > 
where 


Kl = exp(CST1B*H) 
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Recess oat * ( 1 - exp(CST1B*H)) 


Calculations were performed by subroutine BLUE of 
program COLOR. 
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Appendix D 


List of Variable Records D- 2 
Process Variable Table D- 3 to D- 6 
Controller Configurations D- 7 to D-12 
Controller Constants D-13 
Controller Tuning D-14 to D-23 


Polling Sequence D-24 
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| D= 2 
COMPLETE LIST OF VARIABLE RECORDS 


COS = CURRENT OUTPUT STATION 
MPX = RELAY MULTIPLEXER POINT 


LOOP NAME COS MPX VARIABLE DESCRIPTION 

Oigerily5 28 seGOnuw EFFECT CHEST LEVEL 

Ob2e% Fr 7 42 SECOND EFFECT OVERHEAD FLOW 

0127 DATA ACCUMe FROM LOOP 0152 

0128 DATA ACCUMe FROM LOOP 0151 

0129 23 ARTIFICIAL READING OF PRODUCT CONCe 
0130 Je ARTIFICIAL READING OF FEED CONC. 
0131 DATA ACCUMe FROM LOOP 0132 

0332" Fi Tii=6 0 STEAM FLOW TO FIRST EFFECT 

G4a35 <F9 T1i6 16 COOLING WATER FLOW TO CONDENSER 


0134 P22 Tit=3 a3 SECOND EPFECT PRESSURE 
G233 7"P 20 T1I=4 4] FIRST EFFECT PRESSURE 


0136 L14 22 FIRST EFFECT LIQUID LEVEL 

GL37- 1.11 19 SEPARATOR LIQUID LEVEL 

0138 TT1O0 II=5 84 TEG FEED SOLUTION TEMPERATURE 
0139 F12 [1-4 9 TEG FEED SOLUTION FLOW 

0140 ; RATIO OF FEED FLOWS 

0141 Fill lI<3 S FEED WATER FLOW 

0142 TT11 Ire 83 FEED WATER TEMPERATURE 

0143 ¢1 32 FEED CONCENTRATION 

0144 C6 35 PRODUCT CONCENTRATION 

0145 32 FEEDFORWARD COMPENSATION FOR CONC. 
0146 DATA ACCUMe FROM LOOP 0137 

0147 DATA ACCUMe FROM LOOP 0144 

0148 . DATA ACCUMe FROM LOOP 0136 

0149 F8 12 TOTAL FEED FLOW 

0150 F5 38 FIRST EFFECT OVERHEAD FLOW 

0151 F2 IIlI=2 39 FIRST EFFECT BOTTOMS FLOW 

0152 F6 IlI=1 40 PRODUCT FLOW FROM SECOND EFFECT 
0153 RBT2 65 BRIDGE UNBALANCE MEASUREMENT 

0154 RBTl 64 REFERENCE VOLTAGE 

0155 Tl 100 COOLING WATER AT CONDENSER OUTLET 
0156 T2 90 FIRST EFFECT VAPOUR SPACE 

0157 T4 92 SOLUTION TO SECOND EFFECT 

0158 T5 88 STEAM CONDENSATE FROM FIRST EFFECT 
0159 T7 91 FEED TO FIRST EFFECT 

0160 710 97 STEAM TO SECOND EFFECT 

G16) 9 ‘¥H1 103 CONDENSER CONDENSATE 

O162°°T12 96 SEPARATOR VAPOUR 

0163 T15 86 STEAM SUPPLY TO EVAPORATOR 

0164 719 89 LIQUID IN FIRST EFFECT 

0165 T1728 98 SECOND EFFECT STEAM CONDENSATE 
0166 T30 101 RUNDOWN TANK COOLING WATER OUT 
0167 T34 79 PRODUCT FROM SECOND EFFECT 

0168 12 FEEDFORWARD COMPENSATION FOR FLOW 
0169 DUMMY LOOP FOR ADDITION OF SIGNALS 
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EXAMPLE OF PROCESS VARIABLE TABLE 


COLUMN 1 - VARIABLE RECORD IDENTIFICATION 
COLUMNS 2 TO 10 = VARIABLE RECORD CONTENTS 


012501 E21E 374D+00028 3C00 9540+03920-00840+15944415731 
012502+32767+32767 32204+32767-32768 2205+32767=32768 3220 
012503 0150419648 COAl 0000+18720-17920 0001 0200 0000 
012504 0000+00000 

012601 E0O0D 4F5D+00042 3200 84404+07500+00000+14754+00000 
012602+32767+32767 3220 

012701 8625 6D08 0152 C024 O03B2 0202 4ACF 48AF 49CF 
012702 4AA3 4983 49A0 4904 47EA 4820 480B 4802 4890 
012703 48A7 4905 49F1l 4A4A 4903 4960 494F 49A6 4A30 
012704 4935 4910 4A31 4957 4970 49BE 49C4 49C7 4830 
012801 8625 6E€08 0151 CO00 0382 0203 4749 4782 4747 
012802 4732 45€3 45FC 4662 470A 46A6 46E1l 46FA 4749 
012803 4664 4711 4689 45AE 4583 4570 4565 466E 4674 
012804 4604 4778 474D 4602 4682 467D 4613 4643 467C¢ 
012901 D21E 3740400035 0200 9141+04000+00000+16963+00000 
012902+32767+32767 3220+32767=32768 3220+32767=-32768 3220 
012903 0000+16850 1000 0000+000004+00000 0000 0080 0000 
012904 0000+00000 

013001 D21E 374D+00032 0200 91414+01200+00000=-00009+00000 
013002+32767+32767 3220+32767-32768. .3220+32767-32768 3220 
013003 0000416384 1000 0000+00000+00000 0000 0080 0000 
013004 0000+00000 

013101 8625 6808 0132 ,CO0B O3B1l 0239 52¢CC 531B 52D5 
Gaga02e  SOSA.-52CE S206 52€6A,2593€.<52C2 5289 5285 52F6 
Mapeme Bebesngse65 32590 3S26F 5308 5326 5353 5318 52DB 
Meeeee Bef? $247..5342 5360 5286.05290.052Ad.;5296,.>52BD 
013201 C21E 245D+00000 Q000 84C€1+07160+00000+20398+20142 
013202+32767+32767 3220+32767—-32768 3220+32767=-32768 3220 
013203 1078+24200 C€C3CO 0000+24328+16384 0200 0040 00000 
013204 0000+00000 

013301 E21E 235D+00016 3400 94E6€1+10080+00000412279+32507 
013302+32767+32767 3220+32767-32768 3220+32767=-32768 F220 
013303 1578+32767 COCO 00004+32767=32768 0200 0040 0000 
013304 0000+00000 

013401 E21E 1140+00025 3000 96C€0+06800+03000+18263+18311 
013402+32767+32767 3220+32767+00000 3A20+32767+00000 3220 
013403 1A78+26560 COCO 0000419264414336 0001 4C00 000 
013404 0000+00000 

013501 C41E 374D+00041 0000 96C1+10960=-01830+20096+20210 
013502+32767+32767 3220+32767+00000 3A20+32767+00000 3220 
013503 1878+27518 O0CO 0000+24640=-32768 0080 EDO 0000 
013504 8010-00314 
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EXAMPLE OF PROCESS VARIABLE TABLE 


013601 E21E 314D+00022 3C€00 95€1+05500+00000+13411+13107 
013602+32767+32767 3220+32767+00000 3A20+32767+00000 3220 
013603 0151+18080 C4A1l 0000+17815+10560 0002 0070 0000 
013604 0000+00000 

013701 E21E 354D+00019 3C€00 95€1+04920+00000+09283+09151 
013702+32767+32767 3220+32767+00000 3A20+32767+00000 3220 
013703 0152+18944 COA1l 0000+18549+23168 0002 0180 0000 
013704 0000+00000 

013801 C21E 33¢C€+00084 O000 A2C1+20000+00000+07014+06389 
013802+32767+32767 3220+32767+00000 3220+32767+00000 3220 
013803 1478+29175 COCO 00004+29800=—32768 0100 oOvu8sd d000 
013804 0000+00000 

013901 D21E 005D+00009 0400 84€0+4+06100+00000+12899+11393 
013902+32767+32767 32204+32767+00000 3A20+32767+00000 3220 
013903 1378413256 C3CQ 0000+14277+10112 0007 0034 0000 
013904 0000+00000 

014001 8A1E 007€+00313 O000 84€04+06100+000004+15037+00000 
014002432767+32767 3220+32767+00000 3A20+32767+00000 3220 
014003 0141+21700 80A1l 0260400000400000 0000 2200 0000 
014004 0000+00000 

014101 D21E 005D+00006 0400 84C€0+10200+00000+22746+20938 
014102+32767+32767 3220+32767+00000 3A20+32767+00000 3220 
014103 1278421415 ¢3CO 0000+21671+27712 0003 0068 0000 
014104 0000400000 

014201 C21E 32€C+00083 0000 A2C0+20000+00000+06405+06389 
0142024+32767+32767 3220+32767+00000 3220+32767+00000 3220 
014203 1178+25576 COCO 0000+25592=32768 0200 008d 0000 
014204 0000+00000 

014301 821E 366C€4+00304 0000 9100+01200+00000+16384+21845 
0143024+32767+32767 3220+32767=32768 3220+32767=32768 3220 
014303 0140+08704 O0A3 0000+08436=01536 0004 0080 90000 
014304 0000+0000 

014401 821E 376C+00297 O000 9140+04000+00000+16850+16842 
0144024+32767+32767 3220+32767=32768 3220+32767-32768 3220 
014403 0132+19948 COEl FAB1+21251+06400 o000C 0028 0000 
014404 0000+00000 

014501 A21E 336C+00304 3680 91004+01200+00000+16342+16384 
0145024+32767+32767 3220432767"32768 3220+32767-32768 3220 
014503 0169+00022 DOE1l 0000+00000+00000 0000 0044 0000 
014504 0000+00000 

014601 8625 6E£08 0137 C014 0383 0207 2464 2420 240A 
014602 2427 2407 23DE 23FA 240F 2446 2452 2450 2448 
014603 242E 2475 2448 2466 246D 2451 2430 2401 23F9 
014604 2360 23D1l 23C1 2417 2435 2450 2464 2486 2468 
014701 8625 6108 0129 C015 O3B82 0236 4216 4250 4214 
014702 4214 4228 4225 4223 4216 4200 41F8 41F3 4216 
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EXAMPLE OF PROCESS VARIABLE TABLE 


014703 4220 4237 3FE4 4007 4045 4018 4009 4045 402C¢ 
014704 402F 402C 4O0E 4004 4018 4018 403E 3FFF 420F 
014801 8625 6FO8 0136 C013 0383 0208 3461 347E 3524 
014802 358B 3591 350A 343C 3390 3333 342D 3464 3460 
014803 3634 35CD 3488 334€ 342B 3519 3566 34F6 3509 
014804 3536 359A 3544 3494 33F8 32FB8 3258 331F 33FC 
014901 D21E 105D+00012 0600 8441+17440+00000+17465+16910 
014902+32767+32767 3220+32767=32768 3220+32767=32768 3220 
014903 0139412548 COE1l 0000+17543+32192 0002 0240 0000 
014904 0000+00000 

015001 O21E 115D+00038 0800 8441+06460+00000+20028+21015 
015002+32767+32767 3220+32767=32768 3220+32767-32768 3220 
015003 1€78+09396 €3C0 0000+08415=20416 0001 0080 0000 
015004 0000+00000 

015101 D21E 005D+00039 0800 8441+12480+00000+17898+18025 
015102+32767+32767 32204+32767=32768 3220+32767-32768 3220 
015103 1978427521 €380 0000+27575+22080 0002 0030 v000 
015104 0000+00000 

015201 D21E 105D+00040 0800 8440+05340+000004+18839+18841 
015202+32767+32767 32204+32767=32768 3220+32767-32768 3220 
015203 1878+11822 C380 0000+12471+06656 0003 0280 0000 
015204 0000+00000 

015301 E£O0D 6C0C+00065 3200 9700+04000+00000+09440+00000 
015302+32767+32767 3220 

015401 EF16 LBOC+00064 3200 9720+04000+00000+13186+00000 
015402+32767+32767 3220+32767=32768 3220+32767-32768 3220 
015403 0153+00000 F320 

015501 EOOD 40CC+00100 3A00 A200+20000+00000+03641+00000 
015502+32767+32767 3220 

015601 EO0D 41CC+00090 3A00 A200+20000+00000+07597+00000 
015602+32767+32767 3220 

015701 EO0D 42CC+00092 3A00 A200+20000+00000+06266+00000 
015702+32767+32767 3220 

015801 E00D 43CC+00088 3A00 A200+20000+00000+08126+00000 
015802+32767+32767 3220 

015901 EO00D 44CC+00091 3A00 A200+20000+00000+06569+00000 
015902+32767+32767 3220 

016001 E£00D 45¢C+00097 3A00 A200+20000+00000+07585+00000 
016002+32767+32767 3220 

016101 £00D 46CC+00103 3A00 A200+20000+00000+04572+00000 
016102+32767+32767 3220 

016201 FOOD 47CC+00096 3A00 ‘A200+20000+000004+05692+00000 
016202+32767+32767 3220 . 
016301 E00D 48CC+00086 3A00 A200+20000+000004+10305+00000 
016302+32767+32767 3220 

016401 E00D 49CC+00089 3A00 A200+20000+00000+07615+00000 
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EXAMPLE OF PROCESS VARIABLE TABLE 
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Figure D-2 Basic Control Loops 












loop 44 


loop 32 


Product Concentration Control Loops 


Figure D-3 
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Notes on Digital Controllers 


All variable records were written in accordance with the 
Loop Record Builder Form available for this particular 
system, (22,35). Their contents at one particular instant 


suetng the final day of experimentation are shown on Table 
D- e 


Feed Concentration 


Loop 30 was inserted for convenience of operation. When- 
ever the refractometer functioned properly, loop 30 was 
placed in the automatic mode with a setpoint of zero 

and a gain of unity, such that the refractometer reading 
was transmitted unchanged to Loop 43. In cases where the 
refractometer did not function, loop 30 was placed in the 
manual mode, The correct feed concentration was then en- 
tered from the process operator's console into the output 
of this loop. In this way, a reading was available to loop 


43 at all times. 


Differential Vapour Pressure © 


Dierienteies were experienced with the differential vapour 
pressure transmitter. Thus it was necessary to operate loop 
35 in the manual mode for a large part of the program, In 
order to make additional use of this loop, it was connected 
to the first effect pressure transmitter, such that this 


pressure reading was available from the console. 
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Feedback Control 

Loop 44 was the basic feedback controller. Loop 29 was in- 
serted for convenience of operation. It served two main pur- 
poses. Firstly, in case of instrument failure, it could re- 
ceive the correct product concentration from the console, si- 
milar to loop 30. Secondly, when the process model was used 
for control purposes, process and model outputs could be re- 
corded separately by loops 29 and 44, respectively. The out- 
put signal of the feedback concentration controller was 


transmitted to the setpoint of the steam flow controller. 


Feedforward Control 
The feedforward controller Sousiated of loops 45, 68, and 
69. Steady state values of feed flow and concentration were 
supplied from the console as setpoints to loops 68 and 45, 
respectively. The measurements obtained from the process were 
filtered, and the steady state values subtracted from them, 
such that the perturbations, including dynamic compensation, 
were contained in the error signals. Dynamic compensation 
was achieved by the use of exponential filters, the filter 
constants being defined by the time constants in equations 
(5-10) and (5-11). Relation between filter factor FF, time 
constant TI , and poll time interval DT (35): 

In(FF)*TI = -DT 
In the present case TI = 5.3 minutes and DT = 8 seconds. 
It follows that FF = 0.976. 


Gain factors required by equations (5-10), (5-11), and (5-12) 
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D-12 
were changed into range factors and entered into loops 45 
and 68 as proportional constants as follows (see calibra- 


tion curves in Appendix B for recorder ranges). 


Feed concentration: gain required 0.314 
range of steam flow recorder 3.58 


range of concentr. recorder 6.0 

93h 32767 
range factor = ae 
6.0 *32767 


= 0.526 


Feed flow: gain required 0.40 
range of feed flow recorder 8.72 


0.40 
359 *32767 - 


range factor = 0.97 


20 
gi77 *32767 


First effect bottoms flow: gain required 0.61 
range of recorder 6.24 


0.61 
ZB *32767 


range factor = = 1.06 


iO 
$37G, 32767 


The two feedforward controller signals were added by loop 
69 and transmitted to the steam flow controller via the 
output bias of loop 44. Thus the steam flow controller re- 
ceived its steady state signal from the steady state out- 
put of loop 44 and the compensation due to load disturban- 
ces from the bias of loop 44. In cases where feedforward 
control alone was used, the error of loop 44 was artifi- 


cially zeroed out with the aid of loop 29. 

















_ 


: y we 
2 aqool otal bexrsdne bns erotost ogn Bt ant be gases 


>. 93 ow 

; \ - ay x 

‘biso 362) ewollo? as eitneteatos isnao Hrraqon gen 8 db uo, S 
(a9gnax Zebs6o9x tol & x tbneqgh ra aavaus nos : 

= 


[Ss 


: i. j “ a" he i ah 
MfE.Q0 beztvpst oiseg ; soitsx2ngone. » bas" 


rial « 






ore 
r -* ~ a - 7 : a: 
T5O01TO09ST WOLZz m5972 Ito sgnmst as =" a : 
. ‘ = > a 
) IsDtOoSz .TI lo ago0Bt : a 
7 > a 
are - a a 
Plt.U  - tay 
"OL ch cee, 7 he 


\ PI yy , - _ a! " a 2 ; : 
82.0 = ne  TOFSEE SBAeF — ” 
_ - i" we = - _ 


- 
- 


.0 beriupst nitag won't best 
ctl a seb rOuS72 wolt b: 997 ro 93 ares 
az . ¥ a - io 

seeps OPane = 


\ re, ‘ r med 


°.0 = een 10238% ‘egae' 


. 
\ ose 
5 e 


a 


| 
/ 
/ 
“1 
i“>} 
m7 
oo 


- : e > =i 
¢.% webeooes 20 sank es ea bbc 
; ‘ 


fa.0 Bberivoss ates :wolt emodtod_ tne: 239 sonst oy 





nie ay Se A. 6-=- : ae ‘ — & « Lt 
\a \ . a x we 7 a 
yf = FF = tosost ogas 
"2c Coe! 3 9 — 
NONE y.0 i oa “ 
i d e - ; ~~ _ - 2) 
x 
- 1 
qoo!l vd Bebbe stew slangtea - tatios an09 bsiinibatbehs roi 
ai i 
3 rsfifoxjnos wolt masija ody oF ‘bexsiminess & ~ xe } 


cy 


rallovgnd> wort masta sit eurit: |e gor to yaa 6 saat 


“Whee 
te 


' 
+ 
\ 
- 
& 
? 


beede sit moxt langle” ‘stale 







‘ -nedrutelb baol oF sub nokssensqmoo ata 


pfs ‘€. 


brawz0tbhes st oxedw 902 89 ‘at Ma Laide c’ 
eet to rors odd. 


x none (Ps) 
: P Lt 

hese dee 

i 2 _ 





& ed 
oe aw : 








gt 


LIST OF CONTROLLER AND FILTER bahia ntscdh’ 


VARIABLE DESCRIPTION 


STEAM FLOW TO FIRST 
EPFECT 

FIRST EFFECT BOTTOMS 
FLOW 

FIRST EFFECT OVERHEAD 
FLOW 


PRODUCT FLOW 

SECOND EFFECT OVERHEAD 
FLOW 

TOTAL FEED FLOW 
COOLING WATER FLOW TO 
CONDENSER 

FEED WATER FLOW 

TEG FEED SOLUTION FLOW 
FEED CONCENTRATION 
PRODUCT CONCENTRATION 
SECOND EFFECT PRESSURE 
SEPARATOR LIQUID LEVEL 
FIRST EFFECT LIQUID 
LEVEL 

TEG FEED SOLUTION 
TEMPERATURE 


FEED WATER TEMPERATURE 


NAME PROPOR= 


TIONAL 

CON= 

STANT 
Fl 005 
F2 00375 
FS 1.0 
Fé 500 
F7 
F8 405 
F9 005 
Fill 0084 
Fl2 0041 
cl 160 
C6 031 
P22 30460 
Lili 300 
L14 06875 
TT10 1.0 
TT11 1.0 


TABLE D=3 


RESET FiL= 
TIME TER 
SECs © TYPE 
8 
500 UNION 
2000 UNION 
667 UNION 
EXPON 
1000 UNION 
8 EXPON 
333 UNION 
143, UNION 
2000 
684 
2000 EXPON 
4000 EXPON 
4000 EXPON 
32 
16 


FILTER 
CON= 
STANT 


0681 
128 


128 
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Typical Response Obtained with Averaging Control 
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Typical Responses of Feed Temperature T7 to 


a +1.0 1b/min Step Change in Feed Flow Rate 


Figure D-13 
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POLLING SEQUENCE 


TIME IN 
SECONDS 
1 34 36 39 40 41 47 50 51 
2 39 40 41 42 49 51 52 57 
3 33 34 38 39 40 41 45 50 
4 32 39 40 41 49 51 52 59 
5 34 37 39 40 41 50 51 60 
6 39 40 41 43 49 51 52 61 
7 25 29 30 33 34 35 39 40 
8 31 32 39 40 41 49 51 52 
9 34 36 39 40 41 50 51 64 
10 39 40 41 42 49 51 52 65 
11 33 34 38 39 40 41 45 50 
12 32 39 40 41 49 51 52 53 
13 27 34 37 39 40 41 50 51 
14 28 39 40 41 43 46 49 51 
15 25 26 29 30 33 34 35 39 
16 32 39 40 41 49 51 52 55 
17 34 36 39 40 41 50 51 56 
18 39 40 41 42 49 51 52 57 
19 33 34 38 39 40 41 45 50 
20 32 39 40 41 49 51 52 59 
21 34 37 39 40 41 50 51 60 
22 39 40 41 43 49 51 52 61 
23 25 29 30 33 34 35 39 40 
24 32 39 40 41 49 51 52 63 
25 34 36 39 40 41 50 51 64 
26 39 40 41 42 49 51 52 65 
27 33 34 38 39 40 41 45 50 
28 32 39 40 41 49 51 52 67 
29 34 37 39 40 41 50 51 
30 39 40 41 43 49 51 52 
31 25 26 29 30 33 34 35 39 
32 32 39 40 41 49 51 52 55 
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Appendix E 


List of Disk Files 
Program Listings 

BEGIN 

THIS 

FDFWD 

INFER 

DIANE 

EDITH 

COLOR 

THAT 
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LOGICAL 


FILE NOs 


20 
21 
25 


WORD 
NOe 


901-902 
903-904 
905-906 
907-908 
909-910 
9il=912 
913-914 
915-916 
917-918 
919-920 
S2i-S22 
923-924 
925-926 
927-928 
929-930 
921-932 
933-934 
935=936 
937-938 
939=940 
941-942 
943—-944 
945-946 


947 
948 
949 
950-951 
9527953 
954=955 
956 


RIGé 
NAME 


DAFLE 
JOSEF 
BINGO 


LIST OF DISK FILES USED 


(uy SS see sey Mt eR Se my em SES Gy ce one Fe me Son A aio ne ees 


DATA STORED 


DATA ACCUMULATION LOOP CONTENTS 
SORTED VALUES FROM DAFLE FOR PLOTTING 
SUCCESSIVE MODEL OUTPUTS 


LIST OF SKELETON COMMON WORDS USED 


VARIABLE 


NAME 


4 

CST2A 
CST2wW 
CcST1B8 
CSTic 
CST9C 
CST9W 
CST5D 
CST5W 
CST4E 
CST4C 
CSTTF 
CSTTW 


*mAmMoanwyYp 


STONE 
BRICK 
WALL 

ERICA 


IRMA 
LOLA 
JOSEY 


WI 
WST 
INES 


VARIABLE DESCRIPTION 


COMPUTATION INTERVAL IN SECONDS 
PROCESS MODEL CONSTANTS (SEE APPe C) 


INTERMEDIATE QUANTITIES ON FIGURE 504 


DUMMY VARIABLES FOR LATER ADDITIONS 


REAL TIME DIFFERENCE BETWEEN FIRST 
VALUES IN DAFLE AND BINGO IN MINUTES 
FLAG 6 ON FIGURE 5¢5 

SECTOR ADDRESS IN FILE BINGO 

FLAG 5 ON FIGURE 5e¢5 

PRESENT VALUE OF FEED CONCENTRATION 
PRESENT VALUE OF FEED FLOW RATE 

PRESENT VALUE OF STEAM FLOW RATE : 
TIME ELAPSED (NUMBER OF 64 SECe INTER= 
VALS)BEFORE CONTROLLER USES MODEL OUTPUT 
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aD 


101 


110 


112 
1121 


116 
121 
126 


YQOAD 


1265 
127 


128 


MAINLINE BEGIN 


MAINLINE BEGIN 


PROGRAM NAME BEGIN 

PROGRAM TYPE PROCESS PROGRAM» MAINLINE CORE 
LOAD 

ABSTRACT THE PROGRAM SELECTIVELY ACTIVATES 


OR DEACTIVATES DDC LOOP RECORDS 
OR INITIALIZES DATA TRANSFER 
FROM RING BUFFERS TO DISK 


OPERATION PROGRAM IS CONTACTED VIA FUNCTION 

NOs 5 OR KEYBOARD REQUEST 
SUBROUTINES FFINP»s DAINLs RSBUFs OPER»s NONOP 
CALLED VIAQ 


DATA JTYPEsJREAD/291/ 
DEFINE FILE 20(4093209UsKK) 


DECIDE ON WHICH FUNCTION BEGIN IS TO PERFORM 


WRITE(UTYPEs101) 
FORMAT(//TSs'ENTER 1 TO ACTIVATE EVAPORATOR LOOPS's/T5 


Ie*ENTER 2 7 
1 DEACTIVATE EVAPORATOR LOOPS's/T5s"ENTER 3 TO ACTIVATE 
1 OISK FILE" 


CALL FFINP(JREADs190sIsTEROR) 

IF(IEROR) 11291169112 

WRITEC(UTYPEs1121) 

FORMAT(T5s'YOU TYPED IN GARBAGEs TRY AGAIN!) 
GO TO 110 

IF(I) 11291129121 

IF(I=3) 12691269112 

GO TO (126591409150)»5 I 


BEGIN IS TO ACTIVATE LOOPS 
OBTAIN ID OF FIRST AND LAST LOOPS IN SEQUENCE 


WRITE(UTYPE 9127) 
FORMAT(T5s'ENTER ID NUMBERS OF FIRST AND LAST LOOPS IN 


1 SEQUENCE T 
1 BE ACTIVATED') 


CALL FFINP(UREADs2909IDFs0sIDLsTIEROR) 
IF(CIEROR) 12819128591281 
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MAINLINE BEGIN 


1281 WRITE(UTYPEs1121) 
GO TO 128 
1285 CONTINUE 
DO 129 JFIDFsIDL 
CALL OPER(U) 
129 CONTINUE 
WRITE(UTYPE91291) 
1291 FORMAT(T5s'LOOP ACTIVATION COMPLETED!) 
CALL VIAQ 


BEGIN IS TO DEACTIVATE LOOPS 


140 WRITE(UTYPE9141) 
141 FORMAT(T5s'ENTER ID NUMBERS OF FIRST AND LAST LOOPS IN 
1 SEQUENCE T 
1 BE DEACTIVATED!) 
142 CALL FFINP(JREADs290sIDFs0sIDL9IEROR) 
IF(IEROR) 14391449143 
143. WRITE(JTYPE91121) 
GO TO 142 
144 CONTINUE 
DO 145 J=IDFsIDL 
CALL NONOP (J) 
145 CONTINUE | 
WRITE (UTYPE 9146) 
146 FORMAT(T5s'LOOP DEACTIVATION COMPLETED!) 
CALL VIAQ 


BEGIN IS TO ACTIVATE DISK FILE 


150 CALL RSBUF 
CALL DAINL 
WRITE(JTYPE 9151) 
151 FORMAT(T5s'DATA TRANSFER TO DISK HAS BEEN 
1 INITIALIZED!) 
~ CALL VIAQ 
END 
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MAINLINE THIS 


MAINLINE THIS 


PROGRAM NAME THIS 

PROGRAM TYPE PROCESS PROGRAMe MAINLINE CORE 
LOAD 

ABSTRACT THE PROGRAM PRODUCES A GRAPH OF 


VALUES COLLECTED IN THE DISK 
FILE VERSUS TIME 


OPERATION THE PROGRAM IS CONTACTED VIA 
FUNCTION NOs 5 OR KEYBOARD 
REQUEST 


SUBROUTINES CALLED FFINPs GTDDL» GBDATs SCALF 
FCHARs MDGP»s FPLOT»s VIAQ 


DIMENSION-IFR(3)9 ICOM(7)s ITO(3) 

DIMENSION VEC(30)s ITHEAD(20)5 LABX(16)s EDIT(16) 
DIMENSION IFUNC(4)5 I1FCB1(323)5 IFCB2(33)% IVEC(30) 
DIMENSION ABMIN(30)s LABY(32)s ILA(16) 


DATMCTTOSIPLUS/! 'etTOts! bgtet/ | 
DATA IFReICOM/'FR' 9 'OM' sg! Potots! totTtotr stutytety,! 
LI*#¢ 


DATA JTYPEsJREAD s9KPLOT sICHARsNAXISsISORT/2Z91l92909ls9l/ 
DATA XMINsNN/060930/ 

DATA LABX/!TotTtotMtotEty! tot rig tuts! tothlytrtyitnt 
lotU's'TI,gtcEt 
1'S',! ‘1 7 

DATA EDIT/100060 810000 e0slO00000 60 81000000 40516091060 
1910009100060 
10200190601 9s0el91s094%020001/ 

DEFINE FILE 20(4093209UsKK)»9 21(509309UsLL) 
IFUNC(2)=1500 


IFCB1(1)=20 
IFCB1(2)=320 
IFCB1(3)=40 
IFCB2(1)=21 
IFCB2(2)=30 


IFCB2(3)#50 
WRITE (UTYPE91401) 
1401 FORMAT(TSs'ENTER ID OF DATA ACQUISITION LOOP» POLL 
1 TIME INTERVAL 
1N SECONDS'/T5s9'AND DATE. AND TIME AT WHICH THIS LOOP 
1 WAS MADE OPER 
1BLE'/TS59'USING THE FORMAT OF THE FOLLOWING EXAMPLE» 
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EATS. 
an “aes 


dais heeds 


AnD 


OM OO ff 


MAINLINE THIS 


1 0148912858 JA 
1 907914:06') 
145 CALL FFINP(UREAD 96909LPIDsOsLPLT 939RMON1 909IDAT190 
lsTHR1LsOsIMIN1lsI 
1ROR) 
IF(IEROR) 14691479146 
146 WRITE(UTYPE91461) 
1461 FORMAT(T5s'YOU TYPED IN GARBAGEs TRY AGAIN') 
GO TO 145 


CONVERT THE HEXADECIMAL LOOP ID INTO A DECIMAL 
NUMBER 


147 LPID1=LPID/100 
LPID3=LPID-1LOO*LPID1 
LPID2=LP1ID3/10 
LPID3Z=LPID3-10*LPID2 
LPID4G=LPID1*256+LPID2*16+LPID3 


LPID4 CONTAINS THE DATA ACQUISITION LOOP ID IN FORM 
OF A 

DECIMAL NUMBER 

BEGIN DATA TRANSFER FROM DAFLE TO TEMPORARY FILE 
NOe 21 


IFUNC(1)=LPID4 
CALL GBDAT(IFUNC sIFCB1lsIFCB2) 


VERIFY SUCCESS OF DATA TRANSFER 


JINT=IFUNC(3) 
GO TO (19915916917918)s JINT 
is WRITE(UTYPE9151) LPID 
151 FORMAT(T5Ss*LOOP! sT1llesI4eT159'"NOT FOUND IN DISK FILE! 
 L/TSeATE YOUCW 
INT TO TRY AGAINs CALL FUNCTION NOe 5') 
GO TO 28 
16 WRITE(UTYPE9161) IFUNC(2) 
161 FORMAT(TS5sI4e9Tl1s9'VALUES HAVE BEEN FOUND IN DAFLE AND! 
1/TS5e9'HAVE B ; 
1EN TRANSFERRED TO THE TEMPORARY FILE') 
GO TO 20 
i? WRITE(UTYPE9171) 
171 iFORMAT(TSs'THE OISK FILE IS EMPTY AT PRESENT» TRY 
1 LATER!) - 
GO TO 28 
18 WRITE(UTYPEs181) 
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MAINLINE THIS 


FORMAT(T5s'RECORD SIZE OF INPUT FILE IFCB1 IS 

1 NEGATIVE'/T5s'TRANS 

1ER INTERRUPTEDs YOU ARE OFF!) 

GO TO 28 

WRITE(UTYPEs191) IFUNC(2) 

FORMAT(T59'DATA TRANSFER TO TEMPORARY FILE COMPLETED! 
1/T59'NUMBER 

1F POINTS TRANSFERRED I1S'9T37914) 


DATA ARE READY TO BE ACCESSED IN TEMPORARY FILEe 
START PREPARING THE GRAPHe 

CALCULATE DATE AND TIME OF LAST POINT IN TEMPORARY 
FILE e 


ITOSE=LPLT*IFUNC (2) 
ITOMI=ITOSE/60 
ITOHR=ITOMI/60 
ITODA=ITOHR/24 
IREMI=I1TOMI=I TOHR*60 
TREHR=ITOHR=ITODA*24 
IFIHR=IHRI+IREHR 
IFIMISIMINI+IREMI 
IF(60"IFIMI) 20282029203 
IFIMI=IFIMI=60 
IFITHR®IFIHR+1 
IF(24=IFIHR) 20492049205 
IFIHR=IFIHR=24 
ITODA=ITODA+1 


THE TIME AT WHICH THE LAST POINT WAS STORED IN 
DAFLE IN DAYS» 

HOURS» AND MINUTES IS CONTAINED IN ITODAs IFIHR» 
AND IFIMI9 

RESPECTIVELY 


DATA WILL BE MOVED OUT OF THE TEMPORARY FILE AT THE 
RATE OF ONE 

RECORD (30 POINTS) AT A TIMEs 

OBTAIN INPUT PARAMETERS FOR MDGPe 


TOMIN=FLOAT(ITOSE) /6000 

WRITE(UTYPE9206) TOMIN 

FORMAT(TS5s'TOTAL TIME TO APPEAR ON ABSCISSA IS'»sT4l 
lsF7elsT509'MIN 

1TES*) 

WRITE(UTYPEs210) 

FORMAT(T5s'ENTER INFORMATION ON GRAPH IN THIS ORDER! 
1/T59'"MAXIMUM 
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MAINLINE THIS 


LALUE OF DEPENDENT VARIABLE(REAL)'/T59"MINIMUM VALUE OF 
1 DEPENDENT 
LARTABLE (REAL) '/T59'NUMBER OF DIVISIONS ALONG 
1 ABSCISSA( INTEGER) '/T 
l»'NUMBER OF DIVISIONS ALONG ORDINATECINTEGER)!) 
211 CALL FFINP(UREAD 94919YMAXs19YMINSOsNX#O09NYSTEROR) 
IF(TEROR) 21292139212 
212 WRITE(UTYPEs2121) 
2121 FORMAT(TSs'YOU TYPED IN GARBAGE® TRY AGAIN!) 
GO TO 211 
213 WRITE(UTYPEs2131) 
2131 FORMAT(TSs'ENTER LABEL FOR ORDINATE (32 ELEMENTS AT 
1 MOST» ENDING 
LETNAS* _ @ACEFTOUUSTIFIEDIS! 
214 CALL FFINP(UREAD919162sTLA(1) sIEROR) 
IF(IEROR) 21592169215 
215 WRITE(UTYPEs2121) 
GO TO 214 


CONVERT LABEL FROM A2 TO Al FORMAT 


216 DO 218 L=1916 
M=L¥2=1 
K1=ILA(L)/4096 
K2=1LA(L)=K1#4096 
K3=K2/256 
IF (K3)21629216392163 
2162 LABY(M)=K1*40964+(K3=1)*256+64 
GO TO 2165 
2163 LABY(M)=K1*4096+K3%256+64 
2165 M=L*2 | 
K4=ILA(L)-K1#4096-K3%*256 
K5=K4/16 
K6=K4=K5*16 
LABY (M)=K5*4096+K6%256+64 
218 CONTINUE 


DEFINE FIRST 30 ELEMENTS OF TIME VECTOR 


IABSC=0 
ABMIN(1)=060 
DO 22 M=2sNN : 
IABSC=IABSC+LPLT 
ABMIN(M)=FLOAT(IABSC)/6060 
22 CONTINUE 


DEFINE FIRST 30 ELEMENTS OF DEPENDENT VARIABLE 
VECTORs 



























Sys rmaunaene & oa 
DVOJA eMOleIVvid 30 “S0WUAYs TAI UABRI SUBALSAL ter) > 
TAYCRSOSTMIYASSTIOSA LC 
(1 (AZOBTVIVETAMIGAO BVOUA eMOLZIVIO 30 ASSMUAMEL =o 
(AOMFL e Ye OeXeGeWiMve LeXAMYe Lede GABALI ANI F4 JIAD tis 

SISeELSeSLS (AOMBINAL 

oA ESTSeR9YTUIBTI AW fia. 
('ALAQA YAT e30AaRAd at “Q39YT UOY"e@T)TAMAOT ISLS — 
“4° oY ¥Gabe, 226 OORT coe 9 

A TOLSeB9YTLISTIOW €. 

TA @TH3M3j3 $8) 3TAMIGRO HOF! 436A) @3TVa'e aT) TAMAOF LEIS 
VIGNE eTeoM I = 
C*COSTRIT2ZUL TAB © f) HTT 


(O81 4 (L)ALTeSOLeLeQASAL) aHTAA dJa> S4EG 
QLSedlSe@LS (ROMBIIAT 
({SiSeB9YTLISTIAN 21S 
& Las pth Ot "e 
- 1% ys FeAl e ‘eZ 


TAMACA [A OT SA MOSS JIZEGAJ TRAVWOD -- <a 
os. ace? lett femese) oe 






a. MER (+r Se hemegt, aa 
| ate ts ers 09 eos 
S0SeSs b s=i* ae 
“geoen tad AdT=Iy sod 
a - a 
 calsekatSeSaiSien a 
Dd+deS"( LOX) +8QORH LT HEIMIYBAD sors 


als of 08 
seegeseneaeose net) v9Ad €6! 


¥ S3W 
a. saraesnarseeoaeens(ivBAd 
20 | SEs ts umacens ee 
singe WE Mie PN 
AOTISV 3aMIT 30° e1aM3.43 of Tenia. ) 
ut See TES ATS 


i eer es 7 7 
eigen 














se oo ROT 
ro. eee E- 
: 2 ne? > 5 ae = 90a oT yer? Pad 
| | Oe0a\(22galie 
an b a aie 
tb Ss aaeatnan TAC 
. one oo Firs fae : 
Rees 


Se 


7 : oe Dee x 

¥. a 
See eee (>) 
en te ate 


KOE 


HANS 


MAINLINE THIS 


CONVERT DATA IN TEMPORARY DISK FILE TO ENGINEERING 
UNITSe 


CALL GTDDL(LPID4sIDSPYsICONAsICONBs ITEDITsIUNITSIER) 

LRUTERV522392229223 
222 WRITE(UTYPEs2221) 
2221 FORMAT(T5s'DONOR LOOP COULD NOT BE FOUND» DATA 

1 TRANSFER INTERRUPT 

1D'/TSes'IF YOU LIKE TO START AGAIN» CALL FUNCTION NOc 

1 CEN} 

GO T6728 
223 CONB=FLOAT(ICONB) *EDIT(IEDIT) 

CONA=FLOAT(ICONA) *EDIT(IEDIT) /6553640 

K=1 

READ(21'K) (IVEC(L) sL=lsNN) 

DO 24 MM=19NN 

VEC (MM) =CONA*FLOAT(IVEC(MM) )+CONB 
24 CONTINUE 


START GRAPHs PLOTTING FIRST 30 POINTSe 


WRITE(UTYPE 9245) 
245 FORMAT(TSs'ENTER TITLE OF PLOT AS A VECTOR OF 40 

1 CHARACTERS!) 
2452 CALL FFINP(UREAD919202sITHEADs ITEROR) 

IF(IEROR) 24692479246 
246 WRITE(UTYPEs2121) 

GO TO 2452 
247 CALL SCALF(1e09160900e09000) 

CALL FCHAR(30587029001590029080) 

WRITE (79248) IHEAD 
248 FORMAT(20A2) 

CALL FCHAR(32e95970e090019006139000) 

WRITE(79249) IFR* RMON1s IDAT1l»s ICOMs IHR1s IMIN1s ITO 

ls RMON1» ID 

~1T1ls IPLUS»s ITODAs ICOMs IFIHR»s IFIMI 

249 FORMAT(3ZA2sA4sI297ALsSI291393A2% sA4sIT29AlsI29/7Als12913) 

CALL FPLOT(4+590009060) 

CALL MDGP(ABMINsVECeNNo» TOMINeXMIN 9 YMAX 9 YMINSNXONY 

lsKPLOT »s ICHAR s NAX 

1SsISORT sLABX sLABY) 


PLOT THE REMAINDER OF THE DATA IN THE TEMPORARY 


FILE ON THE 
SAME GRAPHe 


262 JRE=IFUNC(2)=K*NN 
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263 
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267 


268 


270 
28 


MAINLINE THIS 


IF(URE=NN) 26392649264 
NNFJRE 

DO 265 KL=lo¢NN 
IABSC=IABSC+LPLT 
ABMIN(KL) =FLOAT( IABSC) /600 
CONTINUE 

K=K+1 

READ(21'K) (CIVEC(L) sL=l9NN) 
DO 266 KM=1sNN 

VEC(KM) =CONA*FLOAT(IVEC(KM) 
CONTINUE 

DO 267 KP=loeNN 


)+CONB 


CALL FPLOT(OsABMIN(KP) sVEC(KP)) 


CONTINUE 

IF(URE=NN) 26892689262 
XNEW=TOMIN+4e¢0*( TOMIN@=XMIN) 
YNEW=YMIN@=205%* (YMAX=YMIN) /5 
CALL FPLOT(+5 9XNEWs YNEW) 
WRITE(JTYPE 9270) 
FORMAT(TSs'LIKE IT OR NOT» 
CONTINUE 

CALL EXIT 

END 
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MAINLINE FOFWD 


MAINLINE FDFWD 


PROGRAM NAME 


PROGRAM TYPE 


ABSTRACT 


OPERATION 


SUBROUTINES CALLED 


CALL JACK (UUDY) 


GO TO (293)>» 
CALL QUEEN 
CALL KING 
CALL VIAQ 
END 


JUDY 


FDFWD 


PROCESS PROGRAMs MAINLINE 
CORE LOAD 


THIS PROGRAM PERFORMS 
MATERIAL BALANCES AROUND 
VARTOUS SECTIONS OF THE 
EVAPORATOR BASED ON FLOW AND 
CONCENTRATION READINGS 
OBTAINED FROM THE 

PROCESS VARIABLE TABLE 


PROGRAM IS CONTACTED VIA 
FUNCTION NOw 5 OR 
KEYBOARD REQUEST 


JACK s QUEENs KINGs GTVLU 
VIAQ 


eS 

















dwaGs 


JVIIVIAM «eMAROOARG S29D0HN9 re _39YT } 
QA0J 3703 © 14 
j aes Le ‘aBees 
2MAOAASG MARDORS SIHT 
- QHUORA SADNVAJAB BATASTAM’ St 74 Diener sen 
BHT 30 2nHO1TDBS SUOTAAV 
GAA WOJ9 WO G3e@AG ROTAROSAV] segue: Tai: 
22vVIGASR WOLTAATHIDHOD 0%) |SKinMHRs +) TGUNE RRR 
BHT MORA GSWIATBO © . Rape rete 
B4GAT 4IBATRAV 22320h9 | tie S55 eee Laie BG EBD: 
‘ype bagpeagign ae i Sk Eee 
ATV Q3TIATNOD él “MAROORG. HULA“ Die 
- AO @-sOK KVOLTSIAUD tv oe) +i Reteste ak FAR 
a - TER Sabet ~ Dept tries. 
: is Tro Avil Ak ee : 
UJVTO ONES pneave exDAL : O3JJAD 23¥ITL 
CAL 









DQNAOAN 


7 


g 
91 


he 


e= AizZ 
SUBROUTINE JACK 
SUBROUTINE JACK(UUDY) 


THIS SUBROUTINE TAKES READINGS REQUIRED FOR 
MATERIAL BALANCE 
FROM PROCESS VARIABLE TABLE 


COMMON FleF29F5sF6eF79F8sCI9CO 
JUDY=1 

CALL GTVLU(306915FlsLULUs4) 

GO TO (2999999)s LULU 

CALL GTVLU(33791sF2sLULUs4) 

GO TO (33999999)s LULU 

CALL GTVLU(336915F5sLULU94) 

GO TO (4899999)»s LULU 

CALL GTVLU(33891 sF6sLULU94) 

GO TO (5999999)s5 LULU 

CALL GTVLU( 29491 9F7sLULUs4) 

GO TO (55999999) LULU 

CALL GTVLU(32991sF8sLULUs4) 

GO TO (6999999)s LULU 

CALL GTVLU(32391s9CI sLULUsl ) 

GO TO (73999999)s LULU 

CALL GTVLU (324s1lsCOsLULUs1) 
GO TO (95999989) LULU 
WRITE(2991) LULU 
FORMAT(T5s'GTVLU ERROR CODE'sI39T25s'"EXECUTION 


1 INTERRUPTED!) 


JUDY=2 
RETURN 
END 


“i 33 
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SUBROUTINE QUEEN 
SUBROUTINE QUEEN 


THIS SUBROUTINE CALCULATES ERRORS OF CLOSURE FOR 
PRESENT 
STEADY STATE CONDITIONS 


COMMON FlsF2sF5sF6sF79F8sCI sCOsSTECsE1lsE29E3 9E49E5 
STECS(F8"F8*CI/CO) /F1 

El=(F8—FS=F7=F6) *100e0/F8 

E2=(F8—-F5=F2)#1000e0/F8 

E3=(F2=-F6=F7)*1000e0/F2 

E4=(FE*CI=FE#CO) *1000e0/(F8*CI ) 

E5=(F8*(100e0=CI )—F6*(100e0—CO)MFI*XSTEC#1LOOCO) /(F1 


1*¥STEC) 


RETURN 
END 


2 
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SUBROUTINE KING 
SUBROUTINE KING 


THIS SUBROUTINE WRITES ALL ERRORS OF CLOSURE ON THE 
TYPEWRITER 


DIMENSION FAYE(8) 

COMMON FAYE sSTEC 9E1sE29E3sE49E5 
WRITE(2910) STEC 
FORMAT(//T5s'PRESENT STEAM ECONOMY IS' 5F8e3) 
WRITE(2911) 

FORMAT (/T5s'PRESENT ERRORS OF CLOSURE IN PER CENT 
1 ARE?) 

WRITE(2921) El 

FORMAT(T59'TOTAL PROCESS! 9F8el) 
WRITE(2931) E2 

FORMAT(T5e9'FIRST EFFECT! sF9el) 
WRITE(2941) E3 

FORMAT(T5e'SECOND EFFECT! »F8el) 
WRITE(2951) E4 
FORMAT(TSs'SOLUTE'»T18sF8el) 
WRITE(2961) E5 

FORMAT(T5s "OVERHEAD FLOWS! 9F7e1l) 
WRITE(2971) 

FORMAT(TSs'END!//) 

RETURN 

END 
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MAINLINE INFER 


MAINLINE INFER 


PROGRAM NAME 


PROGRAM TYPE 


ABSTRACT 


OPERATION 


SUBROUTINES CALLED 


EXTERNAL EDITH»s DIANE 
CALL ZERO( JANE) 


GO TO (10920930940)s JANE 


CALL CHAIN(DIANE) 
CALL CHAIN(EDITH) 
CALL CANCL (4) 
CALL VIAQ 

END 


INFER 


PROCESS PROGRAM» MAINLINE 
CORE LOAD 


PROGRAM SUPERVISES 
INFERENTIAL CONTROL BY 
PROVIDING BRANCHES TO VARIOUS 
CORE LOAD 


PROGRAM IS CONTACTED VIA 
FUNCTION NOe 5 OR 
KEYBOARD REQUEST 


ZEROs VIAQs CHAINs CANCL 
FFINP 
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SUBROUTINE ZERO 


SUBROUTINE ZERO( JANE) 
WRITE (2950) 

50 FORMAT(TS5s'THESE CHOICES ARE OPEN TO YOU NOW!/T59!' 1 
1 START COLLE 
1TING INFORMATION FOR INFERENTIAL CONTROLLER'/T5s!' 2 
1 GET TIMER NO 
24 ON ITS WAY'/T59!' 3 CEASE INFERENTIAL CONTROL'/T5Ss'! 
a4 CALL IT 
3A DAY') 

60 CALL FFINP( 19190 9JANEsTRENE) 

IF (IRENE) 70980970 

70 WRITE(2971) 

71 FORMAT(T5s'YOU TYPED IN GARBAGE» TRY AGAIN!) 
GO TO 60 

80 IF( JANE) 70970890 

90 IF(UANE=4) 1209100970 

100 WRITE(2:9101) 

101 FORMAT(T5s'OKAYs YOU ARE OFF!) 

120 RETURN 
END 
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MAINLINE DIANE 


MAINLINE DIANE 


PROGRAM NAME DIANE 

PROGRAM TYPE PROCESS PROGRAMs MAINLINE 
CORE LOAD 

ABSTRACT PROGRAM PROVIDES ALL CON= 


STANTS AND PARAMETERS 
FOR INFERENTIAL CONTROLLER 


OPERATION PROGRAM IS CONTACTED VIA CALL 
CHAIN FROM INFER 


SUBROUTINES CALLED ONEs TWOs THREE» FOURs SIX 
CHAINs FFINPs VIAQs GTVLU 


EXTERNAL EDITH 

DIMENSION GAIL(450) 

COMMON/INSKEL/GAILeHsCST2A sCST2WsCST1BsCST1ICsCSTIC 
l»sCST9IWsCST5D 
COMMON/INSKEL/CST5WsCSTSE sCST4C sCSTTF sCSTTWsAsBsCeDoF& 

loF 
COMMON/INSKEL/STONE sBRICK sWALL 
COMMON HOLD] sHOLD2ZsSTECsCI sCOsw! sWSTsC2sFLOOR SKY sDUST 
10 CALL ONE , 
CALL TWO(UULITA) 
GO TO (20990)s JULIA 
20 CALL THREE( JEAN) 
GO TO (109909903s30)s JEAN 
30 CALL FOUR 
CALL SIX 
CALL CHAIN(EDITH) 
90 WRITE(29901) 
901 FORMAT(T5Ss'OKAYs YOU ARE OFF!) 
CALL VIAQ 
END 
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SUBROUTINE ONE 


SUBROUTINE ONE 


THIS SUBROUTINE REQUESTS CONSTANTS OR PARAMETERS 
FROM POC 


DIMENSION GAIL(450) 

COMMON/INSKEL/GAIL#H 

COMMON HOLD1 sHOLD2sSTEC 

WRITE(2910) 

FORMAT(T5s*ENTER AS REAL NUMBERS'/T59'HOLD=UP IN FIRST 
L> EPFECT (an 
18S'/T5s'HOLD=UP IN SECOND EFFECT IN LBS'/T59'STEAM 
1 ECONOMY! /T59'C 
2MPUTATION INTERVAL IN SECONDS') 

CALL FFINP (19491 sHOLD1s1l1sHOLD2ZslsSTEC91 sHsKATE) 

IF (KATE) 23930923 

WRITE(29231) 

FORMAT(T5s'YOU TYPED IN GARBAGEs TRY AGAIN!) 

GO TO 20 

RETURN 

END 
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SUBROUTINE TWO 


SUBROUTINE TWO(JULTA) 


THIS SUBROUTINE OBTAINS PRESENT VALUES OF PROCESS 
VARTABLES 
FROM LOOP RECORD TABLE 


COMMON HOLD1sHOLD2sSTEC CI sCOow!l swST 
JULTA=1 


OBTAIN FEED CONCENTRATION CI 


10 CALL GTVLU(32391sCIsLAURAs1) 
GO TO (20980980980) s LAURA 


OBTAIN PRODUCT CONCENTRATION CO 


20 CALL GTVLU(32491sCOsLAURAS1) 
GO TO (30%80s80980)s LAURA 


OBTAIN FEED FLOW RATE wi 


30 CALL GTVLU(32991sWI sLAURAs4) 
GO TO (40980380980) s LAURA 


OBTAIN STEAM FLOW RATE WST 


40 CALL GTVLU(30691 9WSTsLAURA 94) 
GO TO (90980980980)s LAURA 
80 WRITE(2981) LAURA 
81 FORMAT(T5s'GTVLU ERROR CODE's13/T5s'ENTER 1 TO REPEAT 
1 GTVLU OR 2 
L@OGIVES YP! ) 
82 CALL FFINP(1lslsOsJULIAsKATE) 
IFCKATE) 859855985 
85 WRITE(29851) 
851 FORMAT(T5s'YOU TYPED IN GARBAGEs TRY AGAIN!) 
GO TO 82 / 
855 IF(JULIA) 85985586 
86 IF(JULIA=2) 87987985 
87 GO TO (10990)s JULIA 
90 RETURN 
END 
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SUBROUTINE THREE 


SUBROUTINE THREE (JEAN) 


THIS SUBROUTINE VERTFIES THE MATERIAL BALANCE 


COMMON HOLD1 sHOLD2sSTECsCI sCOsw! sWSTsC2 
VAPOR=WST*STEC 

C2=CI*WI/(WI=VAPOR/260) 

COCAL=CI*WI/(WI=VAPOR) 

WRITE(2920) CIs COs WIs WST 

FORMAT(TSs "READINGS ARE'sT199'CI=!'sT22sF7Te29T31l9'CO=' 
loT34sF7e2/Tl 
lo'Wl='sT229F7e2sT31le'WST=! sT3259F 602) 

WRITE(2930) C29 COCAL 

FORMAT(TSs*CALCULATION YIELDS! 9T249'C2='sT27sF7e29T37 
ls 'COCAL=!59T4 
leF7e2) 

WRITE(2940) 

FORMAT(TSs*YOU HAVE THE CHOICE TO'/T59!'!1 HAVE READINGS 
1 REPEATED! / 

153*2° ADJUST “READINGS! /TS59"*3 CALLOIT OFF'/T559'4 

1 CONTINUE?) 

CALL FFINP( 19190 sJEANsKATE) 

IF (KATE) 43944943 

WRITE(29431) 

FORMAT(T5s'YOU TYPED IN GARBAGEs TRY AGAIN!) 

GO TO 42 

IF(JEAN) 43943945 

IF(UEAN@=4) 50960943 

GO TO (60952960960) s JEAN 

WRITE(29521) 

FORMAT(TSs ENTER ADUYJUSTED READINGS CI» COs WI» wWST') 
CALL FFINP(1le4elsCIslsCOslsWI slsWSTsKATE) 

IF(KATE) 54910954 

WRITE(29431) 

GO TO 53 

RETURN 

END 
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SUBROUTINE FOUR 


SUBROUTINE FOUR 


THIS SUBROUTINE CALCULATES CONSTANTS AND 
INTERMEDIATE QUANTITIES 


DIMENSION GAIL(450) 
COMMON/INSKEL/GAIL sH»CST2A 9CST2W 9 CST1IB 9CSTIC9CSTIC 
LsCSTOWsCST5D 
COMMON/INSKEL/CST5W sCST4E sCST4C sCSTTF sCSTTW8AsB9CoDok 
lof ; 
COMMON HOLD1sHOLD2sSTEC9CI »COoWI sWST C2 » 
CST2A= (STECHWST=200#W1)/(260*HOLD1) 
CST2W=STEC#C2/(2e0*HOLD1) 
CSTIB=(STEC#WST=200#W1)/(200*HOLD1) 

CSTIC=WI/HOLD1 

CSTOC=(STEC#WST=2e0*W1)/(260*HOLD1) 

CST9W=(CIl=C2) /HOLD1 

CST5D=(STEC#WST=<WI) /HOLD2 

CSTSW= (200*STEC¥CO=STEC#C2)/(260*HOLD2) 
CST4E=(STEC#WST=W1) /HOLD2 

CST4C= (200*WI=STEC#WST) /(2e0*HOLD2) 
CSTTF=(STEC#WST=WI) /HOLD2 

CSTTW=(C2=CO) /HOLD2 

A==<wST*CST2W/CST2A 

B==CI¥CSTIC/CST1B 

C==w1¥CSTOW/CSTIC 

E==-C2*CST4C/CST4E 

Fe-wl*CSTTW/CSTTF 

D==-WST*CST5W/CST5D 

RETURN 


SUBROUTINE SIX 


SPARE SUBROUTINE FOR EVENTUAL ADDITIONS TO PROGRAM 


DIMENSION GAIL(450) sHELEN(19) sOLGA(8) 
COMMON/INSKEL/GAIL»sHELEN»s STONE sBRICKs WALL 
COMMON OLGAsFLOOR»s SKY »s DUST 

RETURN 

END 
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MAINLINE EDITH 


MATNLINE EDITH 


PROGRAM NAME EDITH 


PROGRAM TYPE PROCESS PROGRAMs 
MAINLINE CORE LOAD 


ABSTRACT PROGRAM INITIALIZES 
INFERENTIAL CONTROLLER 


OPERATION PROGRAM IS CONTACTED VIA CALL 
CHAIN FROM DIANE 


SUBROUTINES CALLED FIVEs REPETs QUEUEs VIAQ 
FFINPs DYTIMs SEVEN» INCHR 


EXTERNAL COLORs MVCONs POKUS 
DIMENSION VIOLA(472) 
COMMON/INSKEL/VIOLAs ERICAs IRMAs LOLA» JOSEYs CIs WIs 


1 WSTs INES 


CALL SEVEN(LISE»sMOLLY) 

GO TO (10920930)» LISE 

CALL FSVE 

CALL REPET 

CALL QUEUVEI(COLORs190945160) 
GO TO 40 ‘ 

CALL REPET 

CALL QUEUE(MVCON919094sMOLLY) 
GO TO 40 

CALL REPET 

CALL QUEVE(POKUSs50s094sMOLLY) 
CALL VIAQ 

END 
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SUBROUTINE SEVEN 
SUBROUTINE SEVEN(LISE»MOLLY) 


THIS SUBROUTINE SETS FLAGS LOLA AND JOSEY AND 
SELECTS ~PHE 
CORE LOAD TO BE CALLED BY TIMER 4 


DIMENSION ANNA(473) 

COMMON/INSKEL/ ANNAs IRMAs LOLAs JOSEYs Cl» Wl» WST» 
LF DNES 

LOLA=1 

JOSEY=1 

MOLLY=1500 

WRITE(2912) 

FORMAT(TSs'ENTER NUMBER OF POINTS BY WHICH MODEL 
1 SHOULD LAG PROCE 
1S) 

CALL FFINP(1lelsOsINESsKATE) 

IF(KATE) 15919915 

WRITE(2935) 

GO TO 14 

WRITE(2920) 

FORMAT(TSs'SELECT TIMER CALLED CORE LOAD'/T5,5!' 1 
1 COLOR'sT159' 2 M 
ICONEST25e*%13 POKUS! ) 

CALL INCHR(1lsLISEs193) 

GO TO (60930%930)s LISE 

WRITE(2931) 

FORMAT(T5s'ENTER NUMBER OF TIMER CYCLES(INTEGER) ') 
CALL FFINP( 19190 %MOLLY SKATE) 

IF(KATE) 33960933 

WRITE (2935) 

FORMAT(T5s'YOU TYPED IN GARBAGE®s TRY AGAIN!) 

GO TO 32 

RETURN 

END 
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SUBROUTINE FIVE 
SUBROUTINE FIVE 


SUBROUTINE DETERMINES DUTIES OF TIMER CALLED CORE 
LOAD COLOR 


DIMENSION VIOLA(472) 

COMMON/INSKEL/VIOLA®s ERICA®s IRMA 

WRITE (29101) 

FORMAT(T5s'SHOULD THE INFERENTIAL CONTROLLER! /T59! 1 
1 RECORD ITS O 
1TPUT ONLY'/T59!' 2 ACTUALLY CONTROL STEAM FLOW!) 
CALL FFINP(1l9lsO0sTRMASKATE) 

IF(KATE) 16917916 

WRITE (29161) 

FORMAT(T5s'YOU TYPED IN GARBAGEs TRY AGAIN!) 

GO TO 15 

IFCIRMA) 16916918 

IF(IRMA=2) 23923916 

CALL DYTIM(IDAsIRIS) 

WRITE(29232) IDAs IRIS 

FORMAT(T59'THE TIME IS NOW's15916) 

WRITE (29234) 

FORMAT(TSs*ENTER TIME ELAPSED SINCE ACTIVATION OF DISK 
1 FILE IN MI 
LUTES(REAL)') 

CALL FFINP(1lslelsERICAsKATE) 

IF(KATE) 24592479245 

WRITE(29161) | 

GO TO 24 

GOTO (25926)s5 IRMA 

WRITE (29251) 

FORMAT(TSs'THE INFERENTIAL CONTROLLER WILL BE 

1 ACTIVATED ONE MINUT 

1 FROM NOW AND SEND ITS OUTPUT TO DISK FILE BINGO!) 
GO TO 27 

WRITE(29261) 

FORMAT(T59'THE INFERENTIAL CONTROLLER WILL BE 

1 ACTIVATED ONE MINUT 

1 FROM NOW AND SEND ITS OUTPUT'/T5s'TO MEASUREMENT OF 
1 LOOP 0144 AS 

2WELL AS TO DISK FILE BINGO') 

RETURN / 

END 
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MAINLINE COLOR 


MAINLINE COLOR 


PROGRAM NAME COLOR 

PROGRAM TYPE PROCESS PROGRAMs 
MAINLINE CORE LOAD 

ABSTRACT THIS IS THE INFERENTIAL CON@= 
TROLLERe IT READS PROCESS 
INPUTS AND CALCULATES PRO= 
CESS OUTPUTS AT SPECIFIED 
INTERVALS. 

OPERATION PROGRAM IS QUEUVED BY TIMER 4 


SUBROUTINES CALLED GREEN» BLUEs WHITEs GTVLU 
PTVLU 


DIMENSION GAIL(450)s HELEN(4) 
COMMON/INSKEL/ GAIL» Hs CST2As CST2Ws CST1Bs CST1C»s 
1 CST9Cs CSTIW 
COMMON/INSKEL/ CST5D9 CST5Ws CST4Es CST4Cs CSTTF es 
1 CSTTWs As Be C 
COMMON/INSKEL/ Ds Es Fs HELENs IRMAs LOLAs JOSEYs CI» 
1 WIs WST 
COMMON/INSKEL/ INES 
COMMON CO 
DEFINE FILE 25(480929UsLOLA) 
GO TO (14915)s JOSEY 
14 CALL GREEN 
GO TO 30 
15 CALL BLUE 
CALL WHITE 
30 CALL VIAQ 
* END 
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SUBROUTINE BLUE 
SUBROUTINE BLUE 


THIS SUBROUTINE WILL USE PROCESS MODEL TO CALCULATE 
PRODUCT 
CONCENTRATION FOR PRESENT POLL TIME INTERVAL 


DIMENSION GAIL(450)s HELEN(4)s MARIE(3) 
COMMON/INSKEL/ GAILs Hs CST2As9 CST2Ws CST1Bs CST1Cs 
LCESTIGs VCSTSA 

COMMON/INSKEL/ CST5Ds CST5Ws CST4Es CST4Cs CSTTF»s 
1 CSTTWs As B 

COMMON/INSKEL/Z Cs De Es Fs HELENs MARIEs CIs WI» WST 
COMMON CO 

XK=CST1C#CI 

XA=CST1B 

POWER=XA#*H/6000 
SECND=XK#(1LleO-EXP( POWER) )/XA 
B=R#EXP (POWER) =SECND 
XK=CST2W*WST 

XA=CST2A 

POWER=XA*H/6000 
SECND=XK*(1e¢0=EXP (POWER) )/XA 
A=A*EXP (POWER) =SECND 
XK=CSTIW*WI 

XA=CST9C 

POWER=XA#H/60 00 
SECND=XK*(1e0=EXP (POWER) )/XA 
C=C#EXP (POWER) =SECND | 
C2=A+B+C 

XK=CST4C#C2 

XA=CSTGE 

POWER=XA#H/6000 
SECND=XK*#(1L1eO=EXP (POWER) )/XA 
E=E*EXP (POWER) =SECND 
XK=CSTTWeWwl 

XA=CSTTF 

POWER=XA#H/60e0 
SECND=XK*(LeO=EXP(POWER) )/XA 
F=F#EXP (POWER) =SECND 
XK=CSTSW*WST 

XA=CST5D 

POWER=XA#H/60 020 j 
SECND=XK#(1leO=EXP (POWER) )/XA 
D=D#EXP (POWER) =SECND 
CO#E+F+D 

RETURN 

END 
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SUBROUTINE GREEN 
SUBROUTINE GREEN 


THIS SUBROUTINE READS PRESENT VALUES OF PROCESS 
VARTABLES FROM 

LOOP RECORD TABLE AND STORES THEM IN INSKEL COMMON 
FOR LATER USE 


DIMENSION ANNA(473)9 JOYCE(2) 

COMMON/INSKELZ ANNA JOYCEs JOSEYs CIs WI» WST 
CALL GTVLU(323s1lsCIsLAURA9s1) 

GO TO (298%9898)s LAURA 

CALL GTVLU(329591leWI sLAURA 94) 

GO TO (3989898)s LAURA 

CALL GTVLU(30691sWST »sLAURA 94) 

GO TO (998%98%8)s LAURA 

WRITE(2981) LAURA 

FORMAT(TSs'GTVLU ERROR CODE'sI3/TSs'PREVIOUS LOADS 


1 WIbe BE USED*) 


JOSEY=2 
RETURN. 


SUBROUTINE WHITE 
THIS SUBROUTINE WRITES CALCULATED VALUES OF PRODUCT 
CONCENTRATION INTO A DISK FILE AND TRANSFERS THEM 
TO THE DDC LOOP RECORD 


DIMENSION ANNA(473) 
COMMON/INSKEL/ ANNA®s IRMAs LOLA® JOSEYs CIs WIs WST» 


1 INES 


COMMON CO 

WRITE(25'LOLA) CO 

GO TO (392)9 IRMA 

IF( INES) 239239201 

LINDA=LOLA@INES 

IF(LINDA) 393921 

READ(25'LINDA) CODL 

OUTP=CODL*5e0 

LOLA=LINDA+INES 

GO TO 24 

OUTP=CO*5 20 

CALL PTVLU(29793 sOUTP»MARY) 

GO TO (3922922922922s22922822922)s MARY 
WRITE(29221) MARY 

FORMAT(T5s'PTVLU ERROR CODE'sI139T269'NO PUTTING THIS 


LOT IME*) 


JOSEY=1 
RETURN 
END 
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MATNLINE THAT 


MAINLINE THAT 


PROGRAM NAME 


PROGRAM TYPE 


ABSTRACT 


OPERATION 


SUBROUTINES CALLED 


DIMENSION VIOLA(472) 


THAT 


PROCESS PROGRAM» 
MAINLINE CORE LOAD 


THIS PROGRAM PRODUCES A POINT 
PLOT OF THE PRODUCT CONCEN= 
TRATION CALCULATED BY THE 
INFERENTIAL CONTROLLER VERSUS 
TIMEe THE PLOT OVERLAYS THE 
LINE PLOT OF THE MEASURED 
CONCENTRATION 


THE PROGRAM IS CONTACTED VIA 
FUNCTION NOe 5 OR 
KEYBOARD REQUEST 


HAZEL» RENEEs FFINPs SCALF 
FPLOT»s POINT» VIAQ 


COMMON/INSKEL/ VIOLAs ERICAs IRMA» LOLA 


COMMON XOLDs YOLDs TOMIN 


DEFINE FILE 25(48092sUsLOLA) 


CALL HAZEL 
CALL RENEE 
CALL VIAQ 
END 
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SUBROUTINE HAZEL 
SUBROUTINE HAZEL 


THIS SUBROUTINE DETERMINES THE PRESENT PLOTTER PEN 
POSITION AND 
MOVES THE PEN TO THE ORIGIN OF THE PROPOSED PLOT 


COMMON XOLDs YOLDs TOMIN 
XMIN=0«0 
WRITE (2941) 
41 FORMAT(TSs'ENTER INFORMATION FOR SUPERIMPOSED PLOT'!/T5 
ls'TOTAL MIN 
1TESs YMAXs YMIN (ALL REAL)!') 
42 CALL FFINP(193s1lsTOMINs ls YMAXs1lsYMINSKATE) 
IF(KATE) 43944943 
43 WRITE (29431) 
431 FORMAT(TSs'YOU TYPED IN GARBAGE» TRY AGAIN!) 
GO TO 42 
44 XOLD=TOMIN+4e0*( TOMIN@=XMIN)/8e0 
YOLD=YMIN@=2065%#(YMAX=YMIN) /560 
XSCAL=8e0/(TOMIN™@=XMIN) 
YSCAL=5e0/(YMAX=YMIN) 
CALL SCALFI(XSCAL sYSCAL sXOLDsYOLD) 
CALL FPLOT(+39XMINsYMIN) 
RETURN 


SUBROUTINE RENEE 


THIS SUBROUTINE PLOTS THE PRODUCT CONCENTRATION 
CALCULATED BY THE INFERENTIAL CONTROLLER OVER THE 
PLOT OBTAINED FROM THE RING BUFFER 


DIMENSION VIOLA(472) 
COMMON/INSKEL/ VIOLAs ERICAs IRMA» LOLA 
COMMON XOLDs YOLDs TOMIN 
X=ERICA=240 
XINCR=3 00#64e0/6000 
NULLY=LOLA 
DO 59 JOAN=lsNULLY 93 
X=X+XINCR 
READ(25'JOAN) Y 
CALL FPLOT(+39X9Y) 
CALL FPLOT(=29XsY) 
CALL POINT(O) 
IF (X=TOMIN) 59962962 

59 CONTINUE 

62 CALL FPLOT(+3sXOLDsYOLD) 
RETURN 
END 
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